See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/315472151

The Simulation of Tire Dynamic Performance Based on qMagic Formulaq

Conference Paper · January 2017

DOI: 10.2991/amcce-17.2017.121

CITATIONS READS

3 138

4 authors, including:

Zhang Rongyun
Anhui University
38 PUBLICATIONS 135 CITATIONS

SEE PROFILE

All content following this page was uploaded by Zhang Rongyun on 11 September 2018.

The user has requested enhancement of the downloaded file.

 Advances in Engineering Research, volume 118
2nd International Conference on Automation, Mechanical Control and Computational Engineering (AMCCE 2017)

The Simulation of Tire Dynamic Performance Based on “Magic Formula”

Peicheng Shi, Qi Zhao, Rongyun Zhang and Li Ye

Anhui Polytechnic University, Anhui Engineering Technology Research Center of Automotive New
Technique, Wuhu 241000, China
[email protected], [email protected]

Keywords: Magic formula; Tire model; Matlab/Simulink; Dynamics

Abstract. Based on “Magic Formula”, the tire dynamics simulation model was established in the
Matlab/Simulink, contraposed many conditions like braking, steering and steering-braking
combination. The main load of tire and the wheel aligning torque were simulated, then got the
relationship curve between the longitudinal force/lateral force and slip rate and the relationship curve
between the wheel aligning torque and the side slip angle in corresponding conditions. The research
results show that the tire dynamic characteristic curve is in good agreement with the actual situation,
which not only verifies the correctness of the tire model, but also lays the foundation for the
subsequent vehicle dynamics simulation based on the tire model.

Introduction
Tire, the transfer components between the road and the vehicle, has a great impact on the vehicle’s
handling stability, and both vehicle motion and active safety control are achieved through the control
of tire force, so the simulation accuracy of the full vehicle’s dynamics characteristic is directly
influenced by the accuracy of tire model[1].
Theoretical model and empirical-semi-empirical model are the two main construction methods of
the tire models[2]. Generally speaking, theoretical model uses specific analytical expressions to
represent the mechanical properties of tires, but it get complex form, difficult calculation and poor
accuracy, which partly limit its application. Compared with the theoretical model, through the
analysis of abundant test data which reflect the tire strength characteristics. The empirical model or
semi-empirical model can effectively expresses the mechanical properties between ground and tire by
using a mathematical formula which includes fitting parameters. Empirical-semi-empirical model is
broader used due to its higher accuracy.
As one of the semi-empirical models, "Magic Formula" tire model is widely recognized as the
most outstanding tire-computing model in the simulation analysis of vehicle handling stability for its
good calculation accuracy [3]. For the "Magic Formula" model applied better during the analysis of
vehicle handling stability, a tire dynamic simulation model is established in Matlab / Simulink based
on "Magic Formula" model, and the simulation of braking, turning and turning braking combination
conditions are carried out, turning out the tire dynamic characteristic curves are in good agreement
with the actual situation.

Tire model
The concrete form of the “Magic Formula”shows in Eq. (1). On the basis of it, the force and torque of
the tire can be obtained by the fitting calculation from the experiment parameters. Each fitting
parameter is shown in Table 1.

Y = D sin{C arctan[Bx − E (Bx − arctan Bx )]}+ Sv (1)

Where Y are respectively the F x (lateral force), F y (longitudinal force), and M z (aligning torque), x
represents the side slip angle or longitudinal slip ratio according to different environments, D = the

Copyright © 2017, the Authors. Published by Atlantis Press.

This is an open access article under the CC BY-NC license (http://creativecommons.org/licenses/by-nc/4.0/). 699
 Advances in Engineering Research, volume 118

crest factor, B = the curve origin slope, E = the curve form factor, C = the shape characteristic factor,
S h , S v are respectively the Horizontal drift and the Vertical drift.
The “Magic Formula” uses combination formula of trigonometric function to fit the experimental
data of the tire, and it perfectly expresses the tire's longitudinal force F x , lateral force F y , aligning
torque M z , turning torque M x , resistance torque M y and the combined effect condition of longitudinal
force and lateral force[4,5] with the same formula. In this formula, the relation between input variables
and output variables of the tire is shown in Fig.1.

Fig. 1. Diagram of tire relationships

The longitudinal force of the steady-state braking environment

The slip ratio has a direct influence on the angular speed of the tire while the vehicle is moving. The
angular speed also changes under the same driving force or brake force[5]. With Eq. (1), the
relationship between tire longitudinal force and vertical load and slip ratio is shown as follows.

Fx 0 = Dx sin{Cx arctan[Bxλ − Ex (Bxλ − arctan Bxλ )]} (2)

2 2
Where C x = 1.65; D x = a 1 F z +a 2 F z , B x = B x C x D x /C x D x ; E x = a 6 F z +a 7 F z +a 8 , B x C x D x =
a 5 Fz
(a 3 F z 2+a 4 F z )× e ; a 1 ,··· ,a 8 are the fitting parameters shown in Table 1.

The lateral force of the steady-state steering environment

With Eq. (1), the relationships between lateral force F y0 and side slip angle and vertical load F z are
shown as follows.

 Fy 0 = Dy sin{Cy αrctαn[Byx − Ey (Byx − αrctαn Byx )]}+ Sv

 (3)
 x = α + Sh
Where C y = 1.3; D y = a 1 F z 2+a 2 F z ; B y C y D y = a 3 sin[a 4 arctan(a 5 F z )](1-a 12 |γ|); B y = B y C y D y /C y D y ;
E y = a 6 F z 2+a 7 F z +a 8 ; S h = a 9 γ; S v = (a 10 F z 2+a 11 F z )γ; γ = the vehicle camber angle; a 1 ,··· ,a 8 are fitting
parameters. Ignoring the effects of Sh and Sv, all fitting parameters are shown in Table 1.

Steering-braking in the environment of swerving with braking

According to 1.1 and 1.2, the tire force is solved in each single working condition, which turns out the
relations among different tire forces in the combined condition of turning and braking are as follows.

 σx
 Fx = σ Fx 0
 σy
(4)
 Fy = Fy 0
 σ
λ tαn α
Where σ = σx 2 + σy 2 ; σx = − ; σx = −
1+ λ 1+ λ

Aligning torque
On the basis of Eq.1, the mathematical model of tire aligning torque can be represented as follows.

700
 Advances in Engineering Research, volume 118

Mz = Dz sin{Cz αrctαn[Bzx − Ez (Bzx − αrctαn Bzx )]}+ Sv

 (5)
 x = α + Sh
a 5 Fz
Where C z = 2.4; D z = a 1 F z 2+a 2 F z ; B z C z D z = (a 3 F z 2+a 4 F z )×(1-a 12 |γ|) e ; B z = B z C z D z /C z D z ; S h
= a 9 γ; S v = (a 10 F z 2+a 11 F z )γ; E z = a 6 F z 2+a 7 F z +a 8 .

Table. 1. The fitting parameters

a1 a2 a3 a4 a5 a6 a7
Fx -21.3 1144 49.6 226 0.069 -0.006 0.056
Fy -22.1 1011 1078 1.82 0.208 0 -0.354
Mz -2.72 -2.28 -1.86 -2.73 0.11 -0.07 0.643
a8 a9 a 10 a 11 a 12 a 13
Fx 0.486 - - - - -
Fy 0.707 0.028 0 14.8 0.022 -
Mz -4.04 0.015 0.066 0.945 0.03 0.07

Tire dynamics model on the basis of Simulink

According to Eq.(2) to (5), the tire model and aligning torque model in the combined working
condition are built based on the MATLAB/Simulink , which is as shown in Fig.2 and Fig.3.

Fig. 2. Tire simulation model under the combined working condition

Fig. 3. Aligning torque computational model

Simulation analysis
In order to clearly show the relationship between the tire force / torque and slip rate, side slip angle,
the vertical loads are respectively set up as 1KN, 3KN, 5KN, 7KN to simulate the tire models in
various working conditions, and the curves got from simulation are as shown from Fig.4 to Fig.7.

701
 Advances in Engineering Research, volume 118

Fig. 4. Longitudinal force- Slip rate curve Fig. 5. Lateral force- Slip rate curve

Fig. 6. Aligning torque- Side slip angle curve Fig .7. Longitudinal force & Lateral force- Slip
rate curve

Fig.4 shows the response curves of the tire longitudinal force in pure braking working condition.
With the slip rate changes gradually, longitudinal force increases sharply. The longitudinal force
reaches the peak when slip rate is about 10%, and then it decreases gradually when it beyond 10%. It
concluded that the braking performance of the wheel is best when the slip rate is 10%. Fig.5 shows
the response curves of the lateral force in pure steering working condition. It is observed that the
lateral force increases linearly when the side slip angle in 0~5°, and then the lateral force changes
nonlinearly when the side slip rate is beyond 0~5°. Fig.6 shows the response curves of aligning
torque. It is visible the aligning torque reaches the peak when the side slip angle is about 5°, and with
the increasing of side slip angle, the aligning torque decreases gradually, and it tends to zero finally.
Generally, the influence of the aligning torque will not be taken into consideration during the analysis
of the vehicle dynamics. Fig.7 shows the relationship between the lateral force/ longitudinal force and
the side slip ratio in the combined working condition of turning and braking. It is visible with the
increasing of side slip ratio, the tire longitudinal force changes slightly, but the lateral force changes
greatly. From the above analysis, the tire model simulation of the force and torque changes have a
high degree of consistency with the actual situation, which verifies the correctness of the tire model.

Conclusion
The tire nonlinear characteristic has an important effect on the steering characteristic and the vehicle
riding stability. Therefore, it is necessary to study the tire dynamics. This paper carried on the
simulation modeling to tires based on “Magic Formula”. The models built in this paper can show a
variety of tire states changes such as load, side slip, roll and longitudinal slip rate, and it also can
effectively analyze the tire mechanics characteristics for any loads, any tire movement states.

702
 Advances in Engineering Research, volume 118

Through the analysis of the simulation results, it can be seen the tire dynamics model based on
"Magic Formula" exhibits a high fitting accuracy for the longitudinal force, the lateral force and the
aligning torque, which can fully reflect the vehicle tire dynamic characteristics. It also can be a
reference in future study.

Acknowledgements
The authors would like to thank Anhui province science and technology research key project
(Grant No. 1604a0902158)and Anhui university scientific research platform innovation team
building projects (2016-2018).

References
[1] Li Gang. Study on Stability and Energy Saving Control for X-by-wire Electric Vehicle with
Four-wheel Independent Drive In-wheel Motors [D]. Jilin University, (2013). (In Chinese)
[2] Hiroshi Fujimoto, Kiyoshi Fujii, Naoki Takahashi. Vehicle stability control of electric vehicle
with slip-ratio and cornering stiffness estimation [C]. Proceedings of the 2007 IEEE/ASME
International Conference on Advanced Intelligent Mechatronics. (2007) .
[3] Josko Deura, Danijel Pavkovic, Gilberto Burgio, et al. A model-based traction control strategy
non-reliant on wheel slip information [J]. Vehicle System Dynamics.(2011),49(8): 1245-1265.
[4] Li Songyan, Min Yongjun, Wang Liangmo, An Lihua. Establishment and Simulation Analysis of
Tire Dynamic Model [J]. Journal of Nanjing Institute of Technology (Natural Science Edition),
(2009), 07(03): 34-38. (In Chinese)
[5] Zheng Xiangmei, Gao Xingwang, Zhao Zhizhong. Simulation Analysis of Tire Dynamic Based
on “Magic Formula” [J]. Machinery & Electronics, (2012), (9): 16-20. (In Chinese)
[6] Liang Chunlan, Sheng Pengcheng, Liu Weize. MATLAB Simulation Research Based on the
“Magic Formula” Tire Model [J]. Journal of Xingtai Polytechnic College, (2013), 01:63-66. (In
Chinese)
[7] Yuan Zhongcheng, Lu Dang, Guo Konghui. Comparative Study on the UniTire Model and the
Magic Formula Model [J]. Automobile Technology, (2006),02: 7-11. (In Chinese)
[8] WANG Xianyun, LIU Yanhua, MA Bin, LIU Zhimin, LIU Liguo. Tire Characterizing Function
Computation Method Based on Magic Formula [J]. Science Technology and Engineering,
(2016),07: 265-269. (In Chinese)
[9] Wang Qian, Yang Zhijian, Ding Kang. Method in Identifying the Parameters of Magic Formula
Tire Model Based on New Self-adaptive Differential Evolution [J]. Journal of Mechanical
Engineering, (2014),06: 120-128. (In Chinese)
[10] Yin Li, Yu Jiangjiang. Simulation analysis research on tire dynamic characteristics [J].
Automobile Applied Technology, (2012), 09:1-4. (In Chinese)

703
View publication stats