Intro to Power Estimate Modelling

It’s just physics, but it’s so much more.

The original description of a complete mathematical model for estimating cycling power output has to go to Martin. But it was Antoine Vayer’s snarling bulldog sound bites and his collaboration with Fred Portoleau on Not Normal Doping that brought the power estimates and the debate over what the numbers meant into the mainstream’s view.

Although the debate over the interpretation of the numbers is truly just getting started, the method itself is not quite the pseudo-science that some would prefer you believe. Mainly the method is based on relatively straight forward physics; the total power is the sum of the power to overcome gravity, wind resistance, rolling resistance, and drive train drag. Yes that’s basically it.

The details, of slightly simplified version of the equation, look like this:

estimate w/kg equation eW/kg = wind w/kg + gravity w/kg + rolling resistance w/kg
wind w/kg wW/kg =((0.5*CdA*air density)*(velocity^3))/(wieght*0.976)
drive train efficiency 0.976
rider weight kg 67
bike weight kg 8
1/2 draft CdA CdA = full draft CdA + no draft CdA / 2
full draft CdA full draft CdA = 0.75 x no draft CdA
no draft CdA no draft CdA = Cd x area
Cd Cd = 0.84
area area =((0.0276*( heightt (meters) ^0.725))*(weight (meters) ^0.425))+0.1647
air density air density = 1.225 x rho1/rho2
rho1/rho2 rho1/rho2 = (( 1-(((gravity*(n-1))*(mid altitude -1) )/((n*R)*Temp1) ) )^(1/(n-1)))*(Temp1/Temp 2)
n polytropic gass constant 1.235
R Nm(kg*K)^-1 287.1
gravity 9.81
Temp 1K (15 deg C) 288.15
Temp 2 K (20 deg C) 293.15
gravity w/kg gW/kg = ((gravity*(rider wieght+bike weight))*vertical velocity)/(rider weight *0.976)
vertical velocity meters / sec vertical velocity = vertical meters / time
rolling resistance w/kg rrW/kg =(((((velocity * CRR)* gravity)*(rider weight + bike weigth)))+((velocity*(91+(8.7*velocity)))/1000))/(0.976 * rider weight)
Crr 0.004
velocity m/s velocity = distance / time

What is most important though is how the math holds up in the real world. To find out Vetoo collected 250 power meter files from 36 Pro-tour riders over 250 climbs with adequate video footage to mathematically estimate power-outputs. I then took this data set and ran some basic statistics in R and the results were quite good.

Martin 0.35

As you can see in the figure above, the mathematically estimated power correlated well with the directly measured power meter data.

Statistically the estimate are “likely” to be within +/- 1.6% (75% confidence interval) of power meter measurements and “very likely” to be within +/- 2.7% (95% confidence interval). For perspective, the “gold standard” SRM power-meters are only reported to be accurate to within +/- 1-2% themselves.

kernel QQ

A quick quality check, shows that the residuals are have a reasonably normal distribution without much weirdness going on.

So unsurprisingly to me anyways, it turns out that there is nothing magical about the physics of pro-cyclists. They too follow at least some of the fundamental laws of nature.

Now it is true strong consistent head or tail wind can throw an estimate off by 5-10% or more. But as we can see from the 250 climbs above, this scenario looks like a rare event in Pro-tour races.

The last point to keep in mind, is that most Grand Tours are going to have at least 3-4 critical climbs that decide much of the GC standings. Since the error has a fairly normal distribution, it gets increasingly cancelled as the number of climbs increases.

 

 

About