Custom
User-defined tire component
Description
Details
Equations
Connections
Parameters
The Custom tire component provides signal inputs for the user to externally feed in the forces and moments with respect to an ISO axis at the contact patch. To calculate these forces/moments, the user can employ the output signals of the Custom tire component such as longitudinal slip or slip angle.
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.
Normal Force
The normal force exerted by the surface to the tire is calculated using the given compliance parameters and surface geometry.
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.
Slip Calculations
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
Above, re is the tire effective radius and is considered to be equal to the loaded radius, rL for the tire component. Ω is the tire speed of revolution, and Vx and Vy are the speeds of the tire center with respect to ISO X and ISO Y axes, respectively. The component code implementation is such that the longitudinal slip and slip angle are continuous and differentiable in the neighborhood of Vx=0.
Using Time Lags
A first-order dynamics to the longitudinal slip and slip angle calculation can be introduced in the Time Lags section of the component properties. When active, the following slip formulation will be used:
Tlong⁢dκdt=re⁢Ω−Vx−κ⁢Vx
Tlat⁢ddttan⁡α=Vy−tan⁡α⁢Vx
The formulation for resultant forces/moments of tire-surface interaction at the tire contact patch are summarized below for the linear tire component (see LinearTire).
Fx=FxISOFy=FyISO
Fz={FzISO¬calculate Fz internallyFzInternalotherwise
Mx=MxISOMy=MyISOMz=MzISO
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
rc
[1] Vector signal output for tire center position w.r.t. the inertial frame
r_c
enin
[1] Vector signal input for surface normal vector
en_in
rzin
[1] Signal input for tire center distance from the surface
rz_in
FxISO
Signal input for the longitudinal force at the contact patch w.r.t. the ISO axis
FyISO
Signal input for the lateral force at the contact patch w.r.t. the ISO axis
FzISO
[2] Signal input for the normal force at the contact patch w.r.t. the ISO axis
MxISO
Signal input for the overturning couple at the contact patch w.r.t. the ISO axis
MyISO
Signal input for the rolling resistance moment at the contact patch w.r.t. the ISO axis
MzISO
Signal input for the self-aligning torque at the contact patch w.r.t. the ISO axis
[1] Available if Surface parameters Flat surface is false and Defined externally is true.
[2] Available if Settings parameter Calculate Fz internally is false.
Inertia
Default
Units
Use inertia
false
True (checked) means use mass and inertia parameters for tire 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
0,1,0
Tire's spin axis (local)
SymAxis
ISO
0
0: keep ISO, 1: rotate π radians around Z-axis
intISO
Calculate Fz internally
True (checked) means Fz is calculated internally using the compliance parameters
useFz
Radial Compliance
These parameters define the radial compliance of the tire.
Stiffness
3.04·105
Nm
Tire radial stiffness
C
Damping
500
N⁢sm
Tire radial damping
K
Size
R0
0.355
Unloaded tire radius
R_0
Surface
Flat surface
true
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
Time Lags
Use time lags
True (checked) means use time lags in slip calculation and enable the following two parameters
useTimeLag
Tlong
0.3
s
Time lag for longitudinal slip
Tlat
Time lag for slip angle
Visualization
Show tire
True (checked) creates a tire visualization and enables following three parameters
showTire
Dw
0.1
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
1
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
See Also
Tire Kinematics
Tires
Download Help Document