Engee documentation

Tire (Magic Formula)

The tire model determined by the coefficients of the empirical formula X.The boys.

blockType: Engee1DMechanical.Vehicles.Tires.MagicFormula

Path in the library:

/Physical Modeling/1D Mechanical/Tires & Vehicles/Tire (Magic Formula)

Description

Block Tire (Magic Formula) it represents the longitudinal movement of the tire, given by the empirical formula X.PAC [1], which is based on four coefficients. It is possible to simulate tire dynamics under constant or variable road surface conditions.

The longitudinal direction of the tire coincides with the direction of its movement on the road surface.

This block is a component based on the block Tire-Road Interaction (Magic Formula).

To improve the accuracy of the tire model, you can set properties such as ductility, inertia, rolling resistance, and a variable effective rolling radius. However, these properties increase the complexity of the bus model and can slow down the simulation. When calculating the model in real time or when preparing the model for semi-natural modeling, the pliability and inertia of the tire should be neglected.

Bus Model

In this block, a tire is considered as a rigid combination of a wheel and a tire that is in contact with the road and is prone to slipping. When the torque drives the wheel axle, the tire transmits the longitudinal force. , on the road. The tire transmits the resulting reaction in the form of force back to the wheel. This action rotates the wheel, creating a longitudinal motion. If we simulate the pliability of a tire, it will also deform elastically under load. If for the parameter Effective rolling radius model the value is set Load and velocity dependent (Magic Formula), then the radius of the tire will also change depending on the load and rotation speed.

The figure shows the forces acting on the tire.

tire magic formula 1

Variables that determine the bus model:

  • — rolling radius of the tire;

  • — the longitudinal speed of the wheel hub;

  • — longitudinal deformation of the tire;

  • — angular velocity of the wheel;

  • — the angular velocity of the point of contact. If there is no longitudinal deformation of the tire, then ;

  • — the longitudinal speed of the tire tread. Typically, the longitudinal tread velocity of a tire includes a component related to tire rotation , and an additional component related to tire deformation ;

  • — the sliding speed of the contact spot. If there is no longitudinal elastic deformation of the tire, then ;

  • — wheel sliding coefficient for tires without elastic deformation;

  • — the threshold speed of the wheel hub;

  • — the lower bound of the denominator of the slip coefficient;

  • — vertical tire load;

  • — the longitudinal force acting on the tire at the point of contact;

  • — longitudinal stiffness of the tire during deformation;

  • — longitudinal damping of the tire during deformation;

  • — the inertia of the wheel, such that the effective mass is ;

  • — the torque applied to the axle with the wheel.

Kinematics and reaction of the tire

You can simulate rolling, sliding and deformation.

Pumping and sliding_

The equation of translational motion of a non-sliding, non-deformable tire looks like this: . The sliding of the tires leads to a change in the longitudinal force .

The sliding speed of the contact spot is . For a tire without deformation .

The non-smoothed value of the slip coefficient of the contact spot is:

The block defines the denominator of the slip coefficient as:

where — parameter value Lower boundary of slip denominator, VXLOW.

Meaning smoothly changes to in transition areas:

,

.

The unit determines the sliding coefficient according to:

where

  • — parameter value Minimum valid wheel slip, KPUMIN;

  • — parameter value Maximum valid wheel slip, KPUMAX.

Meaning it changes smoothly in transition areas:

,

.

The block defines the smoothing threshold of the sliding coefficient as:

For this equation, the locked, sliding wheel has . For perfect rolling .

The formation

If the check box is selected Compliance, then the block considers the tire as elastic. When the tire is deformed, the point of contact of the tire with the road rotates with angular velocity slightly different from the speed of the wheel , which causes the contact patch to slip. The block defines a deformable tire as a translational spring damper with rigidity and damping .

If the checkbox is not checked Compliance Then , and there is no longitudinal deformation of the tire at any time in the simulation, then .

Dynamics of tires and wheels

This block consists of several sub-components. The equivalent block diagram is shown below.

tire magic formula

The block simulates the transient and stationary modes, as well as the start and stop. Blocks Translational Spring and Translational Damper are equivalent to rigidity and damping tires, respectively.

Block Tire-Road Interaction (Magic Formula) represents the longitudinal force on the bus as a function and using the empirical formula X .Guys, where — an independent sliding variable, and — input signal on port N.

Block Wheel and Axle — this is the rolling radius of the tire . The inertia value is the effective inertia. . Characteristic function of the bus determines the longitudinal force . Together with the torque of the drive shaft applied to the wheel axis, determines the angular and longitudinal movement of the wheel.

If the tire deformation is not modeled, then the block does not take into account the sub-components. Translational Spring and Translational Damper in an equivalent circuit, and the contact variables are returned to the wheel variables. In this case, the tire actually has infinite stiffness, and the P port of the block Wheel and Axle connects directly to the C port of the unit Force Sensor.

If in the block Rolling Resistance specified other than Pressure and velocity dependent (Magic Formula) rolling resistance model, then the equivalent circuit will not use the block Force Sensor and the T block port Rolling Resistance. If the checkbox is not checked Compliance, then the P port of the block Wheel and Axle connects directly to the C port of the unit Force Sensor or the T port of the unit Tire-Road Interaction (Magic Formula).

Effective rolling radius

This block determines the effective rolling radius, taking into account its increase or decrease due to centrifugal forces, so that

where

  • — effective rolling radius;

  • — increase in the free radius of the tire under the action of centrifugal forces;

  • — vertical tire load;

  • — parameter value Tire nominal vertical load, FNOMIN;

  • — parameter value Vertical stiffness;

  • — the speed of rotation of the wheel;

  • — parameter value Low load stiffness effective rolling radius, BREFF;

  • — parameter value Peak value of effective rolling radius, DREFF;

  • — parameter value High load stiffness effective rolling radius, FREFF;

  • — parameter value Ratio of nominal tire radius with non-rolling free tire radius, Q_RE0;

  • — parameter value Tire radius increase with speed, Q_V1.

Assumptions and limitations

  • The block assumes only longitudinal movement and does not take into account camber, turn or lateral movement.

  • Tire compliance implies a temporary delay in the tire’s response to forces acting on it. Time-delayed simulation increases the accuracy of the model, but increases the simulation time.

Ports

Input

# N — normal force, N
scalar

Details

Input port related to the normal force acting on the tyre, in N. The normal force is positive if it acts on the tyre in a downward direction, pressing it against the road surface.

Data types

Float64.

Complex numbers support

No

Conserving

# A — axis
`rotational mechanics

Details

A mechanical rotational port associated with an axis.

Program usage name

axle_flange

# H — hub
translational mechanics

Details

A mechanical progressive port connected to the hub of a wheel.

Program usage name

hub_flange

Output

# S — slip
scalar

Details

output port related to the slip coefficient, , between the tyre and the road.

Data types

Float64.

Complex numbers support

No

Parameters

Main

# Parameterize by — parameterization method
Peak longitudinal force and corresponding slip | Constant Magic Formula coefficients | Load-dependent Magic Formula coefficients | Physical signal Magic Formula coefficients

Details

Select how the block parameterises the bus using an empirical formula:

  • Peak longitudinal force and corresponding slip - parameterization of the empirical formula using the physical characteristics of the tyre.

  • Constant Magic Formula coefficients - set the parameters defining the constant coefficients , , and as scalars.

  • Load-dependent Magic Formula coefficients - set the parameters defining the load-dependent coefficients , , , , , and as vectors, one parameter for each coefficient.

  • Physical signal Magic Formula coefficients - set the coefficients of the empirical formula through the M port as a four-element vector .

Values

Peak longitudinal force and corresponding slip | Constant Magic Formula coefficients | Load-dependent Magic Formula coefficients | Physical signal Magic Formula coefficients

Default value

Peak longitudinal force and corresponding slip

Program usage name

friction_model

Evaluatable

No

# Rated vertical load — rated load force
N | nN | uN | mN | kN | MN | GN | dyn | lbf | kgf

Details

Rated vertical load force .

Dependencies

To use this parameter, set the Parameterize by parameters to . Peak longitudinal force and corresponding slip.

Units

N | nN | uN | mN | kN | MN | GN | dyn | lbf | kgf

Default value

3000.0 N

Program usage name

F_vertical_load

Evaluatable

Yes

# Peak longitudinal force at rated load — maximum longitudinal force at rated load
N | nN | uN | mN | kN | MN | GN | dyn | lbf | kgf

Details

Maximum longitudinal force , which the tyre exerts on the wheel when the vertical load is equal to its nominal value .

Dependencies

To use this parameter, set the parameters Parameterize by to . Peak longitudinal force and corresponding slip.

Units

N | nN | uN | mN | kN | MN | GN | dyn | lbf | kgf

Default value

3500.0 N

Program usage name

F_longitudinal_load

Evaluatable

Yes

# Slip at peak force at rated load (percent) — slip coefficient in per cent at maximum longitudinal force and rated load

Details

Slip coefficient , expressed as a percentage (%) when the longitudinal force is equal to the maximum value , and the vertical load is equal to the nominal value .

Dependencies

To use this parameter, set the parameters Parameterize by to Peak longitudinal force and corresponding slip.

Default value

10.0

Program usage name

percent_slip

Evaluatable

Yes

Geometry

# Rolling radius — tyre radius under load
m | um | mm | cm | km | in | ft | yd | mi | nmi

Details

Tyre radius under load .

Dependencies

To use this parameters, set:

  • for the parameters Parameterize by to be set to. Load-dependent Magic Formula coefficients`and for the parameters Effective rolling radius model to be set. `Constant radius;

  • for the parameters Parameterize by the value of. Peak longitudinal force and corresponding slip, Constant Magic Formula coefficients or Physical signal Magic Formula coefficients.

Units

m | um | mm | cm | km | in | ft | yd | mi | nmi

Default value

0.3 m

Program usage name

tire_radius

Evaluatable

Yes

Rolling Resistance

# Model rolling resistance — roll resistance

Details

Select this check box to take tyre rolling resistance into account in the simulation.

Default value

false (switched off)

Program usage name

enable_resistance

Evaluatable

No

Dynamics

# Compliance — consideration of tyre pliability

Details

Whether tyre pliability consideration should be included:

  • if the check box Compliance is not selected, the unit neglects dynamic deformation.

  • if the check box Compliance is selected, the unit treats the tyre as a spring-damped system under load.

Default value

false (switched off)

Program usage name

enable_compliance

Evaluatable

No

# Inertia — inertia model

Details

Select this check box to take the tyre inertia into account.

Default value

false (switched off)

Program usage name

enable_inertia

Evaluatable

No

Advanced

# Lower boundary of slip denominator, VXLOW — lower limit of the denominator of the slip coefficient
m/s | mm/s | cm/s | km/s | m/hr | km/hr | in/s | ft/s | mi/s | ft/min | mi/hr | kn

Details

Lower boundary of the slip factor denominator .

VXLOW - TIR file identifier.

Units

m/s | mm/s | cm/s | km/s | m/hr | km/hr | in/s | ft/s | mi/s | ft/min | mi/hr | kn

Default value

1.0 m/s

Program usage name

v_slip_threshold

Evaluatable

Yes

# Velocity threshold — tyre slip mode threshold
m/s | mm/s | cm/s | km/s | m/hr | km/hr | in/s | ft/s | mi/s | ft/min | mi/hr | kn

Details

Threshold speed , which the unit uses to transition between slip modes.

For details, see [Rolling and Sliding]. Rolling and Sliding for details.

Units

m/s | mm/s | cm/s | km/s | m/hr | km/hr | in/s | ft/s | mi/s | ft/min | mi/hr | kn

Default value

0.1 m/s

Program usage name

v_threshold

Evaluatable

Yes

# Minimum valid wheel slip, KPUMIN — minimum value of wheel slip coefficient

Details

The minimum permissible value of the wheel slip coefficient. A negative value means that the wheel slides in the opposite direction to the rotation.

Default value

-1.5

Program usage name

min_wheel_slip

Evaluatable

Yes

# Maximum valid wheel slip, KPUMAX — maximum value of wheel slip coefficient

Details

Maximum permissible value of the wheel slip coefficient.

Default value

1.5

Program usage name

max_wheel_slip

Evaluatable

Yes

Literature

  1. I.J.M. Besselink, A.J.C. Schmeitz, H. B. Pacejka, An Improved Magic Formula/Swift Tyre Model That Can Handle Inflation Pressure Changes, Vehicle System Dynamics 48, no. sup1 (December 2010): 337–52. https://doi.org/10.1080/00423111003748088.

  2. H.B. Pacejka., Tire and Vehicle Dynamics, Elsevier Science, 2005.