Saturday, January 26, 2013

What's the Virtual CdA and Crr of the LeMond Revolution Trainer?



...and, more importantly, what's the "equivalent mass"?

Back in 2011, I was given access to a LeMond Revolution trainer.  One of the selling points of this trainer is it's "road feel", and knowing that the load was produced by air drag, I decided to see if I could figure out what the "virtual" CdA (drag coefficient) and Crr (rolling resistance coefficient) was for the trainer, just to see if it reasonably represented a typical rider and equipment.  Another part of the "road feel" is the inertial mass of the trainer's flywheel, and I was hoping I could figure out the "equivalent mass" that the flywheel represented.  In other words, what mass of a rider (plus bike) traveling down the road does the spinning flywheel mass represent?  Some of you reading this may recognize the figures from a thread I started back then on the Slowtwitch.com forum, and also one on the wattagetraining.com forum.

The majority of the load experienced by a cyclist when riding outside on level ground is produced by aerodynamic drag and the rolling resistance of the bicycle equipment (tires mostly, but bearings as well...) and is very simply represented by the following equation (assuming still air as well).

Power = Aero drag power + rolling resistance power, or

Power = (1/2 x air density x CdA x V^3) + (Mass * g * Crr * V),
 where V is the velocity of the bike and g is the gravitational constant.

This equation is handy, because if you divide both sides by V, you end up with a linear equation with respects to the variable V^2.  In other words, the equation of a line is:

Y = mx + b, or dependent variable = (slope * independent variable) + (y-intercept)

What that means is: if you are able to measure the power it takes to travel at various bike velocities, then if you plot the Power/Velocity on the y-axis and the Velocity^2 on the x-axis of a chart, you should see a line where the slope is equal to the CdA and the y-intercept is the Crr.  In fact, this technique is a common one used in field testing and is sometimes known as the "regression method", since it allows one to easily do a linear regression fit to the data.  I just thought I'd apply it to the wind trainer as well...

Here's the result of a constant cadence cassette sweep (only 9 of the 10 cogs of a 12-25 cassette, 60 rpms, with a repeat of the 53-14 at 75 rpm for some slightly higher power.) After the sweep, I did some accelerations/decelerations for the purposes of doing some inertial mass estimate.


Taking the average of the power and "virtual speed" (i.e. the result of the gearing, cadence, and assumed wheel rollout) over the last 2 minutes of each step (and the last 1 minute of the 75 rpm step), I then plotted P/V vs. V^2 for the following:



That looks pretty good!  Nice and linear, which means it basically "behaves" like a rider on the roadAssuming an "all-up" mass of 85kg, and a rho (air density) of 1.2 kg/m^3, that y-intercept works out to represent a Crr = .0051 and the slope of the line works out to represent a CdA = .350 m^2. Sounds like a fairly "normal" road bike position (on the hoods) and Crr.  Nice.

Now...about that estimate of inertial mass. My intent was to plug the file and the calculated Crr and CdA into my VE ("virtual elevation") spreadsheet and then modify the mass entry until the small "hills" formed by the accelerations/decelerations "flatten out"...I tried that, but I realized that I needed to reconsider how varying the mass entry affects the calculated rolling resistance force.  See below.  In a perfect representation, that trace shown below would be completely level.  The fact that it's not means that something is amiss with the modeling, probably the acceleration terms.

BTW, "virtual elevation" is a technique of analyzing power files created by Robert Chung which allows one to estimate drag coefficients given a certain Crr, or vice-versa. See a presentation on the subject by Robert Chung here.


At about this time in the Slowtwitch thread, Robert asked to take a closer look at the acceleration sections near the end.  Here's what that looked like:

  

He then asked to see the VE trace in an overlay on the speed trace and so this is what I showed him:


At this point he said: "I presume you were using an all-up mass of around 80kg? So cut that in half to 40kg and double the Crr to .01. That'll keep the steady state power the same but should improve the modeling for the KE component."  

Here's what it looked like with the mass decreased by half to 42.5 kg and the Crr increased to .0105 (I played with the Crr to get the best "fit"):   

The VE changes between the gear changes at least seemed to have been "smoothed" over...but the accels/decels at the end were still a bit "spikey".
 
To that, his response was: "Cool. In the original, the VE was "countercyclical" to the speed. Now it's cyclical. I think that means halving the mass was an overcorrection. So the inertial mass is somewhere between 42.5 and 85."

I then played around with the mass and Crr a bit and here's the "best fit" I could find...and I mostly judged it off of the final "tail" which was a coast down to zero rpm. I adjusted the mass (and Crr accordingly) until the final steady-state leg and the coast down to zero had a minimum of "disjoint". This case was actually 100 lbs (45.5 kg) and a Crr of .0094 


Now...we'll see how close that is to a calculated mass based on the flywheel geometry! 

I then pulled the flywheel cover off and measured up the flywheel on the trainer. To calculate the mass moment of inertia, I decided to just throw the dimensions into Pro/E and let it do the calculations:  


Here's the mass properties output:

VOLUME = 8.0428750e-04 M^3
SURFACE AREA = 2.0893183e-01 M^2
DENSITY = 7.8887728e+03 KILOGRAM / M^3
MASS = 6.3448414e+00 KILOGRAM

CENTER OF GRAVITY with respect to _GL_FLYWHEEL coordinate frame:
X Y Z 0.0000000e+00 0.0000000e+00 2.4937030e-02 M

INERTIA with respect to _GL_FLYWHEEL coordinate frame: (KILOGRAM * M^2)

INERTIA TENSOR:
Ixx Ixy Ixz 4.2069766e-02 0.0000000e+00 0.0000000e+00
Iyx Iyy Iyz 0.0000000e+00 4.2069765e-02 0.0000000e+00
Izx Izy Izz 0.0000000e+00 0.0000000e+00 7.3130437e-02

INERTIA at CENTER OF GRAVITY with respect to _GL_FLYWHEEL coordinate frame: (KILOGRAM * M^2)

INERTIA TENSOR:
Ixx Ixy Ixz 3.8124192e-02 0.0000000e+00 0.0000000e+00
Iyx Iyy Iyz 0.0000000e+00 3.8124191e-02 0.0000000e+00
Izx Izy Izz 0.0000000e+00 0.0000000e+00 7.3130437e-02

PRINCIPAL MOMENTS OF INERTIA: (KILOGRAM * M^2)
I1 I2 I3 3.8124191e-02 3.8124192e-02 7.3130437e-02

ROTATION MATRIX from _GL_FLYWHEEL orientation to PRINCIPAL AXES:
1.00000 0.00000 0.00000
0.00000 1.00000 0.00000
0.00000 0.00000 1.00000

ROTATION ANGLES from _GL_FLYWHEEL orientation to PRINCIPAL AXES (degrees):
angles about x y z 0.000 0.000 0.000

RADII OF GYRATION with respect to PRINCIPAL AXES:
R1 R2 R3 7.7515746e-02 7.7515748e-02 1.0735906e-01 M

The value we're interested in is the the Izz value of .073 kg*m^2.

So...how do we equate this to an "equivalent mass" translating down the road? I like to do it by equating the kinetic energy of the flywheel to the kinetic energy of a bike+rider moving down the road.

Kinetic Energy of rider = Kinetic Energy of the Flywheel
KErider = KEf
1/2 x mass of rider x (velocity of bike)^2 = 1/2 x Izz x (flywheel rotational speed)^2
Mr x Vb^2 = Izz x (Wf)^2

OK...to solve this, I need to equate the flywheel rotational speed (in radians per second) to the equivalent bike velocity (in meters/second). Well, the assumption above was that the wheel rollout was 2080mm, or 1 revolution = 2*Pi radians is equivalent to 2080mm of rollout.

Wheel rotation rate = Ww = Vb x (2*Pi radians/2.080m), where Vb is in m/s, so Ww = 3.0208 x Vb.

Now, to figure the flywheel rotational rate, we need to know the pulley ratio of the drive pulley and the flywheel pulley. By my measuring, this ratio is 8:1. So, the flywheel rotational rate, Wf = 8 x Ww = 8 x 3.0208 x Vb = 24.166 x Vb.

Lastly, I'll plug this relation into the simplified KE equation above along with the calculated Izz from the solid model.

Mr x Vb^2 = Izz x (24.166 x Vb)^2
Mr = Izz x 584 = .073 x 584 = 42.6 kg

42.6 kg is equivalent to a rider weight of ~94 lbs...now, that's also not including the other rotating bits (like the pulleys, cassette, etc.)...but that's pretty darned close to the "equivalent mass" determined using the coastdown and VE above :-D  


So...in the end, what does that all mean?  Well, I'll end with one more quote by Mr. Chung: 

"Aha. Excellent. It may "coast" better than regular trainers -- but not as well as a bike on a flat road. Likewise, it appears to accelerate faster than a bike on a real road. Steady state load is reasonable, though."

Sounds good to me.

No comments:

Post a Comment