The method of Virtual Elevation and all its variants are derived from Newton’s 2nd Law for linear motion, neglecting rotational inertia terms. Do those additional terms play a significant role in elevation-matching? Can we accurately measure the moment of inertia (MOI)? In this post, I’ll describe the terms and provide a method for testing MOI.
[edit: Thanks to Nicko on the SlowTwitch forum for pointing out that the added mass terms in the equations are I/R^2 instead of I/R]
Slope: The Critical Calculation
The method of Virtual Elevation is a methodology for determining the Crr and CdA of a bike-rider combination by matching a derived elevation profile against some features of a known elevation profile. Under ideal conditions involving little wind and steady rider position the method works well. But the subtlety of the calculation lies in determining the slope angle, theta, from the linear equations of motion:
F - Crr m g cos(theta) - 1/2 CdA rho v^2 - m g sin(theta) = m a
Slope angle can be computed in closed form by making a small-angle approximation:
sin(theta) = tan(theta) = theta
cos(theta) = 1
The resulting calculation for slope angle is:
theta = ( F - Crr m g - 1/2 CdA rho v^2 - m a ) / ( m g ) .
An iterative method for slope angle calculation can also be used, however it is questionable whether this is necessary for most elevation profiles.
Newton’s Second Law With Rotational Inertia
The above equation of motion neglects the rotational inertia of the front and rear wheels. Whether this is an issue is for others to decide. My aim here is to add a correction term to the existing slope calculation to incorporate rotational inertia terms. Newton’s Second Law with rotational inertia is:
F - Crr m g cos(theta) - 1/2 CdA rho v^2 - m g sin(theta) = ( m + I1/R^2 + I2/R^2 ) a
I1 and I2 are the moments of inertia (MOI) of the front and rear wheels. This modifies the small-angle slope angle formula as follows:
theta = ( F - Crr m g - 1/2 CdA rho v^2 - ( m + I1/R^2 + I2/R^2 ) a ) / ( m g ),
where I1 and I2 are the MOI of the front and rear wheels, and R is the rolling radius of the wheels.
The big question is whether I1 and I2 are large enough to significantly impact the slope calculation. On one hand, we’d like to make use of as much information as possible to stabilize the slope calculation. On the otherhand, we want to make sure not to encumber the process of applying Virtual Elevation with unnecessary steps. The test protocol is hard enough as it is, after all.
Measuring MOI: The Drop Test
A long time ago, in a galaxy far far away I had a race engineering site for competitive karting called kartuning.com. Don’t bother looking for it — I closed the site quite some time ago due to work pressures. The only remnant of it is this article on the wayback machine that I wrote about the non-insignificant effect of rotational inertia on the dynamics of the kart (WARNING: slow to load).
I introduced the Drop Test as a way to measure the MOI of the rear axle and wheel combination of a kart. The same method could be used on a bicycle wheel. In short, the method goes like this: wrap a string with a weight at one end around the wheel axle. Let it drop and bounce back. Measure the bounce height. The Drop Test code I wrote would compute the MOI based on the original height, the bounce height, and the weight used.
Some Questions For Experienced Virtual Elevation Folk
So, can the rotational inertia terms be used to improve VE calculations? Is there any MOI data for wheelsets out there? Is it worth dusting off the original Drop Test code to compute MOI? Just let me know.
MOI of bicycle wheels can be measured easily when applying the parallel axis theorem (http://en.wikipedia.org/wiki/Parallel_axis_theorem). You use the wheel as pendulum (e.g. pivot point at rim) and measure the periodic time for the pendulum, the mass of the wheel and the distance from the wheel to the pivot point (d). The you ge MOI from
MOI = (T/(2⋅π))**2 ⋅ m_wheel ⋅ g ⋅ d – m_wheel ⋅ d**2
I don’t have data for racing bike wheels but here are at least some values:
MOI = 0.1488 kg·m2 for a touring bike wheel with Schwalbe Marathon 700×35C HS308, Butyl-tube, 36×2mm spokes, Mavic M3CD rim, Deore DX hub (mass of wheel: 1.740 kg)
MOI = 0.0622 kg·m2 for MTB-size wheel with Schwalbe Stelvio 28×559, Butyl-tube, 32×2mm spokes, Mavic X717 rim, Shimano 105 hub, Shimano Ultegra 11-23t-9speed-sprocket
Thanks a lot for that. What a great theorem!
So the effective additional mass added to the acceleration term is:
I/R^2 = ( 0.1488/0.3^2 + 0.0622/0.3^2 ) = 2.34 kg .
Cool! That’s starting to sound like a lot of extra momentum.
The first wheel is a front wheel and the second a rear wheel. In the MOI formula it is devided by 2*Pi (the pi-sympol looks like n). g is the acceleration of gravity.
Three years back, I did the pendulum test on some wheels me and my buddy had. Good wheels have 0.03-0.05 kg m^2.
http://cozybeehive.blogspot.com/2007/11/wheel-rotational-inertia-testing-1.html
http://cozybeehive.blogspot.com/2007/11/wheel-rotational-inertia-testing-part-2.html
Hope that helps.
-Ron
Cozy Beehive bike blog