Hans B. Pacejka
Hans Bastiaan Pacejka (born 1934; died September 2017) was an expert in vehicle system dynamics and particularly in tire dynamics, fields in which his works are now standard references. He was Professor emeritus at Delft University of Technology in Delft, Netherlands.
The Pacejka "Magic Formula" tire models
Pacejka has developed a series of tire design models over the last 20 years. They were named the 'magic formula' because there is no particular physical basis for the structure of the equations chosen, but they fit a wide variety of tire constructions and operating conditions. Each tire is characterized by 10–20 coefficients for each important force that it can produce at the contact patch, typically lateral and longitudinal force, and self-aligning torque, as a best fit between experimental data and the model. These coefficients are then used to generate equations showing how much force is generated for a given vertical load on the tire, camber angle and slip angle.
The Pacejka tire models are widely used in professional vehicle dynamics simulations, and racing car games, as they are reasonably accurate, easy to program, and solve quickly. A problem with Pacejka's model is that when implemented into computer code, it doesn't work for low speeds (from around the pit-entry speed), because a velocity term in the denominator makes the formula diverge. An alternative to Pacejka tire models are brush tire models, which can be analytically derived, although empirical curve fitting is still required for good correlation, and they tend to be less accurate than the MF models.
Solving a model based on the Magic curve with high frequency can also be a problem, determined by how is the input of the Pacejka curve computed. The slipping velocity (difference between the velocity of the car and the velocity of the tire in the contact point) will change very quickly and the model becomes a stiff system (a system, whose eigenvalues differ a lot), which may require special solver.
The general form of the magic formula is:
where b, c, d and e represent fitting constants and R is a force or moment resulting from a slip parameter k.
- Pacejka, H. B., The wheel shimmy phenomenom: A theoretical and experimental investigation with particular reference to the nonlinear problem (Analysis of shimmy in pneumatic tires due to lateral flexibility for stationary and nonstationary conditions), Ph.D. Thesis, Delft University of Technology, Delft, 1966.
- Bakker, E. ; Nyborg, L. ; Hans B. Pacejka Tyre modelling for use in vehicle dynamics studies 1987 Jan. Society of Automotive Engineers, Warrendale, PA.
- PWA Zegelaar, HB Pacejka The in-plane dynamics of tires on uneven roads[permanent dead link] Vehicle System Dynamics, 1996 - ocp.tudelft.nl
- Pacejka, H. B. Tire and Vehicle Dynamics, Butterworth-Heinemann, Oxford, 2002.
- "Tribute: Hans Pacejka 1934-2017". Tire Technology International. September 19, 2017. Retrieved 1 October 2017.
- Ruud van Gaal (December 23, 2010). "Pacejka's Magic Formula". Racer. Retrieved 2011-03-19.
- Cossalter, Vittore (2006). Motorcycle Dynamics (2nd ed.). Lulu.com. p. 39. ISBN 978-1-4303-0861-4.
- Pacejka, Hans B. (2006). Tyre and vehicle dynamics (2nd ed.). SAE International. pp. back cover. ISBN 978-0-7680-1702-1.
- Brian Beckman (2001). "The Physics of Racing, Part 21: The Magic Formula: Longitudinal Version". Retrieved 2011-03-26.
- http://www.tut.fi/plastics/tyreschool/moduulit/moduuli_10/hypertext/index.html MF is "easy to handle, accurate, low effort"
- Beckman, Brian (2007) Brian Beckman: The Physics in Games - Real-Time Simulation Explained, at channel9.msdn.com, Jun 08, 2007, min. 29:53-33:45 quotation:
- Jacob Svendenius, Björn Wittenmark (2003) Brush Tire Model with increased Flexibility, In European Control Conference, September 2003.
- "Modelling of tyres".
- http://www.tut.fi/plastics/tyreschool/moduulit/moduuli_10/hypertext/index.html "These models are not expected to give a very accurate correspondence with measurement but should predict the qualitative trends."
A typical video game runs its physics engine at 50 Hz (50 calculation cycles per second). The slip formulas require a very high calculation rate (~500 – 700 Hz) in order to provide acceptable results. One of the causes is that calculating ω involves the wheel’s inertia and the torque from the vehicle’s engine. This generates rapidly changing values which cause huge slip ratios at low calculation rates.