Brian Beckman, PhD
©Copyright February 2001
In this instalment of the Physics of Racing, we complete the program begun last time to combine the magic formulae of parts 21 and 22, so that we have a model of tire forces when turning and braking or turning and accelerating at the same time. Parts 21 and 22 introduced the magic formulae. The first one takes longitudinal slip as input and produces longitudinal grip as output. The other one takes lateral slip as input and produces lateral grip. Slip depends primarily on driver inputs, grip is force generated at the ground. Longitudinal means in the straight-ahead direction. Lateral means sideways, as in the forces for turning. Since the magic formulae work only in isolation, we have work to do to model turning and braking at the same time and turning and accelerating at the same time.
Last time, we vectorized slip - the input - to come up with combination slip, captured in the vector slip velocity. That vector measures the velocity of the contact patch with respect to (w.r.t.) the ground in one, handy definition. This time, we first boil down combination slip to new inputs for the old magic formulae. In the old magic formulae, we measure longitudinal slip as a percentage of unity, that is, as a percentage of breakaway sliding; and we measure lateral slip as an angle in degrees. These are not commensurable, meaning that we do not use the same units of measurement for both kinds of slip. That's why there was a big, fat question mark in the vector slot for combination slip in one of the tables in part 24. Once we make them commensurable, then we stitch the magic formulae together to get one vector gripping force as a function of one vector slip. This finally allows us to compute the forces delivered by a tire under combination control inputs.
Once again, we are in uncharted territory, so take it all in the for-fun spirit of this whole series of articles. I don't represent anything I do here as authoritative racing practice. I only claim to be bringing the fresh perspective of a stubbornly naïve physicist to the problems of racing cars as an amateur. The standard practice of the professional racing engineering community may be completely different. This is the Physics of Racing, not the Engineering of Racing. I'm after the fundamental principles behind the game. I use techniques that may be foreign to the engineers that build and race cars professionally. My results may not be precise enough for final application. I may take approximations that simplify away things that are actually critically important. On purpose, I'm figuring things out on my own. Often, this helps me understand published engineering information better. Just as often, it helps me debunk and debug the conventional wisdom. If you find mistakes, gaffs, or laughable dumb stuff, or if you know better ways to do things, I encourage you to fire up debate, publish rebuttals, or write to me directly. I've done my best to track down the latest and greatest information, but I've found lots of errors, ambiguities, and inexplicabilities in the open literature. I also suspect a conspiracy, meaning that I'd bet that the tyre manufacturers and pro racing teams don't publish their best information-I certainly wouldn't if I were they.
Disclaimers out of the way, we now have enough tools on the table to combine
the two magic formulae. Recall the formulae from parts 21 and 22: and
for the longitudinal and lateral forces. Here
they are, in isolation:
There are a lot of ways we could stitch them together. This is not the kind
of situation where there is one right answer. Instead, in the absence of hard
theory or experimental data, we have the freedom to be creative, with the
inevitable risk of being wrong. We pick a method that satisfies some simple,
intuitive, physical requirements. First, we must put the inputs on the same
footing. Ask "what is the value of
for which
has its maximum, and what is the value of
for which
has its maximum?" Call these two values
and
. They are constants for given Fx and
: characteristics of a particular tyre and car and surface. So, we can
finesse the notation and just write
and
. The maxima identify points on the rim or edge
of the 'traction circle'. The grip decreases when
exceeds
and when
exceeds
. Let's illustrate with
,
= 0, and the constants from Genta's
alleged Ferrari. Once we substitute all that in (and we'll let you check our
arithmetic from the data in prior articles), we get
We evaluate these equations for
= 0,
= 0, getting
,
, and showing a small lateral force (about 16
lbs) due to conicity and ply steer. The source of that problem is the constant
offset in S, which results from a9 and
a10's being non-zero. We just set them to zero for now. Let's
plot
, slip on the horizontal axis and grip on the
vertical:
The maximum positive grip occurs, just by eyeball, around
= 0.08. To the left of the maximum, adding more slip - more throttle -
generates more grip. To the right of the maximum, adding more slip generates
less grip. That's where we've lost traction. We can find the maximum
precisely by plotting the slope of this curve, since the slope is zero
right at the maximum:
Using secret physicist methods, I've found that this curve crosses the
horizontal axis - that is, goes to zero - at precisely
= 0.0796. This was so much fun that we'll just do it again for
. First, the curve proper:
Notice the same kind of stability situation as we saw before. To the left of the maximum, more slip - more steering - means more grip. To the right of the maximum, more slip means less grip. Here's the slope:
We find that the maximum of the original curve, the zero-crossing of the
slope, occurs at = 3.273°
Once we find the maxima, we can create new, non-dimensional quantities by
scaling and
by these values, namely
. These are pure numbers, so they're commensurable. They are unity when
and
have the values of maximum traction in isolation of one another. We can then
write new functions
and
which have their maxima at
s = 1 and a = 1. We seek a vector-valued
function
of s and a whose longitudinal x
component
expresses the longitudinal force component and
whose lateral y component
expresses the lateral force component under combination slip. Build
this up from
and
so that it satisfies the following requirements:
The following table fleshes out requirement 3 for the cases of braking ( < 0 ) or turning left (
< 0 ). The essential idea is that if the magnitude of either parameter
increases, then the magnitudes of the inputs to the old magic formulae must
increase, but honouring the algebraic signs. If a parameter is positive, it
should get more positive as the magnitude of the other parameter increases.
Similarly, if a parameter is negative, it should get more negative as the
magnitude of the other parameter increases.
sgn( ![]() |
sgn( ![]() |
![]() |
![]() |
input to Fx | input to Fy |
---|---|---|---|---|---|
+ | + | increasing | fixed | increasing | increasing |
+ | + | fixed | increasing | increasing | increasing |
+ | - | increasing | fixed | increasing | decreasing |
+ | - | fixed | decreasing | increasing | decreasing |
- | + | decreasing | fixed | decreasing | increasing |
- | + | fixed | increasing | decreasing | increasing |
- | - | decreasing | fixed | decreasing | decreasing |
- | - | fixed | decreasing | decreasing | decreasing |
Without further ado, here's our proposal for the combination magic grip formula:
|
Using as the input, with the appropriate algebraic signs,
satisfies requirements 1. Multiplying the outputs by the ratio of s to
and a to
magically satisfies requirements 2, 3, and 4. There is, in fact, plenty of
freedom in the choice of the outer multiplier: strictly speaking, any power of
the ratios would do for requirements 2 and 4, and some care will be required to
get the signs right for requirement 3. Until we have a good reason to change it,
we'll just go with the ratio straight up, especially since it
automatically gets the signs right. We close this instalment with a plot of the
magnitude
showing the traction circle very clearly:
The stability criteria are visually obvious, here. If the current, commensurable slip values, s and a, are inside the central "cup" region, then increasing either component of slip increases grip. If they're outside, then increasing slip leads to decreasing grip and the driver is in the "deep kimchee" region of the plot.
ERRATA: The Physics of Racing series has been fairly error-free over the years, but I caught three small errors in part 22 whilst going over it for this instalment. The good news is that they did not affect any final results. I defined the WHEEL frame at the wheel hub but later I implied that it is centred at the contact patch (CP). In fact, the frame at the CP is the important one, and we call it TYRE from now on, avoiding the ambiguous "WHEEL". We never actually used the improperly defined WHEEL frame, so, again, final results were not affected. Also, the dimensions for a3 in Part 22 should be N/Degree, not just N, because a3 furnishes the dimensions for B, which always appears in the combination SB, and has dimensions of degrees. Finally, the dimensions for a6 are 1/KN, not KN. |