Friday, April 19, 2013

More Continental GP4000S testing...including a 20C


I recently had the opportunity to test additional 23C Continental GP4000S tires, along with a retesting of my original sample after having been ridden ~200 miles as a rear wheel.  I figured this would help give a good indication of both the repeatability of the roller testing and also an idea of the consistency across different tires of the same models.  Here's how it went:

  • 04/05/13 - New tire with ~20 miles of use -    Crr = .00336
  • 04/14/13 - Same tire after ~200 miles of use - Crr = .00343
  • 04/14/13 - Tire used in Flo aero tests -            Crr = .00344
  • 04/17/13 - New tire, fresh out of box -            Crr = .00334

So, across those 4 samples, we get an average of  .00339 (I'd round to .0034, which happens to be the result and number of digits I report in the spreadsheet) and a standard deviation of .00005.

If I'm doing my stats right, then this means there's a 99% confidence range of .0033-.0035.

Granted, this is a fairly small sample set, but it matches pretty well with my "gut feel" that the measurements reported in my Crr spreadsheet should be considered to have a tolerance of around +/- .0001, and that tires listed within .0001 of each other are basically "tied".

I also acquired 20C Continental GP4000 in the black color.  My intention there was to first confirm that the black color GP4000 20C tires have the "Black Chili" tread compound (They do...it says so right on the package), and additionally to see how well it rolls.  The idea was that since it has a similar shape and tread markings as the 23C tire, then it possibly would work as well aerodynamically on narrow rims as the 23C tire appears to do on the wider rims.

The result?  20C Continental GP4000 (Black) - Crr = .0041

That's basically the same as what I found the old 19C Bontrager AeroWing TT tire to exhibit (.0043), in which case, I think I'd still prefer the 20C Continental SuperSonic (Crr = .0034) tire for narrow rims, especially for front wheel uses. As we learned in my last blog post, it would take a LOT of aerodynamic advantage to make up for that much of a Crr difference.

Monday, April 8, 2013

Why Tire Crr matters...

...and why you need to look at BOTH the aero drag of a wheel/tire combination AND the tire rolling resistance to help determine what is "fastest" (i.e. Low Crr can make up for a lot of aero "sins")

(update 04/14/13: Added Michelin Pro4 Service Course to chart after roller testing.  See chart and discussion below)

Five years ago, Damon Rinard (when he was working at Trek) made a post to the Slowtwitch triathlon forum where he described calculating a rough average of aero drag combined with rolling resistance.  The data he used was from some wind tunnel testing he had done with various tires on the same rim (a Bontrager ACC - 50mm deep) and he combined this with the expected "on road" rolling resistance from roller testing.  He did this because he found that different tires made a large difference in the aerodynamic drag of the wheel/tire system, and that simply choosing tires based on low rolling resistance OR low aero drag performance might not be the right approach.  

I thought that was a neat approach and reverse engineered some of his data and expanded the idea to look at the effects of varying wind speed.  In order to do so, I needed to scale the aero drag (taken at 30mph tunnel speed) to different apparent wind speeds using the ratio of V^2/(30mph^2) - since drag force varies with the square of wind velocity.  That value was then summed with the rolling resistance force, which is constant, and then the total combined aero and rolling drag was plotted as a function of expected apparent wind speed. 




Interestingly, when Trek/Bontrager released the R4 Aero tire a few years later (after Damon had left Trek to work for Vroomen/White Design) they included a plot very similar to the one I had created back in 2008:





To use that chart, you need to have an idea of what your general "race speed" is going to be and then figure out what the maximum apparent wind speed (i.e. the vector sum of the ambient wind speed at wheel level and the bike ground speed) you'll encounter.  Then you can see what the average total force you'll be expecting from a particular combination.  Now then, that's talking about the retarding force...if you want to know the power required for the different combinations, or the power differences between combinations, then you need to multiply the drag force value by the ground speed to get the rate of doing work (i.e. power) on that total drag force.  Make sense?

One of the main takeaways from that exercise above I got was that looking at JUST aero drag or JUST Crr wasn't telling the whole story.  It's very easy to "waste" a wheel/tire combination's low aero drag by using a slow rolling tire...and vice versa.  But, the other takeaway I got is that really low Crr can "make up" for a lot of less than ideal aero drag performance

Well, for the plot above, the aero drag component was taken as just a simple average of the 5, 10, and 15 degree yaw angle drags.  After seeing the plot below last fall (taken from the Mavic material given out at the CXR80 tire/wheel system introduction), I realized this type of estimate could be updated using a weighted average instead of a straight average.  The weighting would be from the expected % time spent at each yaw angle.  Obviously, we want to choose a combination based on it's performance under the conditions we expect to mostly see in our races.  If we don't see large yaw angles very often, it might not be worth it to worry about differences in performance at those higher yaw angles between the equipment choices we are contemplating.


According to Mavic, the data taken above is from actual measurements (they built a "wind vane" type rig and attached it to a bike) from a large number of rides under varying conditons and courses.  Although it may not represent the actual yaw angle distribution for any one particular ride (those will likely be skewed one direction or the other, depending on the course and conditions) but it does represent what one would expect over a large number of rides.  As such, it should be a good tool for determining a good "all around" wheel/tire system choice.

Now then, what we need to update things further is some good aero data showing the drags at various yaw angles for different tire/wheel combinations. It would be really helpful if the data happened to be for tires which I've already roller tested and have an idea of the predicted "on road" Crr.

Well, we're in luck.  The guys at Flo wheels went to the A2 tunnel last week and did just that on their new Flo30 wheel.  They tested the Michelin Pro4 Service Course, the Conti GP4000S, the Bontrager R4 Aero, and the Vittoria Open Corsa EVO Tri.  The last 2 were tires that I actually loaned them after Chris Thornham had contacted me asking if I knew of a good place to find the Bontrager tires.  I happened to have a nearly new one handy and also offered to loan one of my Vittoria tires as well.  Here's the blog post describing their tunnel visit and the aero data they took:  Flo30 Aero

As can be seen in their data, the GP4000S was the clear winner aerodynamically, with the R4 Aero close behind it.  The one tire that doesn't look too hot is the Vittoria.  Although it stays fairly close to the other tires up to ~7.5 degrees of yaw angle, after that the drag goes way up.  However, we know from my roller testing (Crr chart) that it's Crr is slightly lower than the other 2 tires, so it might be able to make up for that, especially at lower yaw angles.


 So, let's take a look. What is shown below is the result of taking a weighted average (using the Mavic probabilities for the weighting) of the drag values reported for the 3 tires that I have Crr data on (I'm working on getting a Michelin to roller test as well) combined with the Crr of each tire.  I've used the values of Newtons for the drag force (since it's an actual force unit, as opposed to grams) so that the drag force results can be simply multiplied by the expected ground speed (in meters/second) to quickly calculate the power (W).  (If you want to convert the values to grams, then just divide by 9.81 m/s - gravity - and multiply by 1000)



(Note: the above chart assumes a wheel loading of 38kg and represents a single front wheel)

There's some interesting things going on in that chart.  To understand what's going on there, it's helpful to realize that the "steepness" of each curve is controlled by the aero drag (it varies with the square of the wind velocity), while where the line sits vertically in the chart is controlled by the rolling resistance values (a constant force).

Despite the apparently poor showing of the Vittoria EVO Tri tire aerodynamically, at lower expected apparent wind speeds it actually performs slightly better than the other 2 tires shown; up until ~27 km/hr where it's curve crosses the GP4000S curve.  As compared to the R4 Aero tire, that crossover doesn't happen until expected apparent wind speeds of ~35 km/hr.  At the apparent wind speeds that I expect to see during my TT'ing (~45km/hr) the GP4000S is obviously the leader, with the R4 Aero in second, but the Vittoria still has a predicted average combined drag force that's only .21N higher than the GP4000S.  At the expected ground speed of ~42km/hr (i.e. 11.7 m/s), that results in a total power difference of just 2.5W.  That's really not a very large amount, especially considering how much worse the Vittoria appeared to perform aerodynamically.

Now then, you might be asking "why is that?"  Simply put, a lot of it has to do with the fact that the yaw angles where the largest differences in aero performance occur are also the yaw angles that are weighted less in the aero drag average due to their lower probability of being experienced.  Another interesting thing is that the differences in Crr between the Vittoria tire and the other 2 isn't very large (all within .0004 of each other), but that's enough to overcome a seemingly significantly worse aero performance.

So far we've been talking about this subject in terms of races like TTs and triathlon bike legs where the front wheel is seeing "free air"...but, what about other types of bike racing?  Well, when you're in a group and drafting, the apparent wind speed is going to be lower, along with the fact that the yaw angle distribution will narrow as well...and thus, that makes the Crr component all that more important.  If you ever find yourself in the situation of riding up a false flat in a group while using poor rolling tires, you'll understand what I'm talking about here.  Rolling resistance is actually a higher priority than wheel aerodynamics in road racing, in my humble opinion.

For one last plot, I took the data from the rest of the Flo wheels that were subsequently tested with the GP4000S tire.  That plot is shown below, and as expected at lower expected apparent wind speeds, the wheels are closer together in performance than with higher expected apparent wind speed.



There we are...I hope that helps to understand some of the tradeoffs involved with choosing tires and wheels for cycling. Of course, there are other properties further involved in these sorts of tradeoffs, such as durability and "grip", but those other properties are tough to combine into a chart like the above...and so are left up to the user to weigh separately.

Update 04/14/13:
The guys at Flo were nice enough to send me the Michelin Pro4 Service Course tire used in the aero testing above and so I got a chance to run it on the rollers.  The resultant Crr was .0043.  They also sent the Continental GP4000S used in their testing and it rolled identically to my own version of that tire at .0034.  Using that, I was able to update the Total Drag Force chart to include the Michelin on the Flo 30.  The interesting thing about the curve for the Michelin is that despite it's aerodynamics being fairly close to the GP4000S and the Bontrager R4 Aero (especially with the yaw weighting) on the Total Drag Force chart it NEVER overcomes the "hit" it takes on rolling resistance, even out to expected apparent wind speeds of 60kph.  Once again, the importance of Crr comes to the fore...

 

Friday, April 5, 2013

More Roller Testing Results

Well, I happened to recently acquire a Continental GP4000S tire for a very reasonable price, and in light of the recent Flo aero testing, I decided it would be a good time to throw it on the rollers and see how it does.  I also had a couple of well used "training/all condition" tires a friend loaned me so I tested those as well.  Wow...some of those training tires REALLY sap the power!

Anyway, here's the predicted on-road Crr for the 3 tires I tested today:

Continental GP4000S 23C = .0034
Specialized All Condition Pro II 23C = .0062
Specialized All Condition Armadillo Reflect 25C = .0077

The overall Crr chart now looks like this:


And here is the predicted power for 2 tires at 40kph (85kg total load):


As always, the entire spreadsheet is saved here:  Crr Spreadsheet

Tuesday, February 19, 2013

Tire Crr Testing on Rollers - The Chart..and a "how to"





In my last post I outlined the "math behind the madness" of testing the rolling resistance of bike tires on home rollers. In this one, I'll be showing the results of some of that testing I've personally done over the past year or so. I'll also go through a few tips and tricks I've learned in doing this sort of testing...just in case anyone else is crazy enough to try some of this potentially mind numbing testing.

OK, I know a few of you are out there are "champing at the bit" to see the results, so without any further ado, here's a chart showing my estimates for the power to roll a pair of various tires I've tested (on a "real road" and for an 85kg bike plus rider mass):




Here's the same chart, but showing the estimated "on road" Crr values:

As can be seen, that's a fairly wide range of power requirements.  The wrong tire choice can easily "cost" a rider 10-15W of power to go a given speed.  When choosing tires, I commonly think of a scene from "Indiana Jones and the Last Crusade" where the Templar Knight guarding the Holy Grail says "...you must choose, but choose wisely...".  After all, when you're done with a race and you lose by seconds, or inches...you would hate to have the following be said about your tire choice:



The Setup:

OK, now that we've got that all out of the way, I thought I'd describe a bit about the setup I use for doing this sort of testing.  As you can see in the pic at the top of this blog post, it's a fairly simple affair consisting of a set of 4.5" diameter Kreitler rollers, a front fork stand, and a bike equipped with a power meter (and a power meter head unit).  That's really basically it.  A couple of other pieces of equipment that are crucial for getting consistent results, in my experience, are:
  • A means to measure ambient temperature
  • A means to measure rear wheel load
  • A separate speed sensor and magnet (NOT on the wheel, but on the roller - see below)
  • A notebook and pen
For measuring the ambient air temps during the test, I use my trusty Brunton ADC Summit, which I place at about axle level somewhere near the side of the rear wheel of the bike during the testing.


To measure the rear wheel load of the setup, I actually just use a digital bathroom scale which I've checked against known weights and typically is within 0.5 lbs of the actual weight.  In order to make that rear wheel load measurement, I mount the fork in the fork stand, and in stead of placing the rear wheel on the rollers, I stack the scale on top of some wood scraps and place the rear tire on the scale.  I stack it so that the rear axle is the same height off the ground as it would be on the rollers.  As it turns out, a couple of scrap pieces of 3" square wooden post and the thickness of my scale put that measurement to within 1/4" of what it is on the rollers.  Perfect.

Now then, let's talk about the speed sensor I mentioned above.  One of the most important things to get an accurate measure of in this testing is the actual "ground speed" during the test.  This can be done with a wheel mounted magnet and speed sensor, BUT that requires determining and changing the wheel rollout number for EACH tire tested.  In my experience, that can be a bit problematic...especially due to the curved contact patch that is present on the rollers.  It's very hard to get an accurate and consistent measure of wheel rollout that way.  To solve that problem, I realized that instead of triggering a speed sensor on the wheel, I could instead attach a magnet to the end of one of the rear drums and then use a speed sensor triggering off of the drum.  All I needed to do was to carefully measure the diameter of the metal drum (which will NOT be changing from test to test) and use THAT as the ground speed measurement for the testing.  Here's what that looks like:


That's just a small rare-earth magnet attached to the end cap of the roller with double-sided tape, and a Garmin ANT+ speed/cadence sensor taped to the roller frame.

Lastly, the notebook and pen are where I write down the date, what tire I'm testing, the size, the ambient temps during the testing, the power meter zero offset numbers, and the actual measured tire width as mounted.


The Protocol:

Alright, so everything is set up and we've gathered all the stuff needed.  What's next?  How do I do this? Well, here's a quick rundown of how I go about doing a tire Crr test (in "10 easy steps"!) This isn't the only way to do it, but it's how I've settled on things after doing this testing for a while:
  1. Mount the tire on the test wheel - Most of my testing is done with my old yellow-cap PT wheel with a Mavic Open Pro rim.  I started out testing with this wheel in order to get both a hub and crank power to determine the level of typical drivetrain losses in the setup.  I wanted to know that for the occasions when I would be testing tires (such as tubulars) which I couldn't mount to the clincher PT wheel. Since I'm mainly interested in tires for time trials and road racing, I'll test them using a latex inner tube.  Testing by others has shown that using a butyl tube can cause 10-15% higher Crr than with latex.
  2. Pump the tire to the test pressure - What pressure to use is really up to you.  I chose to do all of my testing with 120psi.  The reason I chose that value was mainly so that I could compare my results more easily to the results of others, most notably the testing done by Al Morrison.  Understand that on a perfectly smooth surface, the higher the pressure you pump tire up to, the lower the measured power requirements will be...however, that will only be true on the rollers, or on flat surfaces that are just as smooth.  On "real roads", i.e. roads with typical roughness, that isn't necessarily the case and there will tend to be a pressure above which higher pressures actually will make you slower overall.  Anyway, the key here is to pick a pressure and stick with it through your testing so that you are comparing tires on an equivalent basis.  
  3. Place the wheel in the test bike - Install the rear wheel in the test bike and place the chain in the chosen gear for the testing.  I do my testing in a 53x13 gear for consistency. If it's not in there already, install the fork into the fork mount.
  4. Measure the rear wheel load - This is done how I described above.  I'll usually only do this once during a session, and for me I've found it's typically within a pound or two each time (my body weight tends to be fairly stable).  This doesn't seem to be a super-critical measurement either, since a 1 or 2 lb. difference will only result in ~1-2% error in the final calculation.
  5. Place rollers under rear wheel - At this point, move the scale out from under the rear wheel and slide in the rollers. To get a consistent placement of the rear wheel on the rollers, I'll lift the fork mount slightly off the ground while I allow the rear wheel to spin as it touches the 2 rollers and then I carefully place the fork mount on the ground.
  6. Climb on board - It's now time to saddle up.  I usually approach the bike from the non-drive side and put my left foot on the pedal and then carefully swing my right leg over taking care not to disturb anything in the setup.  I'll then spin the cranks to make sure the PM is awake and the speed reading is working, at which point I'll clip out and zero the PM through the head unit.  I'll note the offset value (from my Quarq) in the notebook along with the ambient temperature.
  7. Tire warmup - Now it's time to warm up the tire to working temperature.  Since my tests are done at 90 rpm (I find it's easier to hold a consistent rpm rather than focusing on the wheel speed) I'll warm up the tire at 95 rpm for 5 minutes.  At the end of the 5 minutes, I'll stop and quickly check the PM zero offset and write the value down in my notebook along with the ambient temperature reading.
  8. The Test - Now it's time for the test.  I'll bring the cadence up to 90 rpms and once that is steady, I'll start a 4 minute interval in the PM head unit.  I'll concentrate on keeping a steady cadence through the whole interval, trying to be especially steady through the final 2 minutes since that is the section of the data I take the average power and ground speed from.  At the end of the 4 minute test interval, I'll again note the PM offset (to make sure it hasn't moved appreciably during the test for some reason) and write down the ambient test.  That's it.  Test over...either I stop there, or if I have more to tires to test, I'll start back at step one (skipping the load measurement for repeat tests) and on through the remaining steps.
  9. Download Data - Now it's time to get the average power and ground speed values from the head unit.  I'll typically load the file into Golden Cheetah and then highlight the final 2 minutes of each test session and read off the averages as calculated.
  10. Calculate the Crr - The final step is to take the average power and speed values, along with the wheel load and ambient temperature taken at the end of the test interval into a spreadsheet I've written to quickly do the calculations.  I've loaded the spreadsheet onto Google Drive and it can be accessed here: Crr Spreadsheet
Other Notes:
  • After doing a number of runs using both the crank-based power meter and the PT wheel, I was consistently finding that for the gearing chosen and the lower power levels (typically 50-100W) seen on the 4.5" rollers, the drivetrain losses were on the order of 5%.  That's the value I enter in the spreadsheet.
  • I did a fair number of runs with the exact same tire and at different ambient temperatures to determine what I should use as the temperature compensation value.  In my case, I found it to be ~1.36% change per deg C (lower Crr with higher temperatures).  Here's the plot of those tests:

  • I normalize the Crr values to 20C.  If you want to know what the Crr would be at various temperatures, you can just enter those temperatures into the appropriate cell on the first sheet of the spreadsheet.
  • Added 2/19/12 - I realized I forgot to point out that I use a "smooth to real road" factor of 1.5X to account for the higher energy dissipation requirements of typical road roughness.  This value is based on comparisons of roller based Crr measurements ("translated" to flat surface) and actual "on road" Crr derived from field tests and other means (i.e. iAero coast down values) for the same tires.

Well...that's about all I can think of for now.  Hopefully that will help encourage others to give this a try.  It's really not that difficult to do and is a good way to help you to "choose wisely" when it comes to tires for your "go fast" bike setup.

Saturday, February 9, 2013

Tire Crr testing on Rollers - The Math


I've been doing a little bit of tire testing lately.  But, before I reveal any results, I thought it would be good to go over some math.  I know, I know...(I can hear the groans already), but I think it's important to review so people understand why it's reasonable to equate power to move a tire on a roller to power on flat ground.  

It's long been known that bicycle rollers act as a sort of rolling resistance "amplifier".  In other words, the differences in the rolling resistance between tires is magnified when riding on the rollers.  It's usually a fairly subtle thing to try to "feel" the difference in rolling resistance in tires when riding outside, but it's pretty easy to tell the fast tires from the slow tires on rollers just by the exaggerated effort it takes.  But, the question has long been "how much" of an amplifier are they? Well, back in 2006 I was discussing this with a few folks and realized that the equations to make that comparison between rollers and a flat surface Crr (Coefficient of rolling resistance) were already available...they just needed to be combined.  Then, it was pointed out to me that the particular geometry of a typical roller setup needed to be accounted for as well.  The normal "dual roller" setup on the rear of a roller set results in a geometric effect that actually increases the normal force on each roller.  In other words, you can't just take the rear wheel load as if it was a single roller. So, I added that to the equations as well.

Anyway, what you see below is the short "paper" I sketched up back then on the subject:

Flat Surface RR from Roller Testing – Tom Anhalt – 5/2/06



The power required to turn a wheel on a drum at a specific speed is governed by the equation:



PDrum = CrrDrum x VDrum x M x g (a)



Where,

PDrum = Power required to turn drum (Watts)

CrrDrum = Coefficient of Rolling Resistance of the tire on the drum (unitless)

VDrum = The tangential velocity of the drum (m/s)

M = The mass load of the wheel on the drum (kg)

g = gravitational constant = 9.81 m/s2





Rearranging equation (a) to solve for the Crr of the tire on the drum results in:



CrrDrum = PDrum / (VDrum x M x g) (b)





Then the contact patch deformation of a tire of a specific diameter and a roller of a specific diameter can be equated to the deformation of an equivalent diameter tire on a flat surface using the following equation [Bicycling Science, 3rd edition, pg 211]:



1/req = 1/r1 + 1/r2 (c)



Where,

req = equivalent wheel radius

r1 = tested wheel radius

r2 = tested drum radius



For convenience purposes, this equation can be rewritten using the appropriate diameters (r x 2) and is then:



1/Deq = 1/Dwheel + 1/DDrum (d)



For a tire of a given construction, it has been shown that the Crr varies inversely proportionally to the wheel radius, and thus the wheel diameter, in the range of Dwheel0.66 to Dwheel0.75 [Bicycling Science, 3rd edition, pg. 226]. To simplify for this purpose, the assumption is made that the Crr varies inversely proportionally to Dwheel0.7



From this, it can be then written that:



Crrflat / CrrDrum = Deq0.7 / Dwheel0.7 (e)





Equation (e) can be combined with (d) and rearranged to give:





Crrflat = CrrDrum x [ 1 / (1 + Dwheel/DDrum)]0.7 (f)







Substituting equation (b) for CrrDrum in equation (f) results in:





Crrflat = [PDrum / (VDrum x M x g)] x [ 1 / (1 + Dwheel/DDrum)]0.7 (g)







Mass Correction Factor:



When doing Crr testing on rollers, the mass loading of the wheel or wheels will need to be corrected due to front-rear loading ratio and the fact that 2 offset rollers contact the rear wheel, thereby increasing the normal force on the rollers due to geometry effects.



Rear Wheel Only Case - When the test is done using a front fork mount and only the rear wheel contacting the rear rollers of the test setup, the following “effective mass” (Meff) needs to be calculated and substituted for M in equation (g) :



Meff = Mrear / cos [arcsin (X/(Dwheel + DDrum))] (h)



Where:

X = separation distance of rear roller axles (consistent units with Dwheel and DDrum)

Mrear = vertical mass load on rear wheel (kg)





Front and Rear Rollers - When the test is performed using both the front and rear rollers, the following Meff needs to be calculated and substituted for M in equation (g) :



Meff = Mfront + Mrear / cos [arcsin (X/(Dwheel + DDrum))] (i)



Where:

Mfront = vertical mass load on the front wheel (kg)





Power Correction:



Depending on the method of power measurement, the following offsets can be used to account for drivetrain and drum rotation losses in the calculation of PDrum for use in equation (g):



For Powertap - PDrum = PPowertap – 5W (accounts for drum bearing losses) (j)



For SRM - PDrum = PSRM – 15W (accounts for drum bearings and driveline losses) (k)



Where:

PPowertap and PSRM are the power readouts (W) from the appropriate power meters.



These power offsets are somewhat arbitrary and should be modified if better data is known about the particular test setup.




That's basically it.  I'd like to note a couple things about the last section on "Power Correction".  First, after doing a bunch of testing since then, I don't bother accounting for the drum bearing losses.  Also, when using a crank-based power meter, like an SRM or a Quarq Cinqo, I've found (after using a PT in conjunction for quite a few tests) that it makes more sense to account for the drivetrain friction with just a straight percentage.  A typical value taken for drivetrain losses on bicycles is ~2.5%, but that is typical at higher power levels (i.e. higher chain tensions).  For the lower power levels I usually see in tire testing with a rear only roller setup (usually ~100W or less, with the better tires closer to ~50W) I typically see ~5% drivetrain losses, so that's the figure I use.

It's important to remember that the point of this is to get a "ballpark" feel for the difference between tires, not necessarily an absolutely accurate value.  It's been shown that percent difference in power requirements on the rollers equate very well to percent differences on the road.  What we're really looking for is a sort of "scaling factor" to put the differences seen on the rollers in perspective as to what to expect on the road.

There you go...equation (g) is easily written into a spreadsheet.  After that, it just takes a few measurements of the roller setup, weighing the rear wheel load, and some time on the rollers with a power meter equipped bike and nearly anyone can "test tires".










Friday, February 1, 2013

LeMond Power Pilot - Does it give good numbers?




Considering that it was announced today that Greg LeMond had formed a new venture to sell the Revolution trainer, and that my last post had dealt with the Revolution, I thought it would be a good time to also talk about the Power Pilot device that LeMond sells for use the trainerIf someone already has a non-wheel based power meter on their bike, then the Power Pilot would be a bit redundant.  But, for those who use a PT wheel primarily, or don't have another form of power measurement, the Power Pilot could be a good alternative for determining the "load" during a trainer workout.  The following is a brief look at the power reporting and recording performance of the LeMond Power Pilot. In particular, it is compared to the output of a “known good” Quarq CinQo crank-based power meter.  "Known good" in this sense is a power meter which has had it's torque slope checked and adjusted and has a zero offset that is stable.



BACKGROUND

The LeMond Power Pilot is a device designed to be used in conjunction with the LeMond Revolution trainer to primarily monitor and record the training efforts of the rider. The Revolution trainer is a high inertia wind trainer with the somewhat unique configuration whereby the rear wheel of the attached bike is removed and is not a part of the driven assembly. The chain of the bike drives a rear cassette which is attached to a relatively high mass flywheel through a belt drive gear reduction system. In my last blog post ("What's the Virtual CdA and Crr of the LeMond Revolution Trainer"), it was shown that this system “mimics” the aero drag of a typical sized rider on a road bike (CdA = ~0.35 m^2) and the rolling resistance of average tires (Crr = ~ .005). The inertial mass was also found to be equivalent to a rider mass of ~45kg, which although it is less than the mass of a typical rider, it is far higher than the much lower inertial masses of most indoor trainers on the market today. This accounts for the often reported excellent “road feel” of the Revolution trainer.


The Power Pilot uses an ANT+ sensor to read the drive pulley rotational velocity, and by extension, the flywheel speed during operation. Along with the known aero properties of the flywheel fan, the aero drag power is calculated using that speed measurement and an estimate of the air density (based upon the user entered altitude, an internal temperature sensor, and an internal humidity sensor.) The speed sensor is also used to determine the acceleration/deceleration of the flywheel mass to account for that in the power reporting. Additionally, the Power Pilot firmware allows for a coastdown calibration to be performed which then accounts for any unit to unit variation in the “fixed” losses of the mechanism.



THE TEST AND RESULTS

In order to determine how well the Power Pilot performs these calculations, it was decided to compare the Power Pilot output to the power values measured using a Quarq CinQo crank-based power meter. The particular CinQo used in this testing is a “known good” unit which has been recently calibrated for torque slope and demonstrates a very stable zero offset. To make the comparison, a ride was undertaken whereby a “cassette sweep” was performed. The ride started with the bike in a gear selection of 53/25 and then progressed down the cassette every ~1.5 minutes until a gear selection of 53/12 was used, all the while keeping a constant 60 rpm cadence. Then, the chainring was shifted to a 39T and the progression was repeated partially up the cassette. Finally, two additional runs were taken at a higher cadence (and thus power). A plot of both power traces vs. time is shown below:





Looking closely, one can see that aside from a slight offset, the values reported by each power meter “track” very closely to each other. This can further be seen if the power values relative to each other are plotted on a point-by-point basis. The plot below shows the power reported at each point in time with the CinQo power values on the x-axis and the Power Pilot values on the y-axis.






Typically, due to factors such as variations in recording rates, calculation algorithms, etc. the sort of plot above doesn't turn out very well as a comparison tool for power meters. However, in this case, the point-by-point power reporting is fairly good and it appears that the Power Pilot reports ~94% of the power reported by the CinQo. Considering that the Power Pilot estimate is taking place “downstream” of the bicycle drivetrain losses (much like a PowerTap wheel) in comparison to the crank-based location of the CinQo, a ~6% drivetrain power loss is reasonable, and typical of an average drivetrain.


Rather than looking at the point-by-point plot and it's curve fit, often it's more useful to look at a plot of the averages reported by the power meters over constant power sections. The plot below is the same as the previous plot, but the averages over each of the “intervals” of the cassette sweep are plotted and a curve is fit to them. As can be seen the fit is fairly “tight” to the data and the slope of the fit (i.e. the drivetrain loss) is similar to that reported above, albeit slightly lower. This appears to point out that the correlation between the 2 devices may be slightly better during constant power efforts. Alternatively, it could reveal that the time correlation of the point-by-point data looked at above may be slightly “off”.







However, understanding where the power meters don't agree can sometimes be much more enlightening in that it can reveal systemic difference between the devices. One tool for doing this is something called a “Mean – Difference” Plot. This plot shows the difference between the power values reported by the 2 devices plotted against their average. The Mean-Difference plot of this ride is shown below:







The things to look for in this type of plot that reveal systemic “issues” are characteristics like the spread of the points getting larger or smaller with the mean power values or the difference not trending in a monotonic fashion. Neither of these types of issues appear to be present above, implying that the Power Pilot “agrees” with the CinQo output in a consistent and predictable manner.


Another area of interest for comparisons of this type is how the 2 devices respond to sprint type efforts. An additional ride was undertaken where a pair of short sprints were performed in order to see if there were any anomalies between how each device reported these efforts. As can be seen below, the shapes of the curves are very similar with the peak values also being very close to each other.







CONCLUSIONS

This brief examination of the power output of a LeMond Power Pilot on a Revolution trainer shows that it does an acceptable job at determining and recording the rider's power output when compared to a “known good” crank-based power meter. The power output is similar to what one would expect to see from a hub-based power meter, which is favorable for the likely customer base for this product, namely users with no other power meter or Power Tap owners who want to train indoors with and/or by power.

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.