Facts and myths on the Pacejka curves

When developing vehicle physics the tire is the single component that has the most influence in the entire system. When designing a tire simulation the Pacejka Curves arise sooner or later.

What’s “Pacejka” exactly?

First of all, the Pacejka Magic formula, also called simply Magic Formula “MF” or “MF-tyre”, is a function y=f(x) which is used to predict and simulate the forces developed by a tire. Given an input x and a list of coefficients (b0..b12 in the picture), the function calculates the force y developed by the tire. Three different functions are used, one for longitudinal forces, other for lateral forces, and the last one for the self-aligning torque (the feedback force you experiment on your car’s steering wheel). Each function has its own set of coefficients.

The coefficients affect the final shape of the curve in a wide variety of ways (see the pic). The goal is fitting the curve with the results obtained in the empiric experiments on the specific tire the curve is being calculated for. Once the proper coefficients for the curve are found, the behavior of that tire can be easy and realistically predicted without having to actually use the real tire. One can imagine using a tire testing facility, testing a specific tire in a wide variety of conditions and measuring the results. Then, based on those results, one could start toying with the coefficients of the Pacekja curves until getting a curve that acceptably fits the experimental results.

What about the coefficient sets for Pacejka curves?

Pacejka and other models are developed by and for the automotive industry. Real coefficient sets are scarce and hard to find, as they are heavily protected intellectual data. Tire manufacturers sell the coefficient sets to interested customers who can then simulate how a specific tire would behave in a vehicle (the MSC Adams software is an industry-standard for modelling this kind of simulations). A Pacejka coefficient set for a single tire model can easily cost around $1000-1500.

You may find many Pacejka sets over there (for instance as part of the car add-ons for the Racer.nl simulator). But most of them are not real Pacejka sets as they haven’t been calculated based on measures on the real tires. Anyone can use a Pacejka editor and get an arbitrary set of coefficients. For that, the Pacejka formula is a realistic and flexible method for simulating an imaginary tire based on the behavior expected in a simulated vehicle.

I’ve written a document which describes how each Pacejka coefficient affects the resulting curve. This is useful for designing tire friction curves matching a specific behavior in vehicle simulation games.

Is Pacejka a necessary or best method for simulating tires?

Surprisingly, the curves themselves are not the most important part of a tire simulation. Remember, a Pacejka curve is simply a function y=f(x). Any function with shape similar to Pacejka (peak – asymptote) will work, no matter how it’s calculated. It could result more or less realistic, but even just  a flat, clamped slope will do a decent work.

The critical part for simulating a tire is how are you calculating the parameters you will feed the curves with. The “x” in y=f(x). And here, my friend, is where things get serious.

What are the problems with the Pacejka-based methods?

Pacejka and other models are based on the amount of slip of the tire over the ground. The longitudinal version uses the difference between the tire’s angular velocity and the actual velocity over the ground (slip ratio). The lateral version is based on the angle between the wheel’s heading direction and the true movement direction (slip angle).

However the slip ratio is a tricky concept. Having slip doesn’t mean that the tire is actually sliding over the ground. The slip ratio reflects the fact that the drive wheels compress the tire at the front of the contact patch where the radius of the wheel is smaller, so it effectively travels faster than the ground underneath. A portion of the contact patch is adherent to the surface while the other portion is slipping. This slip ratio generates the force that moves the car. In a car cruising at constant speed the drive wheels will be spinning a bit faster than the non-drive wheels.

There are several different definitions for slip ratio and slip angle available. These ones are a commonly used definition borrowed from the great book from Brian Beckman The Physics Of Racing:

ω    angular velocity
Re  effective radius
V    lineal velocity
Vy  lateral velocity
Vx  longitudinal velocity

The slip ratio σ and slip angle α are the input for the Pacejka formulas, which yield the longitudinal and lateral forces exerted by the tire. Although this kind of methods is known to provide the state-of-the-art tire simulation, they still exhibit severe problems:

  • A quick look at the formulas reveal an evident and inconvenient fact: the slip formulas become numerically unstable when V approaches zero and are not defined at all when the vehicle is stopped (0/0). Other definitions of the slip ratio also suffer of this problem in a similar way, requiring hacks and workarounds at low speeds.
  • The lateral force is based on the slip angle. However the slip angle does not account for any variation with the speed. Intuitively, the forces generated at high speed must be greater than the forces at low speed with the same slip angles. But using these formulas a velocity of V = (0.2, 0.1) = 0.22 m/s, results in the same slip angle (and thus the same lateral force) as V = (20, 10) = 22.36 m/s.

Is there a better method?

Of course, there are always better methods. The challenge here is to design a method that is fully coherent yet realistic enough that works under all conditions of the tire, including the forces at low speeds and with the vehicle stopped.

This is what I’m currently working on 😉



EDIT (March 2015):

…and finally I’ve designed and implemented a fully coherent and accurate tire model and vehicle simulation that resolves all issues mentioned in this article. The tire and vehicle simulation react realistically to all situations with perfect continuity between all states.

https://www.edy.es/dev/2014/11/early-preview-of-vehicle-physics-pro/

http://vehiclephysics.com