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.
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.
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.
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.
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.
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.
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.
#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
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.
H.B. Pacejka., Tire and Vehicle Dynamics, Elsevier Science, 2005.