Tire-road dynamics given by magic formula coefficients

Simscape / Driveline / Tires & Vehicles / Tire Subcomponents

This block models the interaction between the tire tread and road pavement. The longitudinal force arising from this interaction is given by the magic formula, an empirical equation based on four fitting coefficients. Tire properties such as compliance and inertia are ignored.

Port T represents the tire tread. Port H represents the wheel hub that transmits the thrust generated by the tire to the remainder of the vehicle. Port N accepts a physical signal input of the normal force acting on the tire. The normal force is positive if it acts downward on the tire, pressing it against the pavement. Port S outputs a physical signal with the tire slip measured during simulation.

The block provides two friction variants. The default variant, ```
Fixed
friction coefficients
```

, accepts the magic formula coefficients
as block parameters. This variant treats the coefficients as constants
or load-dependent parameters. Use this variant to model tire dynamics
under constant pavement conditions.

The alternative variant, `Variable friction coefficients`

,
accepts the magic formula coefficients as a physical signal input.
Use this variant to model tire dynamics under variable pavement conditions.
Selecting this variant exposes physical signal port M. Use this port
to provide the magic formula coefficients as a four-element vector,
specified in the order [*b*, *c*, *d*, *e*].

To change block variants:

Right-click the tire block.

In the context-sensitive menu, select

**Simscape**>**Block choice**.Select the desired block variant.

The Tire-Road Interaction (Magic Formula) block models the longitudinal forces at the tire-road contact patch using the Magic Formula of Pacejka [7].

The figure displays the forces on the tire. The table defines the tire model variables.

**Tire-Road Contact Variables**

Symbol | Description and Unit |
---|---|

Ω | Wheel angular velocity |

r_{w} | Wheel radius |

V_{x} | Wheel hub longitudinal velocity |

r_{w}Ω | Tire tread longitudinal velocity |

V_{sx} = r_{w}Ω
– V_{x} | Wheel slip velocity = tread velocity T – hub velocity H |

κ = V_{sx}/|V_{x}| | Wheel slip |

F_{z}, F_{z0} | Vertical load and nominal vertical load on tire |

F_{x} = f(κ, F_{z}) | Longitudinal force exerted on the tire at the contact point. Also a characteristic function f of the
tire. |

A tire model provides a steady-state *tire characteristic
function* *F*_{x} = *f*(*κ*, *F*_{z}),
the longitudinal force *F*_{x} on
the tire, based on:

Vertical load

*F*_{z}Wheel slip

*κ*

The Magic Formula is a specific form for the tire characteristic
function, characterized by four dimensionless coefficients (*B*, *C*, *D*, *E*),
or stiffness, shape, peak, and curvature:

*F*_{x} = *f*(*κ*, *F*_{z})
= *F*_{z}·*D*·sin( *C*·arctan[
{ *B**κ* – *E*·[ *B**κ* –
arctan(*B**κ*) ] } ] ) .

The slope of *f* at *κ* =
0 is *B**C**D*·*F*_{z}.

A more general Magic Formula uses dimensionless coefficients that are functions of the tire load. A more complex set of parameters p_i, entered in the dialog box, specifies these functions:

*F*_{x0} = *D*_{x}·sin( *C*_{x}·arctan[
{ *B*_{x}*κ*_{x} – *E*_{x}·[ *B*_{x}*κ*_{x} –
arctan(*B*_{x}*κ*_{x})
] } ] ) + *S*_{Vx} ,

where

*df*_{z} = (*F*_{z} – *F*_{z0})/*F*_{z} ,

*κ*_{x} = *κ* +
S_{Hx} ,

*C*_{x} = p_Cx1
,

*D*_{x} = *μ*_{x}·*F*_{z} ,

*μ*_{x} =
p_Dx1 + p_Dx2·*df*_{z} ,

*E*_{x} = (p_Ex1
+ p_Ex2·*df*_{z} + p_Ex3·*df*_{z}^{2})[1
– p_Ex4·sgn(*κ*_{x})]
,

*K*_{xκ} = *F*_{z}·(p_Kx1
+ p_Kx2·*df*_{z})·exp(p_Kx3·*df*_{z})
,

*B*_{x} = *K*_{xκ}/(*C*_{x}*D*_{x} + *ε*_{x})
,

*S*_{Hx} = p_Hx1
+ p_Hx2·*df*_{z} ,

*S*_{Vx} = *F*_{z}·(p_Vx1
+ p_Vx2·*df*_{z}) .

*S*_{Hx} and *S*_{Vx} represent
offsets to the slip and longitudinal force in the force-slip function,
or horizontal and vertical offsets if the function is plotted as a
curve. *μ*_{x} is the longitudinal
load-dependent friction coefficient. *ε*_{x} is
a small number inserted to prevent division by zero as *F*_{z} approaches
zero.

The block uses a representative set of Magic Formula coefficients.
The block scales the coefficients to yield the peak longitudinal force *F*_{x0} at
the corresponding slip *κ*_{0} that
you specify, for rated vertical load *F*_{z0}.

Numerical values are based on empirical tire data. These values are typical sets of constant Magic Formula coefficients for common road conditions.

Surface | B | C | D | E |
---|---|---|---|---|

Dry tarmac | 10 | 1.9 | 1 | 0.97 |

Wet tarmac | 12 | 2.3 | 0.82 | 1 |

Snow | 5 | 2 | 0.3 | 1 |

Ice | 4 | 2 | 0.1 | 1 |

The Tire-Road Interaction (Magic Formula) block assumes longitudinal motion only and includes no camber, turning, or lateral motion.

Port | Description |
---|---|

A | Conserving mechanical rotational port associated with the axle of the tire |

H | Conserving mechanical translational port associated with the hub of the tire |

N | Physical signal input port associated with the normal force on the tire |

M | Physical signal input associated with the static and kinetic friction coefficients of the tire |

S | Physical signal output port associated with the relative slip between the tire and road |

**Parameterize by**Select how to use the Magic Formula to model the tire-road interaction. The default setting is

`Peak longitudinal force and corresponding slip`

. This and all dependent parameters are visible only when the block variant is set to`Fixed friction coefficients`

. For more information about the block variants, see the block description.`Peak longitudinal force and corresponding slip`

— Parametrize the Magic Formula with physical characteristics of the tire.`Constant Magic Formula coefficients`

— Parametrize the Magic Formula directly with constant coefficients. If you select this option, the panel changes from its default.`Load-dependent Magic Formula coefficients`

— Parametrize the Magic Formula directly with load-dependent coefficients. If you select this option, the panel changes from its default.

**Velocity threshold**The wheel hub velocity

*V*_{th}below which the slip calculation is modified to avoid singular evolution at zero velocity. Must be positive. The default is`0.1`

.From the drop-down list, choose units. The default is meters per second (

`m/s`

).

Was this topic helpful?