Pacejka 2012
Tire component with Pacejka 2012 formulation and visualization
Description
Details
Equations
Connections
Parameters
References
The Pacejka 2012 tire component employs the 2012 formulation of the Pacejka tire model, presented in [1].
The tire geometry is assumed to be a thin circular disk, which is common in automotive applications. A single point contact is considered for the tire-ground interaction.
The tire kinematics used in this component are described in detail in Tire Kinematics.
Several options are available for defining the surface on which the tire is operating. These options are explained in Surface.
Tire Parameters Block
The Pacejka 2012 tire model has about 180 parameters. Unlike the Linear and the Fiala tire components, where the required parameters are defined in the GUI, to facilitate the parameter handling process the Pacejka Parameters App should be used to generate a parameter block which contains the necessary tire parameters. To open this app, browse to Add Apps or Templates > Tires > Pacejka Parameters. The generated parameter block will be located in the Local Components panel on the left side of the GUI.
The user should place the generated parameter block into the workspace at the same or higher level as the Pacejka tire components that it defines.
Override Parameters
There is an Override check box in the Inertia, Radial Compliance, and Scaling Factors sections of the Pacejka tire component properties.
Enabling one of these check boxes allows the user to override the associated parameters otherwise defined in the tire parameters block. For example, the user can override the inertia properties as shown below.
Checking an Override check box also exposes the associated parameters to apps such as the Parameter Sweep app and the FMU Generation app.
ISO Axis
Unlike the Linear and Fiala tire models, the Pacejka tire model is typically asymmetric, that is Fx⁡−κ≠−Fx⁡κ or Fy⁡−α≠−Fy⁡α. To ensure the correct formulation, the ISO X axis of a tire should point towards the heading of the vehicle. The Show ISO axis option in the visualization section of the tire parameters can be helpful to visually confirm that the ISO axes have been assigned correctly.
If not assigned correctly, the user can change the integer parameter of ISO from 0 to 1 to rotate the ISO axis 180 degrees around ISO Z.
Sideness
The Pacejka tire parameters apply to a specific tire side. This denotes the side of the vehicle where the tire should be mounted. The Side parameter in the properties can be used to mirror the tire. For example, if the parameters of the generated parameter block are for a right side tire, then the tire components mounted on the right side of the vehicle model should be used with Side=0, and those on the left side should have Side=1.
Normal Force and Effective Radius
The normal force exerted by the surface to the tire is calculated using the given compliance parameters and surface geometry. There are two implemented formulations in the Pacejka tire component for calculating the normal force: Pacejka formulation and Linear spring-damper.
Pacejka Formulation
The Pacejka formulation option uses the following equation for the normal force [1]
Fz=1+qV2⁢Ω⁢R0V0−qFcx⁢FxFz02−qFcy⁢FyFz02⁢γ2⁢qFz3+qFz1⁢ρzR0+q_Fz2 ⁢ρz2R02⁢dpi⁢ppFz1+1⁢Fz0
Note that with this option selected for the normal force, the Pacejka effective radius formulation will also be used internally. This formulation is as follows [1]
reff=r_Ω−Fz0⁢Freff⁢FzFz0+Dreff⁢arctan⁡Breff⁢FzFz0C
rΩ=R0⁢qreo+qV1⁢R0⁢ΩV02
where the nominal load, Fz0 and the rest of the parameters used in these equations are defined in the tire parameters block.
Linear Spring-Damper
The Linear spring-damper option is the same formulation used for normal force calculation in the Linear and the Fiala tire components as explained below.
The tire loaded radius is calculated using the distance of the tire center from the surface, rz (see Surface), and the inclination angle, γ (see Tire Kinematics).
rL=rzcos⁡γ
Using a linear spring and saturated damping forces based on the tire compliance, the normal force, Fz, is calculated as follows
FzC={C⁢R0−rLrL<R00otherwise
FzK={K⁢VzrL<R00otherwise
Fz={FzC+min⁡FzC,FzK0<FzC+FzK0otherwise
where Vz is the tire center speed with respect to ISO Z, C is tire stiffness, K is tire damping, and R0 is tire unloaded radius. The use of the min function is to ensure that Fz is continuous at rL=R0.
With this option selected for Fz, the user can choose between the tire unloaded radius, R0, and the loaded radius, rL, to assign to the effective radius, reff.
Slip Calculations
Four options are available for tire slip calculation, Quasi-static, Constant time lags, Stretched string, and Damped transient.
Quasi-static
With the choice of Quasi-static, the following equations for longitudinal slip, κ, and slip angle, α, hold true on a flat surface with no inclination angle
κ={Ω⁢re−VxVx|Vx|>Vx⁡minΩ⁢re−Vx⁢2⁢Vx⁡minVx2+Vx⁡min2otherwise
tan⁡α={VyVx|Vx|>Vx⁡minVy⁢2⁢Vx⁡minVx2+Vx⁡min2otherwise
where re is the tire effective radius and considered equal to the loaded radius (rL), Ω is the tire speed of revolution, Vx and Vy are the speeds of the tire center with respect to ISO X and ISO Y axes, respectively, and Vx⁡min is the velocity threshold used for singularity avoidance in the slip calculations.
The longitudinal slip and slip angle are continuous and differentiable in the neighborhood of Vx=0.
Constant Time Lags
A first-order dynamics to the longitudinal slip and slip angle calculation can be introduced using the Constant time lags option. When active, the following slip formulation is used:
Tlong⁢dκdt=re⁢Ω−Vx−κ⁢Vx
Tlat⁢dtan⁡αdt=Vy−tan⁡α⁢Vx
Stretched String
With this option active, the relaxation lengths will be used in slip calculation as follows
σlong⁢dκdt=re⁢Ω−Vx−κ⁢Vx
σlat⁢dtan⁡αdt=Vy−tan⁡α⁢Vx
where
σlong=max⁡Fz⁢pTx1+pTx2⁢dfz⁢exp⁡pTx3⁢dfz⁢R0F′z0⁢LSκ,σlongmin
σlat=max⁡pTy1⁢sin⁡2⁢arctan⁡FzpTy2⁢F′z0⁢1−pTy3⁢|γ|⁢R0⁢LSα,σlatmin
Parameters in the above equations should be inserted using the GUI.
The load ratio, dfz, is defined as
dfz=Fz−Fz0Fz0
Damped Transient
The implementation of this formulation is according to the section 7.2.2 of [1]. It is a semi-non-linear model that covers the non-linear range of the slip characteristics and also prevents undamped oscillations at low longitudinal speeds.
For the longitudinal deflection
dudt={0Vx<Vlow∧Aκ⁢κsl<κ′∧u⋅Vsx+Vx⁢uσκ<0−Vsx−Vx⁢uσκotherwise
where Vsx is the x-component of the slip velocity, Vlow is the velocity threshold, Aκ is the parameter for the inequality, and σκ is the longitudinal slip relaxation length. Moreover,
κ′=u⁢κslσκ
κsl=3⁢DxCF,κ
CF,κ=Bx⁢Cx⁢Dx
The following equation computes the longitudinal slip
κ=κ′−kVlow⁢VsxCF,κ
The low speed coefficient is calculated using the following equation
kVlow={12⁢kVlow,0⋅1+cos⁡π⁢VxVlowVx<Vlow0otherwise
Similarly, for the lateral deflection
dvdt={0Vx<Vlow∧αsl<tanα′∧v⋅−Vsy+Vx⁢vσα<0Vsy−Vx⁢vσαotherwise
tanα′=vσα
αsl=3⁢DyCF,α
CF,α=By⁢Cy⁢Dy
tanα=tanα′−Vsy⁢kVlowCF,α
The slip angle itself can be calculated from tanα. This might not be necessary because the Pacejka 2012 formulation is based on tanα.
The formulation for resultant forces/moments of tire-surface interaction at the tire contact patch are summarized below for the Pacejka 2012 tire component.
The longitudinal force is
Fx=Gxα⁢Fx0
Fx0=Dx⁢sin⁡Cx⁢arctan⁡Bx⁢κx−Ex⁢Bx⁢κx−arctan⁡Bx⁢κx+SVx
Gxα=1Gxα0⁢cos⁡Cxα⁢arctan⁢Bxα⁢αS−Exα⁢Bxα⁢αS−arctan⁡Bxα⁢αS
The lateral force is
Fy=Fy0⁢Gyκ+SVκ
Fy0=Dy⁢sin⁡Cy⁢arctan⁢By⁢αy−Ey⁢By⁢αy−arctan⁡By⁢αy+SVy
Gyκ=1Gyκ0⁢cos⁡Cyκ⁢arctan⁢Byκ⁢κS−Eyκ⁡Byκ⁢κS−arctan⁡Byκ⁢κS
The normal force, Fz, has been discussed in the Normal Force and Effective Radius section.
The overturning couple is
Mx=R0⁢Fz⁢(qsx1⁢λVMx−qsx2⁢γ⁡dpi⁢ppM⁢x1+1+qsx3⁢FyFz0+qsx4⁢cos⁢qsx5⁢arctan⁡qsx6⁢FzFz02⁢sin⁢qsx7⁢γ+qsx8⁢arctan⁡qsx9⁢FyFz0+qsx10⁢arctan⁡qsx11⁢FzFz0⁢γ)⁢λMx
The rolling resistance moment is
My=FzR0⁢qsy1+qsy2FxFz0+qsy3|VxV0|+qsy4VxV04+qsy5+qsy6FzFz0γ2⁢FzFz0qsy7⁢pipioqsy8⁢λMy
The self-aligning torque is
Mz=M′z+Mzr+s⁢Fx
where M′z is the torque due to pneumatic trail, t, Mzr is the residual torque, and s⁢Fx is the longitudinal force contribution to the self-aligning torque. Each of these terms has a specific expression, discussed in [1] in more detail.
All the employed parameters in the equations above need to be defined in the tire parameters block and be accessible to the Pacejka tire components.
Name
Modelica ID
framea
Multibody frame for tire center
frame_a
Fz
Signal output for the normal force
LongSlip
Signal output for longitudinal slip
SlipAng
Signal output for slip angle
SpinRate
Signal output for tire speed of revolution or spin rate
reff
Signal output for tire effective radius
r_eff
IncAng
Signal output for tire inclination angle or camber
enin
[1] Vector input for surface normal vector
en_in
rc
[1] Vector output for tire center position w.r.t. the inertial frame
r_c
rzin
[1] Signal input for tire center distance from the surface
rz_in
pin
[2] Signal input for tire inflations pressure
p_in
[1] Available if Surface parameters Flat surface is false and Defined externally is true.
[2] Available if Pressure Override parameter is true and Constant Pressure is false.
Inertia
Default
Units
Override
false
True (checked) overrides the inertia parameters and enable the following two parameters
useInertia
m
28
kg
Tire mass
Mass
[I]
[1]
kg⁢m2
Rotational inertia, expressed in frame_a (center of tire)
[1] 0.780001.560000.78
Initial Conditions
Use Initial Conditions
True (checked) enables the following parameters
useICs
ICr,v
Ignore
Indicates whether to ignore, try to enforce, or strictly enforce the translational initial conditions
MechTranTree
r&conjugate0;0
0,0,0
Initial displacement of frame_a (tire center) at the start of the simulation expressed in the inertial frame
InitPos
Velocity Frame
Inertial
Indicates whether the initial velocity is expressed in frame_a or inertial frame
VelType
v&conjugate0;0
ms
Initial velocity of frame_a (tire center) at the start of the simulation expressed in the frame selected in Velocity Frame
InitVel
ICθ,ω
Indicates whether to ignore, try to enforce, or strictly enforce the rotational initial conditions
MechRotTree
Quaternions
Indicates whether the 3D rotations will be represented as a 4 parameter quaternion or 3 Euler angles. Regardless of setting, the initial orientation is specified with Euler angles.
useQuats
Euler Sequence
1,2,3
Indicates the sequence of body-fixed rotations used to describe the initial orientation of frame_a (center of mass). For example, [1, 2, 3] refers to sequential rotations about the x, then y, then z axis (123 - Euler angles)
RotType
θ&conjugate0;0
rad
Initial rotation of frame_a (center of tire) at the start of the simulation (based on Euler Sequence selection)
InitAng
Angular Velocity Frame
Euler
Indicates whether the initial angular velocity is expressed in frame_a (body) or the inertial frame. If Euler is chosen, the initial angular velocities are assumed to be the direct derivatives of the Euler angles.
AngVelType
ω&conjugate0;0
rads
Initial angular velocity of frame_a (center of tire) at the start of the simulation expressed in the frame selected in Angular Velocity Frame
InitAngVel
Settings
e^spin
Tire's spin axis (local)
SymAxis
Side
0
0: default, 1: mirrored
ISO
0: Keep ISO, 1: Rotate ISO pi radians around Z axis
intISO
Pacejka formulation
Normal force equation (Pacejka formulation or Linear spring-damper)
FzMode
Loaded radius
Effective radius (Loaded radius or Unloaded radius)
reffMode
Slip
Choose type of slip calculation (Quasi-static, Constant time lags, Stretched string, or Damped transient)
slipMode
Tlong
0.3
s
Time lag for longitudinal slip
TlongIn
Tlat
Time lag for slip angle
TlatIn
Params
[2]
Parameters for stretched-string formulation: [LSkappa, LSalpha, p_Tx1, p_Tx2, p_Tx3, p_Ty1, p_Ty2, p_Ty3]
ssParams
σlongmin
0.1
Minimum longitudinal relaxation length
TlongMin
σlatmin
Minimum lateral relaxation length
TlatMin
Vlow
10
Velocity threshold to add damping at lower speeds
V_low
kVlow0
10000
N⁢sm
Maximum damping value
k_V_low0
Aκ
1
Inequality constant for the longitudinal slip
A_kappa
σκ
0.05
Longitudinal slip relaxation length
sigma_kappa
σα
Side slip relaxation length
sigma_alpha
[2] 1,1,2.3657,1.4112,0.56626,2.1439,1.9829,−0.90729
Pressure
True (checked) overrides override the pressure parameters and enables the following parameters
overridePressure
pio
2.2·105
Pa
Nominal tire pressure
over_p_io
Constantpressure
true
True (checked) uses constant pressure; false provides an input port for the tire pressure
isConstantPressure
pi
2.4·105
Tire pressure
over_p_cons
Scaling Factors
True (checked) override the scaling factors and enables the following parameter
overrideFactors
λFz0
Nominal load scaling factor
over_lambda_Fz0
λmux
Peak friction coefficient (x) scaling factor
over_lambda_mux
λmuy
Peak friction coefficient (y) scaling factor
over_lambda_muy
λmuV
Slip speed decaying friction scaling factor
over_lambda_muV
λKxκ
Brake slip stiffness scaling factor
over_lambda_KxKap
λKyα
Cornering stiffness scaling factor
over_lambda_KyAlp
λCx
Shape factor (x) scaling factor
over_lambda_Cx
λCy
Shape factor (y) scaling factor
over_lambda_Cy
λEx
Curvature factor (x) scaling factor
over_lambda_Ex
λEy
Curvature factor (y) scaling factor
over_lambda_Ey
λHx
Horizontal shift (x) scaling factor
over_lambda_Hx
λHy
Horizontal shift (y) scaling factor
over_lambda_Hy
λVx
Vertical shift (x) scaling factor
over_lambda_Vx
λVy
Vertical shift (y) scaling factor
over_lambda_Vy
λKyγ
Camber force stiffness scaling factor
over_lambda_KyGam
λKzγ
Camber torque stiffness scaling factor
over_lambda_KzGam
λt
Pneumatic trail scaling factor
over_lambda_t
λMr
Residual torque scaling factor
over_lambda_Mr
λxα
Alpha influence on Fx (kappa) scaling factor
over_lambda_xAlp
λyκ
Kappa influence on Fy (alpha) scaling factor
over_lambda_yKap
λVyκ
Kappa induces ply-steer Fy scaling factor
over_lambda_VyKap
λs
Mz moment arm of Fx scaling factor
over_lambda_s
λCz
Radial tire stiffness scaling factor
over_lambda_Cz
λMx
Overturning couple stiffness scaling factor
over_lambda_Mx
λMy
Rolling resistance moment scaling factor
over_lambda_My
λVMx
Overturning couple vertical shift scaling factor
over_lambda_VMx
Surface
Flat surface
True (checked) means the road surface is assumed flat. It is defined by a plane passing through (0,0,0) and the normal vector given by e^g
flatSurface
Defined externally
True (checked) means the road surface is defined external to the tire component. Additional input and output signal ports are activated.
externallyDefined
δL
0.01
Base distance for local surface patch approximation
deltaL
Data source
inline
Data source for the uneven surface. See following table.
datasourcemode
Surface data
Surface data; matrix or attached data set
table or data
Smoothness
linear
Smoothness of table interpolation
smoothness
nIter
2
Number of iterations to find the contact point candidate, recommended value between 1 and 5
Content of Data source matrix.
Surface normal
First Column
First Row
Global Z
x values
y values
Global Y
z values
Global X
Visualization
Show tire
True (checked) creates a tire visualization and enables following three parameters
showTire
Dw
Tire width (for visualization)
D_w
Tire color
black
color00
Band color
yellow
Tire band color
color01
Tire transparency
True (checked) means the tire is transparent
transparent0
Show force arrow
True (checked) display a force vector and enables the following three parameters
showForceArrow
Show components
True (checked) means three arrows for force components in ISO axes will be shown instead of a single total force arrow
showForceComponents
Force arrow color
red
Specifies the color of the force arrow
color1
Force arrow transparency
True (checked means the force arrow is transparent
transparent1
Force arrow scale
Nm
Scales the length of the force arrow
scale1
Show torque arrow
True (checked) displays a torque vector and enables the following three parameters
showMomentArrow
True (checked) means three arrows for torque components in ISO axes will be shown instead of a single total torque arrow
showMomentComponents
Torque arrow color
blue
Specifies the color of the torque arrow
color2
Torque arrow transparency
True (checked) means the torque arrow is transparent
transparent2
Torque arrow scale
N⁢mm
Scales the length of the torque arrow
scale2
Show tangent plane
True (checked) displays the tangent plane of the contact patch and enables the following four parameters
ShowTanSurface
th0
Patch visualization thickness
rp
0.2
Patch visualization radius
r_patch
Patch color
Green
Color of the contact patch
color3
Patch transparency
True (checked) means contact patch is transparent
transparent3
Show ISO axis
True (checked) displays the ISO axes and enables the following two parameters
showISO
Axis scale
Length of each XYZ ISO axis in the visualization
scaleISO
Axis transparency
True (checked) means axes are transparent
transparentISO
Advanced Parameters
Vx⁡min
Velocity threshold used for singularity avoidance in the slip calculations
V_x_min
εsgn
0.001
Used to smooth sign⁡x as tanh⁡xεsgn
epsilon_sign
εnorm
1·10−8
Used to prevent singularity in vector normalization
epsilon_norm
[1] Pacejka, Hans. Tire and vehicle dynamics. Butterworth-Heinemann, 2012.
See Also
Tire Kinematics
Tires
Download Help Document