Tyre Models
Tyre Models
Tyre Models
1
USER MANUAL 2008
Table of Contents
1
OVERVIEW ..................................................................................................................................... 3
1.1
1.2
1.3
1.4
1.5
MODEL USAGE.............................................................................................................................. 7
2.1
2.2
2.3
INTRODUCTION .......................................................................................................................... 3
MF-TYRE .................................................................................................................................. 4
MF-SWIFT................................................................................................................................. 4
NEW FEATURES IN MF-TYRE/MF-SWIFT 6.1 ............................................................................... 6
LICENSING OF MF-TYRE/MF-SWIFT 6.1 ..................................................................................... 6
OVERVIEW .............................................................................................................................. 12
BACKWARD COMPATIBILITY ...................................................................................................... 14
SCALING FACTORS................................................................................................................... 16
PARAMETERS IN THE TYRE PROPERTY FILE................................................................................ 17
ADAMS.................................................................................................................................. 28
MATLAB/SIMULINK/SIMMECHANICS ........................................................................................ 31
LMS DADS ............................................................................................................................ 33
THIRD PARTY SOFTWARE ......................................................................................................... 34
REFERENCES.............................................................................................................................. 35
2/35
1 Overview
1.1
Introduction
The contact interaction between tyres and the road largely affects the driving performance of vehicles.
Automotive engineers are optimising the tyre-road interaction so that the vehicle handles well and
operates both safely and comfortably under any circumstance. To analyse the influence of tyre
properties on the dynamic behaviour of vehicles, the engineer requires an accurate description of the
tyre-road contact phenomena. TNO Delft-Tyre provides a complete chain of tools and services for
detailed assessment and modelling of vehicle-tyre-road interaction.
The tyre models MF-Tyre and MF-Swift can be used in vehicle dynamics simulations in all major
simulation packages to efficiently and accurately represent tyre behaviour for applications ranging
from steady-state to complex high frequency dynamics.
MF-Tyre and MF-Swift contain the latest implementation by Delft-Tyre of Pacejkas renowned Magic
Formula tyre model. With MF-Tyre you can simulate validated steady-state and transient behaviour,
making it a very suitable tyre model for vehicle handling, control prototyping, or rollover analysis. With
MF-Swift you can simulate tyre dynamic behaviour up to about 100 Hz, which is particularly useful for
vehicle comfort, durability, dynamic vehicle control, or vibration analysis.
Special attention has been paid to include behaviour necessary for special applications such as
motorcycles (regular and racing), motorsport (e.g. Formula 1) or aircraft tyres.
TNO Delft-Tyres MF-Tyre and MF-Swift are available for all major simulation packages. TNO DelftTyre makes sure that the tyre model implementation and simulation results are identical and that the
same set of tyre model parameters can be used for all these packages. Further, MF-Tyre and MFSwift are fully compatible with all previous official TNO Delft-Tyre releases.
3/35
1.2
MF-Tyre
MF-Tyre is TNO Delft-Tyres implementation of the world-standard Pacejka Magic Formula tyre
model, including the latest developments by TNO and Prof. Pacejka [1] and [2]. MF-Tyres semiempirical approach enables fast and robust tyre-road contact force and moment simulation for steadystate and transient tyre behaviour. MF-Tyre has been extensively validated using many experiments
and conditions. For a given pneumatic tyre and road condition, the tyre forces and moments due to
slip follow a typical characteristic. These steady-state and transient characteristics can be accurately
approximated by MF-Tyre.
Steady state tyre lateral force as function of longitudinal and lateral slip, calculated using MF-Tyre.
MF-Tyre calculates the forces (Fx, Fy) and moments (Mx, My, Mz) acting on the tyre under pure and
combined slip conditions on arbitrary 3D roads, using longitudinal, lateral and turn slip, wheel
inclination angle (camber) and the vertical force (Fz) as input quantities.
MF-Tyre is valid for large slip angles (typically over 30 degrees), longitudinal slip (100%), large load
variations (including truck tyre loads) and large camber angles (including motorcycle camber angles;
MF-Tyre 6.x includes the functionality of MF-MCTyre). It can handle road undulations that have a
wavelength larger than the tyre circumference and is typically applied for vehicle handling simulation.
1.3
MF-Swift
In addition to the Magic Formula description in the MF-Tyre part of the model, MF-Swift uses a rigid
ring model in which the tyre belt is assumed to behave like a rigid body. This means that the model is
accurate in the frequency range where the bending modes of the tyre belt can be neglected, which,
depending on the tyre type, is up to 60 100 Hz. MF-Swift has been validated using measurements
of a rolling tyre (7 to 40 m/s) containing frequencies up to 120 Hz. The model includes essential
gyroscopic effects.
The tyre model functionality is primarily based on [1] [6]. TNO has made several crucial changes
and enhancements in cooperation with Prof. Pacejka to the models as described in [1] in order to
improve functionality, robustness, calculation times, user-friendliness and compatibility between
various operating modes.
MF-Swift uses an efficient single point contact for slip calculation which results in full compatibility with
MF-Tyre. Due to the introduction of a so-called phase leading network for the pneumatic trail, MFSwift is suitable for path curvature with a wavelength in the order of two times the contact length. For
braking/traction applications, wavelengths as small as half the contact length are well described. The
4/35
transient slip behaviour is well described up to full sliding, due to modelling of decrease in relaxation
length for increased slip levels.
5/35
1.4
With respect to MF-Tyre/MF-Swift 6.0 the following changes have been made:
1.5
Introduction of tyre pressure dependency on the tyre characteristics. This includes the Magic
Formula, tyre stiffness, rolling resistance and other properties.
Replacement of the 2D road contact method using basic functions by the more robust and
accurate ellipse contact method. The ellipse parameters can be used for both 2D and 3D road
contact. Backward compatibility is maintained, so older tyre property files with basic function
parameters will keep on working.
A parameter DRUM_RADIUS has been added to the TNO road surfaces to allow simulations
on a drum surface. The tyre model automatically adjusts tyre properties to account for the
global road curvature.
The licensing system of MF-Tyre and MF-Swift 6.1 depends on the multibody/simulation package in
which it is used and the used operating system. Please read the license manual, license agreement
and terms of use that are supplied with the Delft-Tyre and/or multibody/simulation software. If things
are unclear please contact TNO Automotive (http://www.delft-tyre.nl/). The operating modes that are
supported by MF-Tyre and MF-Swift licenses are discussed in section 2.1.
6/35
2 Model usage
2.1
Operating modes
MF-Tyre/MF-Swift 6.1 is set up in a modular way and allows a user to independently set the operating
mode of the Magic Formula, tyre dynamics and contact method. In some software packages this is
done by defining a four digit value for the parameter ISWTCH in the GUI (DADS); for some other
packages the selections can be made from a menu (e.g. SIMPACK, MATLAB/Simulink). In ADAMS
changes to the operating mode can be made by setting the parameter USE_MODE in the [MODEL]
section of the tyre property file. For details on various implementations see chapter 5.
7/35
Basically USE_MODE (or ISWTCH) = ABCD (e.g. 1134); the following choices can be made:
8/35
NOTE: In principle all combinations are possible, although some make more sense than others.
Typically you do not use road contact for 2D or 3D roads without activating rigid ring dynamics. On the
other hand you may want to use rigid ring dynamics on a flat road surface e.g. in case of ABS/ESP or
shimmy analysis. Obviously the choice of the operating mode will affect the calculation times.
MF-Tyre and MF-Swift
The next table lists the operating modes that are supported by MF-Tyre and MF-Swift licenses.
2.2
MF-Tyre 6.1
0,1,2,3,4
0,1,2
0,1,2,3
0,1,2,3
MF-Swift 6.1
0,1,2,3,4,5
0,1,2,3,4
0,1,2,3,4,5
0,1,2,3
Axis systems
MF-Tyre/MF-Swift 6.1 uses the ISO sign conventions as shown in the figure below.
9/35
Vsx
Vx
tan ( ) =
(note:
Vsy
Vx
In these equations Vx is the x-component (in the wheel centre plane) of the wheel contact centre
horizontal (i.e. parallel to road) velocity V; Vs is the wheel slip velocity, with components Vsx and Vsy,
which is defined as the horizontal velocity of the slip point that is thought to be attached to the wheel
at a distance that equals the effective rolling radius below the wheel centre in the wheel centre plane.
Units
The output of the tyre model is always in SI units (m, N, rad, kg, s).
The tyre property file uses SI units by default (m, N, rad, kg, s); this is always the case when it is
generated by MF-Tool. It is allowed to use a different set of units (e.g. mm or inch for length). The
specification in the [UNITS] section file applies to all parameters in the tyre property file.
The tyre model expects SI units to be passed via the interface between tyre model and the multibody
simulation program, as defined in the specification of the Standard Tyre Interface (STI) [8]. However
many multibody codes do not use units internally and leave the choice of a consistent set of units to
the user. In many cases this implies that the vehicle model has to be defined using SI units to avoid
unit conversion problems.
Please contact TNO if you have special, non-standard requirements with respect to units.
10/35
2.3
Various signals are available for post-processing. Depending on the implementation they are selected
by means of a keyword, signal number or other methods.
[N]
[N]
[N]
[Nm]
[Nm]
[Nm]
slip quantities:
7 kappa
8 alpha
9 gamma
10 phi
[-]
[rad]
[rad]
[1/m]
[m/s]
[m]
[m]
[m]
[m]
[-]
[-]
[m]
(not always available)
[m]
(not always available)
[m/s]
[m/s]
[m/s]
[rad/s]
[m]
(not always available)
[m]
[m]
[m]
[-]
[-]
[-]
[m]
[rad]
[1/m]
Note that the wheel spindle forces and moments are in general obtained from the multibody package.
11/35
3.1
Overview
The tyre property file (*.tir) contains the parameters of the tyre model. Sample tyre property files are
provided with the installation. The file is subdivided in various sections indicated with square brackets.
Each section describes a certain aspect of the tyre behaviour. The next table gives an overview:
Input limitations
[INFLATION_PRESSURE_RANGE]
[VERTICAL_FORCE_RANGE]
[LONG_SLIP_RANGE]
[SLIP_ANGLE_RANGE]
[INCLINATION_ANGLE_RANGE]
minimum
minimum
minimum
minimum
minimum
Magic Formula:
[SCALING_COEFFICIENTS]
[LONGITUDINAL_COEFFICIENTS]
[OVERTURNING_COEFFICIENTS]
[LATERAL_COEFFICIENTS]
[ROLLING_COEFFICIENTS]
[ALIGNING_COEFFICIENTS]
[TURNSLIP_COEFFICIENTS]
and maximum
and maximum
and maximum
and maximum
and maximum
Though at first sight the number of coefficients may seem extensive, Delft-Tyre has established two
methods to significantly facilitate tyre model parameterisation:
1. MF-Tool: this is an automated fitting tool to determine the tyre model parameters and manipulate
the resulting characteristics [8]. Fitting Magic Formula coefficients is a well established process
within the vehicle industry.
Furthermore, MF-Tool features a generic method for identifying MF-Swift parameters from
standardised measurements such as loaded radius, contact length and cleat/drum tests.
2. Reduced input data requirements: if no (or limited) measurement data is available it is also
allowed to omit coefficients in the tyre property file. Built-in procedures will be used to provide a
reasonable estimate for the missing data and only a small number of coefficients are needed. The
next table gives the minimum required coefficients.
When using this reduced parameter file, detailed effects such as combined slip, tyre relaxation
effects and enveloping behaviour on short wavelength road obstacles are included, although the
related parameters are not explicitly specified.
12/35
coefficient
meaning
FITTYP
UNLOADED_RADIUS
MASS
Tyre mass
GRAVITY
FNOMIN
VERTICAL_STIFFNESS
VERTICAL_DAMPING
LONGITUDINAL_STIFFNESS
LATERAL_STIFFNESS
PDX1
PKX1
PDY1
PKY1
PKY2
Tip:
The use of estimated combined slip possibly improves the performance of the tyre model when
extrapolating to (very) low friction values. Estimated combined slip can be turned on by setting
the combined slip coefficients in the tyre property file to zero or by omitting them.
13/35
3.2
Backward compatibility
MF-Tyre/MF-Swift 6.1 is backward compatible with MF-Tyre 5.x, MF-MC-Tyre 1.x , SWIFT 1.x and
MF-Tyre/MF-Swift 6.0.x. Tyre property files generated for these tyre models will work with MFTyre/MF-Swift 6.1 and give the same simulation results as before.
1996
motorcycle tyres
handling
MF-Tyre 5.0
MF-MCTyre 1.0
MF-Tyre 5.1
SWIFT 1.0
2001
MF-Tyre 5.2
MF-MCTyre 1.1
SWIFT 1.1
SWIFT 1.2
2004
MF-Tyre/MF-Swift 6.0
2008
MF-Tyre/MF-Swift 6.1
However some differences may occur at very low speeds when relaxation behaviour is included
combined with a forward velocity below the value specified with the parameter VXLOW in the
[MODEL] section. Due to new formulations the tyre behaviour is much more realistic for these
operating conditions.
In the case of MF-Swift minor differences may occur between the 1.x, 6.0.x and 6.1 versions due to a
different formulation of the contact patch dynamic behaviour. These differences can be observed in
the tyre contact forces and slip values, whereas at wheel axle level the differences remain small.
Due to the built-in estimation procedure it is possible to use for example an existing MF-Tyre 5.2 tyre
property file and perform simulations including turn slip, rigid ring dynamics and tyre enveloping
behaviour, thus already benefiting from the new functionality available in MF-Tyre/MF-Swift 6.1.
14/35
Note 1: the selection of the appropriate set of Magic Formula equations is based on the parameter
FITTYP in the [MODEL] section of the tyre property file. The following conventions apply:
FITTYP=5
FITTYP=6
FITTYP=21
FITTYP=51
FITTYP=52
FITTYP=60
FITTYP=61
MF-Tyre/MF-Swift 6.1 accepts all these values for the parameter FITTYP. It is recommended not to
change the value of the parameter FITTYP unless you are sure that the model parameters in the tyre
property file are meant for that specific Magic Formula version!
Note 2: As described in section 2.1, the modular approach of the tyre model allows a user to select
various combinations of Magic Formula equations, contact methods and dynamics.
Former MF-MCTyre users explicitly will have to select smooth road contact with circular cross
section (B=2) to get the same results using MF-Tyre 6.1 with their MF-MCTyre datasets.
Former SWIFT-Tyre 1.x users will have to select 2D road contact using basic functions (B=4) and
rigid ring dynamics(C=3) to get the same results as before.
Note 3: The camber angle scaling factors LGAX, LGAY and LGAZ are not supported anymore. The
camber influence in MF-Tyre/MF-Swift 6.x can now be more conveniently controlled by the new
parameters LKYC (Fy) and LKZC (Mz). These parameters allow explicit scaling of the camber
stiffness and camber moment stiffness. These new parameters also have to be used in combination
with MF-Tyre 5.x and MF-MCTyre 1.x datasets.
15/35
3.3
Scaling factors
Tyre force and moment testing is often done in a laboratory environment (e.g. using a flat track tyre
tester or a drum). The artificial road surface on the tyre test machine may be quite different from a real
road surface. Combined with other factors like temperature, humidity, wear, inflation pressure, drum
curvature, etc. the tyre behaviour under a vehicle may deviate significantly from the results obtained
from a test machine. Differences of up to 20 % in the friction coefficient and cornering stiffness have
been reported in literature for a tyre tested on different road surfaces compared to lab measurements.
For this purpose scaling factors are included in the tyre model, which allow the user to manipulate and
tune the tyre characteristics, for example to get a better match between full vehicle tests and
simulation model. Another application of the scaling factors is that they may be used to eliminate
some undesired offsets or shifts in the Magic Formula.
The most important scaling factors are:
LMUX longitudinal peak friction coefficient (Fx)
LKX
longitudinal slip stiffness (Fx)
LMUY lateral peak friction coefficient (Fy)
LKY
cornering stiffness (Fy)
LKYC camber stiffness (Fy)
LTR
pneumatic trail (Mz)
LKZC camber moment stiffness (Mz)
LMP
parking moment at standstill (Mz)
Normally when processing the tyre measurements these scaling factors are set to 1, but when doing a
validation study on a full vehicle model they can be adjusted to tune the tyre behaviour. The scaling
factors are defined in the [SCALING_COEFFICIENTS] section of the tyre property file.
16/35
3.4
The following table lists the required and optional parameters for each tyre model version. For
convenience, a comparison is made with the previous model versions.
x: required parameter
MF-MCTyre 1.1
SWIFT 1.2
MF-Tyre 5.2
MF-Swift 6.0
MF-Tyre 6.0
MF-Swift 6.1
MF-Tyre 6.1
[MODEL]
FITTYP
TYRESIDE
61 61 60 60
x
LONGVL
Reference speed
VXLOW
ROAD_INCREMENT
ROAD_DIRECTION
PROPERTY_FILE_FORMAT
USE_MODE
HMAX_LOCAL
TIME_SWITCH_INTEG
21 52
[DIMENSION]
UNLOADED_RADIUS
WIDTH
RIM_RADIUS
RIM_WIDTH
Rim width
ASPECT_RATIO
INFLPRES
NOMPRES
MASS
Tyre mass
IXX
IYY
BELT_MASS
Belt mass
BELT_IXX
BELT_IYY
GRAVITY
M_B
I_BY
I_BXZ
C_GRV
Gravity constant
[OPERATING_CONDITIONS]
[INERTIA]
x
[VERTICAL]
FNOMIN
17/35
MF-Tyre 6.0
MF-Swift 6.0
MF-Tyre 5.2
MF-MCTyre 1.1
MF-Swift 6.1
VERTICAL_DAMPING
MC_CONTOUR_A
MC_CONTOUR_B
BREFF
DREFF
FREFF
Q_RE0
Q_V1
Q_V2
Q_FZ2
Q_FCX
Q_FCY
Q_CAM
PFZ1
BOTTOM_OFFST
BOTTOM_STIFF
LATERAL_STIFFNESS
YAW_STIFFNESS
FREQ_LONG
FREQ_LAT
FREQ_YAW
FREQ_WINDUP
DAMP_LONG
DAMP_LAT
DAMP_YAW
SWIFT 1.2
MF-Tyre 6.1
VERTICAL_STIFFNESS
[STRUCTURAL]
DAMP_WINDUP
DAMP_RESIDUAL
x
x
x
x
DAMP_VLOW
Q_BVX
Q_BVT
PCFX1
PCFX2
PCFX3
PCFY1
PCFY2
PCFY3
PCMZ1
C_BX0
C_RX
C_BT0
C_BY
C_RY
C_BGAM
C_RP
K_BX
18/35
K_BT
K_BY
K_BGAM
MF-MCTyre 1.1
SWIFT 1.2
MF-Tyre 5.2
MF-Swift 6.0
MF-Tyre 6.0
MF-Swift 6.1
MF-Tyre 6.1
[CONTACT_PATCH]
Q_RA1
Q_RA2
Q_RB1
Q_RB2
ELLIPS_SHIFT
ELLIPS_LENGTH
ELLIPS_HEIGHT
ELLIPS_ORDER
Order of ellipsoid
ELLIPS_MAX_STEP
ELLIPS_NWIDTH
ELLIPS_NLENGTH
Q_A2
Q_A1
ELLIPS_INC
Q_LBF
Q_LOS1
Q_LOS2
Q_LIMP1
Q_LIMP3
Q_LIMP4
Q_LIMP2
[INFLATION_PRESSURE_RANGE]
PRESMIN
PRESMAX
FZMIN
FZMAX
KPUMIN
KPUMAX
ALPMIN
ALPMAX
[VERTICAL_FORCE_RANGE]
[LONG_SLIP_RANGE]
[SLIP_ANGLE_RANGE]
[INCLINATION_ANGLE_RANGE]
CAMMIN
CAMMAX
[SCALING_COEFFICIENTS]
LFZO
19/35
MF-Tyre 6.0
MF-Swift 6.0
MF-Tyre 5.2
MF-MCTyre 1.1
MF-Swift 6.1
LMUX
LEX
LKX
LHX
LVX
LCY
LMUY
LEY
LKY
LKYC
LKZC
LHY
LVY
LTR
LRES
LXAL
LYKA
LVYKA
LS
LMX
LVMX
LMY
LMP
SWIFT 1.2
MF-Tyre 6.1
LCX
LKC
LCC
LEC
LSGKP
LSGAL
LGYR
[LONGITUDINAL_COEFFICIENTS]
PCX1
PDX1
PDX2
PDX3
PEX1
PEX2
PEX3
PEX4
PKX1
PKX2
PKX3
PHX1
PHX2
PVX1
PVX2
RBX1
RBX2
20/35
MF-Swift 6.0
MF-MCTyre 1.1
MF-Tyre 6.0
REX1
REX2
RHX1
PPX1
PPX2
PPX3
PPX4
PTX1
PTX2
PTX3
SWIFT 1.2
MF-Swift 6.1
RCX1
MF-Tyre 5.2
MF-Tyre 6.1
RBX3
[OVERTURNING_COEFFICIENTS]
QSX1
QSX2
QSX3
QSX4
QSX5
QSX6
QSX7
QSX8
QSX9
QSX10
QSX11
QSX12
QSX13
QSX14
PPMX1
PCY1
PDY1
PDY2
PDY3
PEY1
PEY2
PEY3
PEY4
PEY5
PKY1
PKY2
PKY3
PKY4
PKY5
PKY6
PKY7
PHY1
PHY2
[LATERAL_COEFFICIENTS]
21/35
MF-Tyre 6.0
MF-Swift 6.0
MF-Tyre 5.2
MF-MCTyre 1.1
MF-Swift 6.1
PVY2
PVY3
PVY4
RBY1
RBY2
RBY3
RBY4
RCY1
REY1
REY2
RHY1
RHY2
RVY1
RVY2
RVY3
RVY4
RVY5
RVY6
PPY1
PPY2
PPY3
PPY4
PPY5
PCY2
PHY3
PTY1
PTY2
PTY3
SWIFT 1.2
MF-Tyre 6.1
PVY1
[ROLLING_COEFFICIENTS]
QSY1
QSY2
QSY3
QSY4
QSY5
QSY6
QSY7
QSY8
QBZ1
QBZ2
QBZ3
QBZ4
QBZ5
QBZ9
QBZ10
QCZ1
[ALIGNING_COEFFICIENTS]
22/35
MF-Tyre 6.0
MF-Swift 6.0
MF-Tyre 5.2
MF-MCTyre 1.1
MF-Swift 6.1
QDZ2
QDZ3
QDZ4
QDZ6
QDZ7
QDZ8
QDZ9
QDZ10
QDZ11
QEZ1
QEZ2
QEZ3
QEZ4
QEZ5
QHZ1
QHZ2
QHZ3
QHZ4
SSZ1
SSZ2
SSZ3
SSZ4
PPZ1
PPZ2
QTZ1
MBELT
SWIFT 1.2
MF-Tyre 6.1
QDZ1
x
x
[TURNSLIP_COEFFICIENTS]
PDXP1
PDXP2
PDXP3
PKYP1
PDYP1
PDYP2
PDYP3
PDYP4
PHYP1
PHYP2
PHYP3
PHYP4
PECP1
PECP2
Camber w.r.t. spin reduction factor varying with load parameter in camber stiffness
QDTP1
QCRP1
QCRP2
Turn slip moment (at alpha=90deg) parameter for increase with spin
QBRP1
QDRP1
23/35
MF-Tyre 5.2
SWIFT 1.2
MF-MCTyre 1.1
Obsolete parameters which may be in a tyre property file, but are ignored by MF-Tyre/MF-Swift 6.x
TYPE
MFSAFE1
MFSAFE2
MFSAFE3
M_A
I_AY
I_AXZ
M_R
I_R
K_RX
K_RY
K_RP
FLT_A
Q_KC1
Q_KC2
description
[MODEL]
[SHAPE]
[INERTIA]
[STRUCTURAL]
[VERTICAL]
BOTTOM_TRNSF
[CONTACT_PATCH]
[SCALING_COEFFICIENTS]
LGAX
LGAY
10
LGAZ
11
in MF-Swift 6.0 and 6.1 a new formulation is used without residual mass
parameter deleted
10
11
24/35
! Comments section
$--------------------------------------------------------------------------UNITS
[UNITS]
LENGTH
= 'meter'
FORCE
= 'newton'
ANGLE
= 'degree'
MASS
= 'kg'
TIME
= 'sec'
$--------------------------------------------------------------------------MODEL
[MODEL]
ROAD_TYPE
= '...'
$---------------------------------------------------------------------PARAMETERS
[PARAMETERS]
...
In the [UNITS] section, the units that are used in the road data file are set. The [MODEL] section is
used to specify the road type, see listing above. The [PARAMETERS] section contains general
parameters and road surface type specific parameters. The general parameters are listed below:
25/35
General
MU
OFFSET
ROTATION_ANGLE_XY_PLANE
DRUM_RADIUS
The road surface type specific parameters are explained in the next sections:
Plank Road
HEIGHT
START
Distance along the X-axis of the road to the start of the cleat.
LENGTH
BEVEL_EDGE_LENGTH
DIRECTION
Rotation of the cleat about the Z-axis with respect to the Yaxis of the road. If the cleat is placed crosswise,
DIRECTION = 0. If the cleat is along the X-axis,
DIRECTION = 90.
DIRECTION
z
+
x
START
LENGTH
HEIGHT
Vx
Polyline
The [PARAMETERS] block must have a (XZ_DATA) subblock. The subblock consists of three
columns of numerical data:
Column one is a set of X-values in ascending order;
Columns two and three are sets of respective Z-values for left and right track.
Example:
26/35
[PARAMETERS]
MU
= 1.0
OFFSET
ROTATION_ANGLE_XY_PLANE
$
$ X_road
Z_left Z_right
(XZ_DATA)
-1.0e04
0.0500
...
...
...
Sine Road
HEIGHT
START
Distance along the X-axis of the road to the start of the sine
wave.
LENGTH
DIRECTION
Rotation of the bump about the Z-axis with respect to the Xaxis of the road. If the bump is placed crosswise,
DIRECTION = 0. If the bump is along the X-axis,
DIRECTION = 90.
N_BUMPS
Finally, sample road data files are provided with the installation.
27/35
5.1
ADAMS
MF-Tyre/MF-Swift 6.1 is offered as a user programmed tyre in ADAMS. To use the TNO tyre model
you need a customised ADAMS solver. These are included in the delivery. The next table gives an
overview of supported ADAMS versions and operating systems.
ADAMS
operating system
version
Windows
Linux
HP-UX
2003
2005
2005r2
2007r1
under development
='USER'
USER_SUB_ID
= 815
This ensures that the TNO MF-Tyre/MF-Swift 6.1 tyre model is called. This can also be checked in the
ADAMS message file (*.msg), the following statement should appear:
TYR815 -> DELFT-TYRE MF-Tyre/MF-Swift 6.1 xxxxxxxx-x
28/35
global integration: the tyre differential equations are solved in the ADAMS solver together
with the multibody equations
local integration: the tyre differential equations are solved locally inside the tyre model
independent of the multi-body model
Local integration can significantly speed up the simulation time when using rigid ring dynamics on an
uneven road surface. For calculations on a level road surface without rigid ring dynamics a global
integration will be faster and more accurate. The parameters for this local integrator inside the tyre
model are set in [MODEL] section of the tyre property file, for example:
HMAX_LOCAL
= 0.00025
TIME_SWITCH_INTEG
= 0.1
HMAX_LOCAL defines the step size of the local integrator, too big values may result in instability and
generally 0.25 ms is a safe value. TIME_SWITCH_INTEG defines the time when the switch is made
from global to local integration. It is possible to have ADAMS calculate static equilibrium for the tyre
model and at a later stage during the simulation switch to local integration to speed it up.
Switching between local and global integration is only possible if a sufficient states are available in the
ADAMS model. The ADAMS message file will provide additional information on this. Some examples:
LOCAL
LOCAL
29/35
Global integration is only possible when the first number is smaller than or equal to the second one.
The number of states available is defined by the tyre GSE.
NOTE 1: when using local integration the maximum step size HMAX of the ADAMS integrator has to
be set to 1 ms or smaller, otherwise the simulation results may become inaccurate or unstable.
NOTE 2: to use global integration (if possible), comment out the line defining HMAX_LOCAL from the
tyre property file by using a $ or ! character.
30/35
5.2
MATLAB/Simulink/SimMechanics
MF-Tyre/MF-Swift 6.1 is offered for MATLAB/Simulink 6.5 and up. The command dteval can be used
to evaluate the Magic Formula model for series of input variables. For more information on dteval,
please type help dteval on the MATLAB command line. For simulation model development in
MATLAB 2006a and up, blocks are available from the library TNO_dtlib.mdl in Simulink.
In addition to the normal functionality, the Simulink and SimMechanics blocks allow a user to change
tyre scaling factors as a function of time or any other signal available in the model. Further, some
blocks are provided to easily model moving and non-moving road surfaces, coordinate system
transformation and animation of the wheel using the Virtual Reality Toolbox. See the help function of
the blocks and the Simulink and SimMechanics demos for more information.
backward compatibility
For older versions of MATLAB (6.5 and up) the library TNO_dtlib_v65.mdl in Simulink can be used.
The only difference with respect to the latest library TNO_dtlib.mdl is that SimMechanics is not
supported.
The MATLAB command line functions mfread and mfeval have been replaced by the new function
dteval.
The sequence of the signals in the output vector (varinf) in the Simulink tyre block has changed.
Please use the help function of this block to learn more about the new definition. In addition a Bus
Selector block may be used to select the appropriate output signals based on their names.
31/35
= 1.391 kgm2
wheel Iyy
= 2.736 kgm2
tyre mass
tyre Ixx
= 0.326 kgm2)
tyre Iyy
= 0.636 kgm2)
Note: When switching on rigid ring dynamics the mass/inertia distribution is adjusted in such a way
that the mass and inertia properties of the complete wheel (rim+tyre) remain unchanged.
Initialisation
When using rigid ring + initial statics the tyre model will give the following messages:
Delft-Tyre
4721.4 N
0.3038 m
angular velocity
You can use this information to set the correct angular velocity of the wheel when specifying the initial
conditions in your model.
32/35
5.3
LMS DADS
MF-Tyre/MF-Swift 6.1 is offered for DADS 9.6. To introduce the tyre model and to change the tyre
model settings (tyre property file, scale factors, etc.) in the DADS GUI, select Force, Tire, STI in the
DADS modelling panel:
To plot the tyre model outputs after having performed a simulation, open the DADSGraph menu and
select tire element and the signal you want to plot:
33/35
5.4
MF-Tyre/MF-Swift 6.x is also available in third party simulation software. Some examples are:
Virtual.Lab (LMS), SIMPACK (INTEC), MADYMO (TASS), CarSim/BikeSim/TruckSim (MSC). Please
contact your simulation package supplier or TNO for more information.
34/35
6 References
[1]
[2]
Pacejka, H.B., I.J.M. Besselink: Magic Formula Tyre model with Transient Properties,
Supplement to Vehicle System Dynamics, Vol. 27, pp. 234-249, 1997.
[3]
Zegelaar, P.W.A., The Dynamic Response of Tyres to Brake Torque Variations and Road
Unevenesses, dissertation, Delft University of Technology, The Netherlands, 1998.
[4]
Maurice, J.P., Short Wavelength and Dynamic Tyre Behaviour under Lateral and Combined
Slip Conditions, dissertation, Delft University of Technology, The Netherlands, 1999.
[5]
[6]
Besselink, I.J.M., H.B. Pacejka, A.J.C. Schmeitz, S.T.H. Jansen: The SWIFT tyre model:
overview and applications, Presented at the AVEC 2004: 7th International Symposium on
Advanced Vehicle Control, 23-27 August 2004.
[7]
A. Riedel, J.J.M. van Oosten: Standard Tyre Interface, Release 1.4. Presented at 2nd
International Colloquium on Tyre Models for Vehicle Dynamics Analysis, February 20-21
1997. Issued by the TYDEX - Working group.
[8]
TNO Automotive: MF-Tool 6.1 Users Manual, TNO Automotive, The Netherlands, 2008.
35/35