Win a stage in the Tour de France using Mathematics
Tour de France Mathematics

June 14, 2022

Share this article:

In seventeen days, it is the most wonderful time of the year for cycling fans: the Tour de France will start again. This year, the Tour will cover a total distance of 3328 kilometers through four different countries over a total of 24 days, in which only three rest days will be held. In most of the 21 stages, cyclists will try to break away from the main group of riders, the peloton, create distance between themselves and the peloton and try to win a stage. It turns out that the optimum time to make such a breakaway can be determined using a mathematical model. In this article, the model will be explained. The main result is that the optimum time to make a breakaway is when cycling uphill.

What is a breakaway?

The mathematical model derived by L.H. Gaul, S.J. Thomson, and I.M. Griffiths in the article ‘Optimizing the breakaway position in cycle races using mathematical modelling’ predicts the optimum time for a rider to break away from the peloton. The peloton is the main group of riders in a cycling race. The reason that cyclists ride in a big group is that it reduces the air resistance experienced by the cyclists. As the energy savings from this are quite big, it is favourable for a rider to stay in the peloton. However, to be able to win a stage, a cyclist has to break away from the peloton at some point. As this costs more energy than staying in the peloton, an interesting question is when it is optimal to make this breakaway. This question is answered in the model of Gaul et al.

The model for a single rider

Gaul et al. start by forming a mathematical model for the one-dimensional motion of a single rider. First, some notation is introduced. The motion of a rider is characterized by the distance travelled along the course, \hat{x} at time \hat{t}. Then, the course undulations are characterized by the angle of incline \theta = \theta(\hat{x}) and to describe the motion of the rider in the direction of travel, Newton’s second law is used: \hat{F}_{tot} = \frac{\hat{m}d^2\hat{x}}{d\hat{t}^2}, where \hat{m} is the mass of the rider. Moreover, the force \hat{F}_{tot} is the total force. The forces considered are as follows:

  • The force due to gravity, which is defined as \hat{m}\hat{g}\sin(\theta(\hat{x})). Here, \hat{g} is the acceleration due to gravity.
  • The frictional force, \hat{F}_f.
  • The phenomenological drag force \hat{F}_D = \frac{1}{2}\hat{\rho}c_D\hat{A}(d\hat{x}/d\hat{t})^2. Here, \hat{\rho} is the density of the air, c_D is the drag coefficient, and \hat{A} is the frontal area of the bike and rider. 
  • The pedalling force, \hat{F}_p.

These forces are schematically shown here:

Forces cycling
Source: Gaul, L. H., Thomson, S. J., & Griffiths, I. M. (2018). Optimizing the breakaway position in cycle races using mathematical modelling. Sports Engineering, 21(4), 297-310

Hence, Newton’s Second Law then gives a second-order differential equation for the rider’s displacement:

    \[\hat{m} \frac{d^2\hat{x}}{d\hat{t}^2} = \hat{F}_p - \hat{F}_f - \frac{1}{2}\hat{\rho}c_D\hat{A}\left(\frac{d\hat{x}}{d\hat{t}}\right)^2 -  \hat{m}\hat{g}\sin(\theta(\hat{x})).\]

Solving the model for a single rider

The model is solved by Gaul et al. for \hat{t} \geq 0, where \hat{t}= 0 corresponds to the time one is interested in. Moreover, they impose two initial conditions: \hat{x}(0) = 0 and \frac{d\hat{x}(0)}{d\hat{t}} = \hat{v}_i, where \hat{v}_i is the velocity of the rider. Then, the second-order differential equation is non-dimentionalized in the paper, where parameters are chosen to obtain the dimensionless governing equation:

    \[m \ddot{x} = F_0 - \dot{x}^2  - m \sin(\theta(x)),\]

where a dot means differentiation with respect to time. Moreover, the time taken to travel the total distance of the course X is obtained in the paper by the following equation:

    \[T_{end} = T_0 + \epsilon T_1 + O(\epsilon^2)\]

,

where the details of the functions can be found in the paper.

Modelling the peloton and the breakaway

These equations are now used by Gaul et al. to seek an optimal strategy for a rider to break away by modelling a group of riders and a breakaway. In this model, there is a single breakaway rider. This rider then stops profiting from the drag reduction of the peloton. Here, Gaul et al. consider a course of length X and find the time it takes the peloton to complete the course. They call this T_{end}. The breakaway rider travels with the peloton until some point X_b before leaving the peloton and travelling alone. This also gives the time for the breakaway rider to travel the course: T_{b,end}. The paper seeks for the maximum difference in completion times: \Delta(T_{end}) = T_{end} - T_{b,end} by varying the breakaway position X_b

Gaul et al. use the following equation to capture the motion of the peloton:

    \[\ddot{x} = F_0 - \gamma \dot{x} - \sin(\theta(x)),\]

where the mass factor m = 1, as the mass is scaled by the average mass of the riders in the peloton. Moreover, they use a drag multiplication factor of \gamma = 0.7, as a rider that cycles directly behind another rider experiences a drag reduction of approximately one-third. As a rider escaping from the peloton needs a higher pedal force of F_b, the total pedal force of the rider will be F_p = F_s + F_b. For a single breakaway rider, the following equation holds:

    \[m_b \ddot{x}_b = F_0 + F_b - \dot{x}^2_b - m_b \sing(\theta(x_b - X_b)).\]

Using these expressions, Gaul et al. found the initial extra force a rider will need to provide to accelerate away from the peloton:

    \[F_b(0) > (1 - m_b \gamma) \dot{x}^2 + F_0(m_b-1).\]

Then, the time difference between the breakaway rider and the peloton completing the course is found by   \Delta(T_{end}) = T_{end} - T_{b,end}, where the details of the functions can be found in the paper.

Accounting for fatigue effects

However, these equations assume that a rider does not suffer from fatigue effects. Gaul et al. accounted for the fatigue effect by including fatigue effects in the pedal force. The excess force F_b that a rider must use to break away from the peloton, has to be high enough to also account for the increase in the drag force of a solo rider. Here, Gaul et al. considered the movement of potassium ions out of the muscle cells. They linked the potassium concentration to the force output and included a stamina limitation. This gave them an equation that describes the decay of the pedal force over time after a breakaway. The details of this equation can be found in the paper.

Results of the model

Using these equations, it was found that it is optimal to perform a breakaway when cycling uphill. The reason for this is that the velocities of the riders are lower and the drag reduction of the peloton is thus less significant. It is also best to break away at the beginning of a long hill climb, as the velocity is then expected to be low for a sustained period of time and the rider can thus create distance between themselves and the peloton. However, if an uphill section is followed by a downhill section, or if the uphill section is quite long, it might be optimal to breakaway halfway up the hill. The reason for this is that the extra force can then be sustained up to the top of the hill.

Can you now win a Tour de France stage?

So, using this mathematical model, can you now win a stage of the Tour de France? The short answer is: only if you are a very good cycler. Unfortunately, Mathematics alone is not enough to win a stage in the Tour de France. You still have to train to be a good cyclist, be able to join a cycling team, be chosen to start The Tour, and then still have to have enough power to break away from the peloton. Moreover, you have to be a man (as only men can start The Tour). So, winning a stage of the Tour de France by performing a breakaway is easier said than done. However, you now know the Mathematics behind such breakaways.

The article is based on the following source: Gaul, L. H., Thomson, S. J., & Griffiths, I. M. (2018). Optimizing the breakaway position in cycle races using mathematical modelling. Sports Engineering, 21(4), 297-310.

Read more

Are you living inside a black hole?

Are you living inside a black hole?

In 1915 Einstein wrote the leading description of gravitation in modern physics, also known as the general theory of relativity. This theory describes the universe so well that it never fails, except for the Big Bang and at centres of black holes, both are...

Why Mathematics is a Language

Why Mathematics is a Language

Mathematics is considered the heart of the beta-field, whereas languages are by definition alpha. One could say, therefore, that mathematics and language are direct opposites. Yet, mathematics is sometimes called the language of science and considering the definition...

The rise of pitch invasions

The rise of pitch invasions

With the conclusion of the English football season, emotions always seem to run high. The fear of relegation or the thrill of a promotion or even a championship tend to do that to football fans. After the final whistle it is not rare for fans to spill out of the...