Downloaded from SAE International by Univ of California Berkeley, Thursday, July 26, 2018

~.A
~
e
The Engineering Society
For Advancing Mobility
~ Land Sea Air and Space 400 COMMONWEALTH DRIVE WARRENDALE, PA 15096

hni
a r rl s
870421

Tyre Modelling for Use In Vehicle

Dynamics Studies
Egbert Bakker and Lars Nyborg
Chassis Engineering
Volvo Car Corp.
Hans B. Pacejka
Vehicle Research Laboratory
Delft University of Technology

International Congress and Exposition

Detroit, Michigan
February 23-27. 1987
 Downloaded from SAE International by Univ of California Berkeley, Thursday, July 26, 2018

The appearance of the code at the bottom of the first page of this paper indicates
SAE's consent that copies of the paper may be made for persona! or intern;)l use, cr
for the personal or internal use of specific clients. This consent is given on the con·
dition, however, that the copier pay the stated per article copy fee through the
Copyright Clearance Center, Inc., Operations Center, P.O, Box 765, Schenectady,
N.Y. 12301, for copying beyond that permitted by Sections 107 or 108 of the U.S.
Copyright Law. 'This consent does not extend to other kinds of copying such as
copying for general distribution, for advertising or promotional purposes, for
creating new collective works, or for resale.

Papers published prior to 1978 may also be copied at a per paper tee of $2.50
under the above stated conditions.

SAE routinely stocks printed papers for a period of three years following date of
publication. Direct your orders to SAE Ordp.r Department.

To obtain Quantity reprint rates, permission to reprint a technical paper or per~

mission to use copyrighted SAE publications in other works, contact the SAE
Publications Division.

All SAF papal a", ahllTt1rl..v am} muuw

;" rh( SAF GI"hal ,\(,,/>;11/" /Jala/>Ql(

SAE GLOBAL MOBILITY DATABASE

No part of this pUblication may be reproduced in any form, in

an electronic retrieval system or otherwise, without the prior
written permission of the publisher.

ISSNOI48· 7191
Copyright © 1987 Society of Automotive Engineers,lnc.

This paper is subject to revision, Statements and opinions ad· Persons wishing to submit papers to be considered for pre-
vanced in papers or discussion are the author's and arc his sentation or publication through SAE should send the manu-
responsibility, not SAE's; however, the paper has been edited script or a 300 word abstract of a proposed manuscript to:
by SAE for uniform styling and fonnat. Discussion will be Secretary, Engineering Activity Board, SAE.
printed with the paper if it is published in SAE Transactions.
For permission to publish this paper in fuU or in part, contact Printed in U.S.A.
the SAE Publications Division,
 Downloaded from SAE International by Univ of California Berkeley, Thursday, July 26, 2018

870421

Tyre Modelling for Use in Vehicle

Dynamics Studies
Egbert Bakker and Lars Nyborg
Chassis Engineering
Volvo Car Corp.
Hans B. Pacejka
Vehicle Research Laboratory
Delft University of Technology

ABSTRACT first results of this project have been presented in

this paper.
A new way of representing tyre data obtained The description and understanding of tyre
from measurements in pure cornering and pure behaviour should be easily obtainable from
braking conditions has been developed in order to measured data. Furthermore, it is desirable that
further improve the Dynamic Safety of vehicles. its parameters characterize in some way the typi-
The method makes use of a formula with coeff- fying quantities of the tyre such as slip stiffnesses
icients which describe some of the typifying at zero slip and peak values of forces and torques.
quantities of a tyre, such as slip stiffnesses at zero This feature would make it possible to investigate
slip and force and torque peak values. The formula the effect of changes of these quantities upon the
is capable of describing the characteristics of side handling and stability properties of the vehicle.
force, brake force and self aligning torque with That is, one might then easily model a hypotheti-
great accuracy. This mathematical representation cal tyre.
is limited to steady-state conditions during either
pure cornering or pure braking and forms the REPRESENTATION OF TYRE DATA
basis for a model describing tyre behaviour during
combined braking and cornering. Three possible ways to represent measured
tyre data are in use:
a) representation by tables,
b) representation by graphs,
SAFETY IS A FIELD which has been given top c) representation by formulae.
priority in developing Volvo. cars, because people The first two possibilities are difficult to handle in
should be able to travel safely and securely. For theoretical studies and therefore they do not fulfil
this reason, the vehicle properties have to be the requirements mentioned earlier.
adapted to the driver, so the vehicle will be Two possibilities may be distinguished for the
controllable and its behaviour predictable in· all representation by formulae:
situations. This is called Dynamic Safety and is a) formulae containing series (Fourier,
obtained after optimizing stability, steering and polynomials),
brake performance by tuning chassis and tyre b) formulae containing special functions.
design to each other. Here, the tyre plays a vital The use of series has some disadvantages:
role and thorough knqwledge of its properties is - relatively many coefficients are needed to get a
required. It is therefore necessary to have a proper close curve fit,
description oftyre behaviour available. In order to the resulting curves have a wavy look, which
achieve this, a joint project between Volvo and the implies that the variation of the slope along the
Delft University of Technology was set up. The curve differs considerably from that of the

0148·7191/87/0223·0421 $02.50
Copyright 1987 Society of Automotive Engineers, Inc.
 Downloaded from SAE International by Univ of California Berkeley, Thursday, July 26, 2018

2 870421

measured data. Y = Dsin(BX) (1)

extrapolation beyond the fitted range often
yields large deviations. with Y standing for either side force. self aligning
in general. the coefficients do not describe the torque or brake force and X denoting slip angle (0)
typifying quantities in a recognisable way. or longitudinal slip (Kl. The longitudinal slip is
which makes it impossible to change these defined as the ratio of the difference between the
quantities in a controlled. simple way. speed of rotation of the dri ven or braked tyre and
A representation with the aid of polynomials of the straight free rolling tyre. and the speed of
is described by Sitchin [I]*' This investigator rotation of the straight free rolling tyre. expressed
attempts to solve some of the problems connected as a percentage. A negative value results from a
with the use of polynomials by dividing the data braking torque. In Eq. (I) D is the peak value and
into five regions and using a variable format the product DB equals the slip stiffness at zero
regression equation for each of them. The first two slip.
disadvantages could be solved in this way. but the
others still remain. A problem appears to show up
for the derivative. which is discontinuous at the Side loren
boundaries of the regions.
The best way to fulfil the requirements is to 8rllkn fOlce

find a special function. which through its partic-

ular structure is capable of describing the meas-
ured data with great accuracy and which has Self nligning lo.que

parameters related to the typifying quantities in a

simple manner.
The formula should. if possible. be able to
describe:
the side force as a function of slip angle.
the brake force as a function of longitudinal
slip.
the self aligning torque as a function of slip
angle.
The side force and self aligning torque are Fig 1. - Steady-state tyre characteristics.
measured during pure cornering (cornering with-
out braking) and the brake force durin g pure
braking (braking without cornering). Both corner- Eq. (1) does not give a good representation for
ing and braking are in steady-state conditions. larger values of X. A gradually increasing
A possible special function used by Ruf [2] is extension of the X axis appears to be necessary. To
a combination of tanh functions. The resulting accomplish this. the arctan function has been
formula can match side and brake force charac- used. The formula (1) now changes into:
teristics rather well. However. the coefficients do
not describe the typifying quantities of the tyre in Y = D sin (C arctan (BX)) (2)
an easy way and therefore these quantities are
difficult to recognize and to vary. In Eq. (2) D is still the peak value. the slip
stiffness at zero slip is now equal to the product
PROPOSED TYRE FORMULA BCD (from now on called the stiffness). The
coefficient C governs the shape of the curve. For
The basic form of each of the characteristics of large values ofX. Eq. (2) reduces to:
the tyre (cf. Fig.I) suggest the use of the sine
function as a first step in developing the final Y = D sin (!nC) (3)
formula as summarized in the Appendix.
Consequently. C defines the extent of the sine that
will be used and therefore determines the shape of
• Numbers in brackets [ 1designate references at the curve. The value of C makes the curve look
the end of the paper. like a side force. brake force or a self aligning
 Downloaded from SAE International by Univ of California Berkeley, Thursday, July 26, 2018

870421 3

torque characteristic (we find for C in those cases icients, which is able to describe all the measured
successively: 1.3, 1.65 and 2.4). With C determined characteristics. The four coefficients are:
by the shape and D determined by the peak value,
only B is left to control the stiffness. B= stiffness factor
Still Eq. (2) is not good enough to describe C= shape factor
every possible measured characteristic. There may D= peak factor
be a need for an additional coefficient which E= curvature factor
makes it possible to accomplish a local extra
stretch or compression of the curve. The coefficient So far, the characteristics are assumed to pass
E has been introduced into the formula in such a through the origin. In reality, however, this will
way that stiffness and peak value remain not always be the case. Due to ply steer, conicity
unaffected. and rolling resistance, the characteristics will be
shifted in the horizontal and/or vertical directions
Y = D sin (C arctan (B¢)) (Fig. 2), In order to fit the measured
with (4) characteristics, these shifts should be included in
¢ = (I-E) X + (EfB) arctan (BX) the equation. We obtain:

The influence of E on the side force characteristic Y = D sin (C arctan (B¢)) + Sy (5)
has been shown in Fig. 2. Similar effects occur with
with the self aligning torque and brake force ¢ = (I-E) (X + Sh) + (EfB) arctan (B (X + Sh))
characteristics. Sh = horizontal shift
The result is an equation with four coeff- Sy = vertical shift

Side force I
-
Fy [NI I
6000 - I Sh
-
T - -.. ;"-;;'~':::'---- --
I
4000 - ' .",-

I
.' /""
/~/
I
0
2000 I
I
-I
BCD
1-" 10 "5 L...
-15 -10 -5 5
° I Slip angle a [deg.J
5v
-2 -I
I
-I
/....::"/~OOO I
__
," ..-0' I
------:: .............
-6000 - I
I _EeoO
__ O<E< 1
•••• E < 0
I
I
I

Fig. 2 - Coefficients appearing in tyre formula.

 Downloaded from SAE International by Univ of California Berkeley, Thursday, July 26, 2018

4 870421

THE MEASUREMENTS - Full scale tyre quadratical correction curves around the meas-
measurements have been carried out on a dry ured points. Finally, the average was taken of two
asphalt road. The experiments were conducted or more a-sweeps.
with a specially built test trailer under steady- The measurements during pure braking were
state conditions. The measured quantities are: carried out by gradually increasing the brake
the side force (Fy), pressure until wheel lock occurs. The data was
- the brake force (Fx), collected while the longitudinal slip was increased
- the self aligning torque (M z)' from 0 to 100% (K-sweep). Here also a
These forces and torque were measured during condensation, correction and averaging process
pure cornering, pure braking or during combin- has been carried out. The experiments on
ations of cornering and braking. During cornering, combined cornering and braking were conducted
the speed was 70 km/h and during braking 60 by doing K-sweeps at fixed values of the slip angle.
krnlh. The tests included variations of slip angle The camber angle was fixed at zero degrees.
(-10 to 14 deg.), vertical load (2, 4, 6, 8 kN) and FITTING MEASURED DATA- The data
camber angle (-5, 0, 5 deg.l. obtained after processing in the way described
The pure cornering tests were carried out by above has been used as raw input data for the
slowly varying the slip angle (a-sweep) from zero fitting process. An optimization program from the
to the maximum value, from there to the NAG library has been used to obtain the four
minimum value and back to zero again. The coefficients and two shifts for every measured
measured data was condensed per 1/4 degree characteristic. For one of the measured tyres, the
intervals to one measured point per interval. After resulting coefficients have been tabulated in Table
that, corrections were made for load variations, 1. .The vertical shift obtained for the brake forces
which occurred during the tests, by making use of represents the value for the rolling resistance.

Table 1
Coefficients for tyre formula (first fit)
(a [deg.], K [%])
LOAD
B C D E Sh Sv BCD
(kNl

2 0.244 1.50 1936 -0.132 -0.280 -118 780.6

Fy 4 0.239 1.19 3650 -0.678 -0.049 -156 1038
6 0.164 1.27 5237 -1.61 -0.126 -181 1091
8 0.112 1.36 6677 -2.16 0.125 -240 1017
2 0.247 2.56 -15.53 -3.92 -0.464 -12.5 -9.820
Mz 4 0.234 2.68 -48.56 -0.46 -0.082 -11.7 -30.45
6 0.164 2.46 -112.5 -2.04 -0.125 -6.00 -45.39
8 0.127 ·2.41 -191.3 -3.21 0.009 -4.22 -58.55
2 0.178 1.55 2193 0.432 0.000 25.0 605.0
Fx 4 0.171 1.69 4236 0.619 0.000 70.6 1224
6 0.210 1.67 6090 0.686 0.000 80.1 2136
8 0.214 1.78 7711 0.783 0.000 104 2937
 Downloaded from SAE International by Univ of California Berkeley, Thursday, July 26, 2018

870421 5

SIde lo.eo Sel! aligning torquo

Fy IN] M~ INm]
0000 150·

10 15
Slip \logla a ldog.] "o Idog.]
S lip ongle

- 100 ~

*F, 6000N
-150 .
.. M."ulom."I'

:-!!L.'':_~

Fig. 3 - Side force characteristic fitted using the Fig. 4 - Self aligning torque characteristic fitted
tyre formula compared with measured using the tyre formula compared with
data. measured data.

Blah lorco
-F>lIN]

oooo:~
- -................J

'000

-~_ ... ----- ---+---_.__._-- ..

-25 0 50 75
- Longiluditllll ~lip
-Jl (%J

Mouutomoni.
!,., F, _ 6000N i

Fig. 5 - Brake force characteristic fitted using the

tyre formula compared with measured
data.
 Downloaded from SAE International by Univ of California Berkeley, Thursday, July 26, 2018

6 870421

Figures 3, 4 and 5 show comparisons between (aligning stiffness) characteristics, we use the
measurements (after processing) and fitted approximation:
characteristics for one specific load. The result is
very satisfactory: the formula is able to describe a3 F z' + a4 F z
all the measured characteristics. BCD=---- (8)
At the moment, the shifts are not of particular
interest for further examination and are therefore
put equal to zero (Le. ply steer, conicity and roll- The shape factor (e) appears to be practically
ing resistance will not be taken into account). If independent ofF z. We take for:
so required, it is possible to add the shifts later on. the side force C= 1.30
INFLUENCE OF VERTICAL LOAD - To the brake force C= 1.65 (9)
reduce the total number of quantified coefficients the self aligning torque: C= 2.40
(Table 1) and to be able to calculate forces and
torques at vertical loads which are different from The stiffness factor (B) is found by dividing the
the values used in the measurements, it is necess- stiffness by the shape and the peak factor.
ary to include the vertical load explicitly in the
formula. To do so, the coefficients have to be B = BCD/CD (IO)
written as a function of the vertical load (F z).
The peak factor (D) as a function ofF z may be Finally, the curvature factor (E) as a function of F z
approximately represented by the relationship: is given by:

(6) (I 1)

For the stiffness (BCD) of the side force character- After processing the raw data again (corrected
istic (cornering stiffness), we take the formula: for the earlier determined horizontal and vertical
shifts) with the aid of the optimization program,
the new coefficients are obtained. Table 2 gives the
coefficients for the same tyre as in Table 1, but
and for the stiffness of both brake force (longitud- now incl uding the infl uence of the vertical load.
inal slip stiffness) and self aligning torque

Table 2
Coefficients for tyre formula (with load influence)
(Fz[kNJ)

al a2 a3 a4 a5 as a7 as
Fy -22.1 1011 1078 1.82 0.208 0.000 -0.354 0.707
Mz -2.72 -2.28 -1.86 -2.73 0.110 -0.070 0.643 -4.04
Fx -21.3 1144 49.6 226 0.069 -0.006 0.056 0.486

. Table 3
Coefficients for tyre formula connected with camber influence
(y [deg.])

a9 aID all al2 a13

Fy 0.028 0.000 14.8 0.022 0.000
Mz 0.015 -0.066 0.945 0.030 0.070
 Downloaded from SAE International by Univ of California Berkeley, Thursday, July 26, 2018
870421 7

Side l(lrc~ Sell al'9"ing torque

Fy IN) O·O-.v O ·1t>4 •••••••••••••• M~ INml
0000 f\'~
....;..-------1
.~ _F,
••• F,

__ F,
8000 tJ
fiOCXlN
':OOON
._F, 2000 'j

10 Hi .'
Slip angln Q Ideg.1
...........
l "" Ideg.1
S lip "ngle

. .>
/
I - - - - -............~... / O~,

*X F,
fiOOOll
~OOON
••• f,
-"
BOOO N
Boool<
'00

/
/' ., ""
x "
8000 rl
6000 rl
4000 r,
F, 4000r~ __ F, 4000 N 150 + F, 2000 r,
.......................0."4>·0 + f, 2000r, ._ F, 20001< ,:;'
i

Fig. 6 - Side force characteristics fitted using the Fig. 7 - Self aligning torque characteristics fitted
tyre formula which includes the influence using the tyre formula which includes the
of the vertical load, compared with meas- influence of the vertical load, compared
ured data. with measured data.

l"uko lo.eo
-F~ INI
t't.".III."

........,••••I"••••" ••• jS .

• •

SO 75
~ longitudinal slip
-~ 1%1

""'" MN,u,.",.n,.
t. "10
•, " BOOO
6000 N
.••• F,
-',
BlXX} N
6000 'I
" 4000 ••
__ F,
4000 'I
""'" + ""
X
2000 II ._F, 2000 N

Fig. 8 - Brake force characteristics fitted using the

tyre formula which includes the influence
of the vertical load, compared with meas-
ureddata.
 Downloaded from SAE International by Univ of California Berkeley, Thursday, July 26, 2018

8 870421

dFy/dp dFy/da
IN/dog,1 "00 IN/dllg·1 I,•••• F, 80001;1
_f, 600011
-_F, .:000 /j
1200 -.
!~:;~F;--~ 6OOOrJ! "00 ~_F: 2000 N
_F, 6000~l
__ F, .10CXl rl
0.:;.• ••••••••••
,.-F, 2000N I "'00
"'\.....

'00
\.
\\,
600 - . . . . "

'00
\
... \
........... '00 \. \ '.
, \. \_\ \...~
"'00 4000 6000

-200 ~
Slip l,lngle (>
"
[dll!!.)
Sido IOlee
Fy IN]

Fig. 9 - Calculated slope of side force characteris- Fig. 10 - Calculated slope of side force characteris-
tics as a function of slip angle for different tics as a function of side force for different
vertical loads. vertical loads.

Pneumatil; trail Pneumatic lIail

l{m] tlml
O,OB ~ OOB ' •••• F, BOOON'
_F, 6000N
--F: .:OOOtl
r.... F:~~'-BOOON I '-F, :1000lj
_F.·6000 N "
__ F, _ 4000 N

//.......•..•............\.
OOt; -
,=--:-F, ~ 2000~_ '" ...•........................................

,""'1-----_.

••••

1-;;;;;-;;;:;;
t::: ....~~~
/I.~·' \\\
••••••

;?>-- I'
,,/

",/./
# 4

I
' , ' ,

"_ '........
••••••

••••

"'-''''''-_~'..;;;.''"",,-~~
'"

, \ ) \
_5 0 5 10
Slip ongJo" [deg.]
-'. J 4000 6000 .J Side lorce
Fy INI
1

Fig. 11 - Calculated pneumatic trail as a function Fig. 12· Calculated pneumatic trail as a function
of slip angle for different loads. of side force for different loads.
 Downloaded from SAE International by Univ of California Berkeley, Thursday, July 26, 2018

870421 9

The new parameter a3 for the side force indicates pure cornering and braking conditions, will be
the maximum possible cornering stiffness gener- used now in a mode which describes the beha viour
ated by the tyre. of the tyre at zero camber angle when a combin-
A comparison between measurements and the ation of cornering and braking takes place. The
newly obtained characteristics has been shown in modelling process implies the use of the basic tyre
Figs. 6,7 and 8. The influence of the vertical load is modelling theory as described in [31.
now included without much loss in accuracy. Figure 13 shows a top view of a tyre during
INFLUENCE OF CAMBER -The total num- combined cornering and braking. The force and
ber of coefficients to be quantified may be reduced torque acting on the tyre and the velocity of the
even more by including camber in the formula as wheel have been indicated:
well. The influence of camber (y) is incluced by the total force with components in lateral and
adding two shifts and a possible change in stiff- longitudinal direction (Fy and F x),
ness. The two additional shifts are: the self aligning torque (M z),
the speed in the direction of tra vel V which is
composed of the slip speed vector (V s ) and the
(2) rolling speed vector (Vr) of the tyre.

The change in stiffness is obtained by multi-

plying B by O-alZlyl). We have for the increment
ofB:

(3)
-
~IV
1-

The values of the self aligning torque at

higher slip angles will change due to this change
in stiffness. To compensate for this effect, the
curvature factor E for Mz must be divided by F- Fy
(l-a13Iyl). For the side force, this compensation is
not necessary because the values at higher slip
I
angles are not changed appreciably by changing Fig. 13 - Top view of tyre during combined
the stiffness. The values obtained have been braking and cornering showing velocity
shown in Table 3. and force vector diagrams.
At this stage about five thousand measure-
ment points have been reduced to three equations BRAKE AND SIDE FORCE DURING COM-
with a total of thirty-one coefficients (coefficients BINED CORNERING AND BRAKING - If the
with value zero not included). It is now possible to tyre is assumed to have isotropic properties then
calculate forces and torques which would arise in the direction of the total force is opposite the slip
circumstances not included in the measurement speed vector. The length of the total force vector is
programme. The equations are summarized in the a function of the total newly introduced slip quan-
Appendix. tity a [31. This theoretical slip has components
Beside side force, self aligning torque and which are defined as follows:
brake force, other information may be obtajned
with the aid of the formula. For instance dFy!da as ax = Vsx/V r
a function of slip angle or as a function of side force ay = Vsy!V r (4)
(Figs. 9 and 10). It is also possible to calculate the
pneumatic trail as a function of slip angle or side
force, merely by dividing self aligning torque by
the side force (Figs. 11 and 12). The vector a has the same direction as the slip
speed vector Vs. Therefore we find:
MODEL FOR COMBINED CORNERING AND
BRAKING F =
-(ii!a)F(a)
F x = -(ax/a)F(a) (5)
The formulae obtained which describe the F y = -(ay!a)F(a)
 Downloaded from SAE International by Univ of California Berkeley, Thursday, July 26, 2018

10 870421

The relations between the theoretical slip and the

F~o
practical slip quantities used earlier (0 and K) are: Fyo .
INI

,,,..-........." " __ F,,,

It'"
(16) I ............ - F yo

,, ,.- -....
= VsylVr = -tana/(I + K) "
6000 fl
ay
I I .......
------ -------------
where we had: II ' - _....
JI
/I
1/ V --_....
----- "
---------
llOOO N

(17) ""
II
II
--- ....
II I
------------------ "
.1000 tl

"""/
'
As stated before, the theory is valid for a tyre
with isotropic properties (the same characteristic
.u'
I'
I'
,I' /""--------- --------- - - " 2000 II

for side and brake force). But in reality a tyre is

not isotropic, so the theory must be adjusted.
"V
1/1
,
II

Instead of a single characteristic F(a), two

different characteristics should be used (cf. Lugner o 01 02 03 04 05 0& 07 OS 09 10

[ 4]): Theoretical slip 0

- for pure cornering: Fyo (a)

- for pure braking: Fxo (a) Fig. 14 - Original calculated basic characteristics.
These characteristics are called the original basic
curves. The total force may now be derived as
follows:

Fx = -(axla) F xo (a)
Fy = -(ayla) Fyo (a) (18)
. . . . -r-- . . __ F.c>
""'" ,---:::=F=:::::---------------,
F = YF x' + Fy' ,I' I --____ _ f
yo

! -L --------_
1 -----_
The original basic curves are obtained by calcula-
ting Fx and Fy in pure conditions and writing them 6000 /
1
/ ........ -T--------___ f,
----
1I000N

as a function of the theoretical slip with the aid of "

1,1 / I
--------------------
Eq. (16). Figure 14 shows some of these character-
istics for different vertical loads.
11'/
II _---
I
_
f, 6000/1

I --------
A problem of physical nature arises when the
values of the slip connected with the peak values
4000

"f(
I'
I'
"" I /"
I
----- _
-----
f, 4000!1
of F xo and Fyo differ considerably from each other. ,1
~
"
At a slip a in between these two peaks, we have the
situation that total slip would occur on one curve 2000: . . . . ------------------- _
".Y
, '/I I
and partial sliding on the other. In reality, the tyre I

has one total condition of sliding. To solve this

I
problem a new slip quantity is defined which ,,
normalizes the theoretical slip. This new slip is ,
called the normalized slip a*' The normalized slip No.mnli~ed slip o'

is obtained by dividing a by the slip value which

occurs at the peak of the characteristic. We desig- Fig. 15 - Normalized basic characteristics.
nate these values of slip as a xm and a ym for the Fx
and the F y characteristics respectively. The
components of the normalized slip read:
 Downloaded from SAE International by Univ of California Berkeley, Thursday, July 26, 2018

11
870421

0*

1
----------------------------------
!t
- - - - - - - - - 1"- I
I
I Ed
I
I
I
I
I
-F x
I
1 I ' ..0*
I
I
I
I
\ \ I
\',. ."..... . . . . Fyo
I
!
, '- ---
_- _
------~_..!-_-----
........
.

Fy
Fig. 16 - Visualization of determination of magnitude and direction of the total horizontal force Fusing
the normalized basic characteristics, starting out from a given normalized slip vector (cf. Figs. 17,
18 for the factor Cd).

Direction focla, DircClion '.:Iclor

" " 2 dog "

F, 6000 N

F: 6000 N

FI ,4000N

F, 2000N
~.:::
F, " 2000N

, L ~ ~ ~ ----J ,t ~ ~ --' ...J

100 1000
" " 100
Normllliled slip a·
1000
" Normalized slip,,'

Fig. 17 - Direction factor as a function of 0* as Fig. 18 - Direction factor as a function of 0* as

derived from measured data for a slip derived from measured data for a slip
angle of 2 degrees. angle of 5 degrees.
 Downloaded from SAE International by Univ of California Berkeley, Thursday, July 26, 2018

12 870421

0,,* = Ox/Oxm tion factor (Cd), when introduced, yields:

Oy* = oy/oym (19)
Fx = -(ax */0*) Fxo (0*) (25)
Fy = -Cd (Oy */0*) Fyo (0*) (26)

The basic curves, but now as a function of 0* f'igure 16 shows the contruction of the total
have been shown in Fig. 15. On both characteris- horizontal force. The magnitude of Cd will be 1 for
tics (Fxo and Fyo), partial sliding occurs for values small values of 0* and Oym/oxm for large values of
of 0* smaller than 1 and total sliding for values 0*. For intermediate values of a* the value for Cd is
larger than 1. With 0* as the new slip the determined from the measured relationships
equations (18) become: during combined cornering and braking. Figures
17 and 18 show e"amples of the variation of Cd
Fx = -(ax */0*) F xo (0*) (o*,F z, a) for different loads and slip angles.
F y = -(Oy */0*) Fyo (0*) (20) Because of the action of the brake force, the
resulting cornering stiffness of the tyre appears to
F = YF/ + Fy" become larger. For this reason, the basic normal-
ized side force characteristic has been given an
The direction of the total force is now based on initial slope which is a function of the brake force.
the direction of a'. For small values of 0*, Eq. (20) This can easily be achieved by adjusting the
is correct since then the forces will be independent coefficient B for the side force. The increment of B
and the linear relations hold: is:

F xo = CPK 0xm 0* LlB = -PFxB (P>O) (27)

Fx = -(ax*/0*) Fxo (0*) (21)
= -(ax*/o*)CFx oxmo* with P denoting the longitudinal force sensitivity
= -CFK ax of the cornering stiffness.
Further study is required to explain the
and for Fy after a similar calculation: physical background of P and of the course of Cd,
specially at high loads.
(22) Figures 19 until 22 show comparisons
between measurements and calculations for a tyre
with CFK and CFQ the longitudinal slip and the in various situations during combined cornering
cornering stiffness respectively. and braking. The graphs depict the variations of
However, for large values of 0* and equal the brake force and the side force as a function of
friction coefficient (11) in " and y direction, the longitudinal slip and also the side force as a
total force should be in the opposite direction of the function of the brake force. To be able to compare
slip speed vector (and thus 0). calculations and measurements, the measured
brake forces were corrected for the rolling
Fxo = Il F z resistance (vertical shift). With the aid of Cd and
Fx = -(oxla)Fxo (0*) (23) adjustment ofB (side force), the model gives a good
= -(ox/o)Il F z description of the measurements.
SELF ALIGNING TORQUE DURING COM-
and similarly for Fy: BINED CORNERING AND BRAKING - To obtain
the self aligning torque during combined
(24) cornering and braking, the side force has first to be
multiplied by the pneumatic trail (t). The
So with small slip, 0* should be used to deter- pneumatic trail has been derived as a function of a
mine the direction of the total horizontal force and (or tan a) during pure cornering. The deflection
with large slip, a should be used. Obviously, a need distribution over the contact length during
e"ists for an additional factor which controls the combined cornering and braking is equal to the
direction of the total horizontal force. This direc- distribution during pure cornering if tan a = a (cf.
 Downloaded from SAE International by Univ of California Berkeley, Thursday, July 26, 2018

870421 13

.= ,---="-.,,----=...,,,,-----------------, ,= ,--------------------------,
Sidtl 10leo F IN]y F, 800011
F, 6000 II Brake fOlco -Fl( [NI F, 6000 t<
~- F, £000" f, <100011
3 F, 4000'1 F, 200011
4_ f: 2000"

.................! .

......

-.•.............
. .
'. '.
....••.......•.......•...•........

~.
............•.......•.
..........•....................:.:::::.:.::::.::.:::
...................
'.
10 110 90 100 o 60 70 80 90 100
-Longitudinal 511p -~ 1%1 -longitudinal 511p -.0;:]%1

Fig. 19 - Side force and brake force vs longitudinal Fig. 20 - Side force and brake force vs longitudinal
slip at 2 degrees slip angle. slip at 5 degrees slip angle .

•= , - - - - - - - - - - - - - - - - - - - - - - ,
Sid~ lorce
.=,-----------------------,
Side force
Fy IN] Fy [NI

,=
~ I--------'-'-"~=-'~:~. ~. ~. ~. "...... ..\

l-===::::::::::::~" -;:;:; ) 1----..

•...•.••.......
!

,= :--.~..~..~..~..-..-.----,,-"~.." ..../ / /
'=r-....~
..-
..........

o 6000 8000 6000 6000

Blllkc fOlce -FJI. INI Brnke fOlce -F.INj

Fig. 21· Side force vs brake force at 2 degrees slip Fig. 22 - Side force vs brake force at 5 degrees slip
angle. angle.
 Downloaded from SAE International by Univ of California Berkeley, Thursday, July 26, 2018
14 870421

[3]). Consequently instead of t (tan 0) we may

write t (0) which can be converted into t(o*). The
self aligning torque now becomes:

Mz' = -t(o*) Fy (28)

",--------------------,
51111 aligning to/que

Due to carcass compliance, the lines of action of

M11Nmi ...

>'*.:::·;::~:·~·F.
both the side force and the brake force will be
shifted which causes a contribution to the self , ., ;;;;,;;··iI'fm ..i INI

.,~ -=~."'.~ /:=<.//

aligning torque:

Mz = Mz' - sFxFy (29)

with \-"-'-"""-'-'
1 1 (30)
s = Co; Cc; F, 6000 ,.

_150 I- ..."..."..'_~
~
Ccx and Ccy denoting the longitudinal and lateral
f, _ Ilooo r.
carcass stiffnesses respectively.
Figures 23 and 24 show comparisons between .,,,
calculations and measurements. The measured
characteristics have been corrected for an appar- .250L... -.J
ently constant offset of the line of action of the
brake force of approx. 15 mm. This means that Fig. 23 - Self aligning torque vs brake force at 2
even with a slip angle of zero, a large self aligning degrees slip angle.
torque is present. The overturning couple (M x) was
not measured, sono information is available about
the offset of the vertical force which may have
caused this. The measured values shown in Figs.
23 and 24 have been corrected for this phenome-
non, which is not yet fully understood. ",--------------------,
Sell aligning la/quo
POSSIBLE EXTENSIONS - The model so far M11Nmi

developed, including interaction between side and

brake force, may be extended with camber and Bloke force
-F~ IN]
positive longitudinal slip (driving), for instance, if
measurements in these conditions are available. -50- f, 4000"

The measurements described earlier did not

include these combinations ofinput. . ,,, ....•........

CONCLUSION AND SUMMARY

.", ........••......................................
The proposed formula turns out to be not only
very accurate in describing the measured data. It
also characterizes some of the typifying quantities -200- •••• MN,u'emonl~

5 ~og
of the tyre such as slip stiffnesses and peak values <>

and permits the calculation of forces and torque in _2S0L... ....J

conditions which deviate from those imposed
during the actual measurements. Compared with Fig. 24 - Self aligning torque vs brake force at 5
the conventional way of representing tyre data by degrees slip angle.
polynomials with coefficients which have no direct
physical meaning, this way of representing tyre
data seems to be a big step in the right direction.
The results of the model which describes the
behaviour of the tyre during combined cornering
 Downloaded from SAE International by Univ of California Berkeley, Thursday, July 26, 2018

870421 15

and braking are very similar to the measurements APPENDIX

which were carried out in these conditions. This
means that the representation of the tyre data PROPOSED TYRE FORMULAE
which is valid in pure conditions can also be used
to provide information for a combined condition if SIDE FORCE:
a new slip quantity is used: the normalized slip.
Because of the physical meaningful coefficients Fy '" D sin (C arctan (Bcjl)) + t.S v
of the representation, it is easy to vary some of the with
typifying quantities of the tyre. Simulations may cjl '" (I-E)(a +t.Sh)+ (E/B) arclan(B(a+t.Sh))
then be carried out with a hypothetical tyre. In
this way, the performance of a specific vehicle can D '" alFz' + a2Fz
be predicted. C '" 1.30
The described mathematical representation is a3sin (a4aretan (a5Fz))
just the beginning of an effort to describe tyre B '" ( - - - - - - - - ) (I-a I2Iyl)
behaviour in a way which is accurate, compact, CD
physically meaningful and easy to use. By contin- E '" asF/ +a7Fz+a8
uing this effort, hopefully, the tyre will be given t.Sh '" agy
a sufficiently complete representation in vehicle t.Sv '" (al0F/ + allFz)y
models as accurate as justified by the great
importance the tyre has in controlling the horizon- SELF ALIGNING TORQUE:
tal motion and therefore in the Dynamic Safety of
the vehicle. Mz '" D sin (C arctan (Bcjl)) + t.Sv
with
REFERENCES cjl '" (I-E)(a + t.Sh) + (E/B) arctan(B(a + t.Sh))

1. A. Sitchin, "Acquisition of Transient Tire D '" alFz' + a2Fz

Force and Moment Data for Dynamic Vehicle C '" 2.40
Handling Simulations", SAE Paper No. a3 F z' +a4 F z
831790, November 7-10, 1983 B '" ( ) (I-a12lyl)
CDe asFz
2. G. Ruf, "Theoretische Untersuchungen des E '" (asFz' + a7Fz + a8)/(l-aI3Iyl)
Federungsverhaltens von Vierrad-Strassen- t.Sh '" agy
fahrzeugen", Fortschr.-Ber. VDI-Z Reihe 12 t.Sv '" (alOF/ + allFz}y
No. 44104-105, Dusseldorf 1983
BRAKE FORCE:
3. H.B. Pacejka, "Modelling of the Tyre as a
Vehicle Component with Applications", Lehr- F x '" D sin (C arctan (Bcjl»)
gang V 2.01, Carl-Cranz-Gesellschaft, 1982 with
cjl '" (I-E) K + (EIB) arctan (BK)
4. P. Lugner, R. Lorenz, E. Schindler, "The
Connection of Theoretical Simulation and D '" a[F/ + a2Fz
Experiments in Passenger Car Dynamics", C '" 1.65
Proc. 8th IAVSD Symp. on the dynamics of a3F/ +a4 F z
vehicles, Cambridge, Mass., August 15-19, B '" - - - - -
1983, ed. J.K. Hedrick, Swets & Zeitlinger CDeasFz
B.V., Lisse Neth., 1984 E '" asF/ + a7Fz + a8

The values for the coefficients al until a8 are

tabulated in Table 2 and for ag until a13 in Table 3.
Downloaded from SAE International by Univ of California Berkeley, Thursday, July 26, 2018
 Downloaded from SAE International by Univ of California Berkeley, Thursday, July 26, 2018

This paper is subject to revision. Statements and opinions ad· Persons wishing to submit papers to be considered for pre-
vanced in papers or discussion are the author's and are his sentation or publication through SAl should send the manu-
responsibility, not SAE's; however. the paper has been edited script or a 300 word abstract of a proposed manuscript to:
by SAl for uniform styling and fonnat. Discussion will be Secretary, Fn/;-tinecring Activity Board, SAl:.
printed with the paper if it is published in SAE Transactions.
For permission to publish this paper in full or in part, conta\.·t Printed in U.S.A.
the SAl Publications Division.
Downloaded from SAE International by Univ of California Berkeley, Thursday, July 26, 2018