# Nonlinear Research and Efficient Parameter Identification of Magic Formula Tire Model

Authors: Zhun Cheng; Zhixiong Lu
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

The Magic Formula tire model can describe the mechanical properties of tire accurately and thus is applied in the research field of vehicle dynamics widely. The Magic Formula tire model has the characteristics of a great number of parameters and the high nonlinearity, so it is hard to identify parameters. Researchers generally use different intelligent optimization algorithms for parameter identification. However, in the process of parameter identification, with a few experimental data, parameter identification results generally have the low accuracy, while, in the case of a large number of experimental data, the amount of work done in the experiment will increase and there will be many experimental errors. To solve these problems, this paper researches the longitudinal force of tire and proposes an interpolation method and a method based on the nonlinear research of the tire force. The results of parameter identification experiments on the two kinds of tire data show that both of the two methods can be used for the parameter identification of Magic Formula tire model fast and accurately with only a few experimental data. In addition, this paper proposes a method estimating the maximum longitudinal force and corresponding slip rate.


---

## Body

## 1. Introduction

Magic Formula tire model is able to well reflect the change rule of the tire force [1–4] and widely used; however, due to the nonlinear property and the required excessive parameters, there are many problems in the current researches on parameter identification of Magic Formula. van Oosten and Bakker [5] proposed that, in order to improve the fitting accuracy, 2000 and 1500 measured data points were needed, respectively, for pure slip the conditions and for combined slip conditions. Cabrera et al. [6, 7] and Ortiz et al. [8] proposed calculating the globally optimal solution based on differential evolution algorithm, but it is still likely to fall into the local minimum. Moreover, for the parameter identification of pure longitudinal force, the number of iterations for genetic algorithm is more than 500, and the elapsed time is 180 s. Wang et al. [9] proposed a new self-adaptive differential evolution algorithm (NSADE) and compared with other two differential evolution algorithms IOA and SSPDE; based on NSADE algorithm, the convergence to the global optimum is high in speed, the number of iterations is about 200, and the elapsed time is about 700 s. Zhang et al. [10] proposed a hybrid optimization method, based on which, firstly, the approximate optimal solution is solved according to the genetic algorithm, and then the accurate parameters are identified by the numerical optimization algorithm; the number of iterations for genetic algorithm is 10000; then the calculation results are well fit with the experimental data.In conclusion, there are some problems in the parameter identification of Magic Formula. If the measured data points are fewer, then the result of the parameter identification is not unique, and there will be many values [11]; due to a small amount of data from the introduced algorithm, many curves may meet the requirements of the objective function; however, for the same tire in the same case, the characteristic curve of the tire force is unique; moreover, the experimental error cannot show and guarantee that several measured data are the highest in accuracy, and the obtained curve is the actual tire force curve. If the measured points are increased, more experimental errors will be introduced, and the workload will be greatly increased. In addition, the interaction of parameter B,C,D,E increases the difficulty in algorithm identification; the researchers usually calculate B,C,D,E, respectively, and substitute the values into Magic Formula, so the number of the iterations for the algorithm increases, and a great number of data need to be calculated for iteration each time.In order to solve these problems, this paper proposes an interpolation method and a method based on the nonlinear research of the tire force. The former is based on 9 real experimental types of data under the same kind of load and uses the Lagrange, Hermite, and Spline interpolation to get more experimental data through simulation and then verifies the validity of experimental data obtained through simulation by calculating the errors which may be produced in simulated experimental data. In this way, we can get more experimental data with fewer experiments and then make parameter identification using the genetic algorithm. The method has the high accuracy and a very fast rate of convergence to the optimal solution. The latter, based on the nonlinear characteristic of Magic Formula tire model and the research on the changing laws for longitudinal force of tire, proposes two parameter identification models ofE>1 and E≤1, respectively, and a method estimating the maximum longitudinal force. The results indicate that under the same kind of load, we can make a fast and accurate parameter identification with a few, sometimes only three, real experimental types of data and get a very accurate estimate of maximum longitudinal force. This paper offers some reference to parameter identifications in other fields.


## 2. Parameter Identification of Magic Formula Tire Model Based on Interpolation Method

The basic Magic Formula of tire longitudinal force is as follows:(1)Yx=Dsin⁡Carctan⁡Bx-EBx-arctanBx,(2)yx=Yx+SV,(3)x=X+SH,where Y represents the theoretical longitudinal force of tire; y represents the real longitudinal force of tire; x represents the real slip rate; X represents the theoretical slip rate; parameters B,C,D,E,SH,SV represent the rigidity factor, the shape factor, the peak factor, the curvature factor, the horizontal displacement, and the vertical displacement, respectively. In addition, B,C,D,E,SH,SV are related to the function with the vertical load as the independent variable, so there are secondary parameters [12]. If we can get primary parameters accurately and then we can get secondary parameters easily, and because there are various tire models based on Magic Formula, this paper does not research solving secondary parameters. SH and SV are small, so the research in this paper neglects their values.The parameter identification method of tire model generally adopts the heuristic intelligent optimization algorithm, like the simulated annealing algorithm, the genetic algorithm, the particle swarm optimization algorithm, and so on. People generally use the coefficient of determinationRSQ as the objective function of iteration of algorithm. The closer the value of RSQ to 1 is, the higher the accuracy of parameter identification is. The calculation formula of RSQ is as follows:(4)RSQ=1-∑i=1nFx_predictedxi-Fx_measuredxi2∑i=1nFx_measuredxi-∑i=1nFx_measuredxi/n2,where Fx_predicted represents the predicted value of longitudinal force of tire; Fx_measured represents the measured value of longitudinal force of tire; xi represents the ith slip rate; n represents the total number of measured values.The number of experimental data of tire force is the key to parameter identification of Magic Formula tire model. If the number of measured data points is limited, many curves may meet the requirements of objective function, and some curves making the objective function optimal may have strange shapes and thus make the results of parameter identification untrustworthy. Although the tire force has small errors in measurement points, it has considerable errors in some unmeasured points locally; that is, there may be many abnormal points. Substantial data points measured may increase the difficulty and time consumption of experiment and also introduce more experimental errors. With the interpolation method, we can obtain the high-accuracy simulation of experimental data and the proper interpolation method can reduce the error of interpolated data points and the errors of first and second derivatives and then ensure the accuracy of interpolation.To further ensure the accuracy of interpolation, this paper uses the piecewise interpolation method and each segment is based on 4 measured data points. This is because, according to the changing laws [13], the changes of longitudinal force of tire are generally divided into the following three segments: the first segment approximates to a line; the second segment approximates to a parabola with the slope decreasing first and then increasing; the third segment approximates to a curve with the slope decreasing continuously. Therefore, the parameter identification based on the interpolation method only requires the measurements of 9 types of data under each vertical load.This paper chooses the measured data of 2 kinds of tire for the research. The first kind of tire’s measured data comes from literature [14]; the second kind of tire’s measured data comes from ADAMS built-in mdi_pac_94 tire data. The numbers of measured data of two kinds of tire are both 101 under each vertical load. According to the previous analysis, we select 9 measured types of data under the same vertical load as the raw data for interpolation method. Figure 1 shows all measured data and selected 9 measured data types of the first kind of tire. Figure 2 shows all measured data of the second kind of tire.Figure 1
The comparison of simulated experimental data and real experimental data of the first kind of tire obtained with the Lagrange interpolation.Figure 2
The comparison of experimental value and predicted value of the second kind of tire.


### 2.1. Simulating Experimental Data Based on Lagrange Interpolation

According to the previous analysis, we divide the 9 measured data types selected into the following three segments by slip rate: segment one[0,2,6,8], segment two [8,12,15,20], and segment three [20,50,80,100]. Then, we make the interpolation processing for each segment with the n-order Lagrange interpolation polynomial [15–18]. After the processing, each segment has 97 tire force data types obtained through simulation and 4 measured data types. Figures 1 and 2 show the simulated data of tire longitudinal force obtained based on Lagrange interpolation.Figure1 shows that the simulated data of interpolation are very close to real experimental data, with the maximum relative error of 6.39% and the average relative error of 0.499% which is very small.Then, we use Lagrange interpolation remainderRx to calculate and simulate the maximum error [16–18] which may be reached by experimental data. The estimation formula of error Rx of Lagrange interpolation is as follows:(5)Rx=fx-Pnx≤Mn+1n+1!ωnx,where x represents simulated slip rate data; it is any insertion point in the interval [x0,xn];xi(i=0,1,…,n) is the real slip rate; f(x) is the Magic Formula; Pn(x) is the Lagrange interpolation function; Mn+1=maxx0≤x≤xn⁡fn+1x;ωn(x)=(x-x0)(x-x1)⋯(x-xn), andn’s value is 3.Table1 shows the maximum errors may be reached by simulated experimental data obtained with Lagrange interpolation. In the table, Fz represents the vertical load of tire, and slip represents the slip rate.Table 1
The maximum errors of simulated experimental data by Lagrange interpolation method.
F
z (kN)
Slip (%)
[
0,2
,
6,8
]
[
8,12,15,20
]
[
20
,
50,80,100
]
2
127.5880 N
46.3727 N
204.7502 N
4
452.0734 N
94.3467 N
271.1843 N
6
1097.03 N
129.7200 N
264.4633 NFigure2 shows the results of the data of the second kind of tire with the vertical load of 2000 N in the same processing mode.Figure2 shows that, for the data of the second kind of tire, measured data tally with the experimental data obtained through simulation. The maximum absolute errors may be reached by the simulated data in the first, second, and third segments which are 0.1682 N, 0.0958 N, and 168.9035 N, respectively. The data’s maximum relative error with real experimental data is 1.08% and average relative error is 0.13% which is very small.


### 2.2. Simulating Experimental Data Based on Hermite Interpolation

In practical problems, the Hermite interpolation [16–18] can not only reduce the error between the simulated value and the real value but also reduce the first-order derivative error and thus makes the curve of interpolation have the first-order smoothness. The estimation formula of error Rx of cubic Hermite interpolation is as follows:(6)Rx=fx-H3x≤h4384M4,where f(x) is the Magic Formula, H3(x) is the cubic Hermite interpolation function, h=max0≤i≤n-1⁡xi+1-xi,M4=maxx0≤x≤xn⁡f4x, andn’s value is 3.We make the Hermite interpolation processing for the data of the first kind of tire and calculate the maximum error which may be reached by the experimental data obtained through simulation. Table2 shows the results.Table 2
The maximum errors of simulated experimental data by Hermite interpolation method.
F
z (kN)
Slip (%)
[
0,2
,
6,8
]
[
8,12,15,20
]
[
20
,
50,80,100
]
2
26.5117 N
1.9490 N
91.3059 N
4
106.4405 N
4.0898 N
117.1631 N
6
224.9895 N
5.8826 N
142.3747 NComparing with the real measured value of tire, the simulated experimental data of the first kind of tire shows the maximum relative error of 3.69% and the average relative error of 0.82% which is small.The second kind of tire’s maximum errors may be reached in the first, second, and third segments which are 0.0048 N, 0.0056 N, and 12.7693 N, respectively. Comparing with real experimental data, the data’s maximum relative error is 2.05% and average relative error is 0.18% which is small.


### 2.3. Simulating Experimental Data Based on Spline Interpolation

The Spline interpolation [16–18] can meet the higher requirement on smoothness relative to the Hermite interpolation. The method has first-order and second-order continuous differential quotients in the joining point. The estimation formula of error Rx of cubic Spline interpolation is as follows:(7)Rx=fx-S3x≤eM4h4,where fx is Magic Formula, S3x is cubic Spline interpolation function, e=5/384,h=max0≤i≤n-1⁡xi+1-xi,M4=maxx0≤x≤xn⁡f(4)(x), andn’s value is 3.This paper makes the Spline interpolation processing for the data of the first kind of tire and calculates the maximum errors which may be reached by experimental data obtained through simulation, and Table3 shows the results.Table 3
The maximum errors of simulated experimental data by Spline interpolation method.
F
z (kN)
Slip (%)
[
0,2
,
6,8
]
[
8,12,15,20
]
[
20
,
50,80,100
]
2
138.8472 N
10.0133 N
406.1219 N
4
591.0820 N
20.4681 N
528.9091 N
6
1343.9231 N
28.6259 N
644.3526 NComparing with the real measured value of tire, the simulated experimental data of the first kind of tire has the maximum relative error of 5.78% and the average relative error of 0.55% which is small.The second kind of tire’s maximum errors may be reached in the first, second, and third segments which are 0.0197 N, 0.0238 N, and 60.3164 N, respectively. Comparing with real experimental data, the data’s maximum relative error is 1.08% and average relative error is 0.13% which is small.


### 2.4. Parameter Identification of Magic Formula Tire Model Based on Genetic Algorithm

According to the content above, the data of tire force obtained based on Lagrange, Hermite, and Spline interpolation almost agree with real data of tire force. Then, we choose the tire force data obtained based on the Lagrange interpolation and use the genetic algorithm [19–21] and consider RSQ as the objective function of algorithm iteration to identify the parameters in Magic Formula tire model. The genetic algorithm’s initial population is generated randomly, and the number of individuals NIND=200 and the hereditary algebra GEN=30. Table 4 gives the results of primary parameter identification of the first kind of tire.Table 4
Primary parameter identification results of the first kind of tire.
Vertical loadFZ (kN)
B
C
D
E
R
S
Q
2
0.1549
1.64
2202.97
0.5643
0.9998
4
0.1836
1.64
4238.65
0.5997
0.9997
6
0.2080
1.64
6099.19
0.5708
0.9996Figure3 shows that the algorithm can get the optimal solution in the 15th generation or so, and the computer consumes the time of about 10 s. Because simulated experimental data are close to real experimental data, the error of any point in the curve can reach the minimum, thus eliminating the possibility of outlier locally. To better compare the results of parameter identification, we compare all measured data (101 data types under the same kind of vertical load) of the first kind of tire with the estimated values of Magic Formula tire mode which have been identified. As shown in Figure 4, the error is very small.Figure 3
Evolution graph of parameter identification of the first kind of tire with genetic algorithm.Figure 4
Comparison of experimental values and predicted values of the first kind of tire.The values of primary parameter identification of the second kind of tire areB=0.0161,C=1.3635,D=3085.20,E=1.4189,RSQ=0.999995.


## 2.1. Simulating Experimental Data Based on Lagrange Interpolation

According to the previous analysis, we divide the 9 measured data types selected into the following three segments by slip rate: segment one[0,2,6,8], segment two [8,12,15,20], and segment three [20,50,80,100]. Then, we make the interpolation processing for each segment with the n-order Lagrange interpolation polynomial [15–18]. After the processing, each segment has 97 tire force data types obtained through simulation and 4 measured data types. Figures 1 and 2 show the simulated data of tire longitudinal force obtained based on Lagrange interpolation.Figure1 shows that the simulated data of interpolation are very close to real experimental data, with the maximum relative error of 6.39% and the average relative error of 0.499% which is very small.Then, we use Lagrange interpolation remainderRx to calculate and simulate the maximum error [16–18] which may be reached by experimental data. The estimation formula of error Rx of Lagrange interpolation is as follows:(5)Rx=fx-Pnx≤Mn+1n+1!ωnx,where x represents simulated slip rate data; it is any insertion point in the interval [x0,xn];xi(i=0,1,…,n) is the real slip rate; f(x) is the Magic Formula; Pn(x) is the Lagrange interpolation function; Mn+1=maxx0≤x≤xn⁡fn+1x;ωn(x)=(x-x0)(x-x1)⋯(x-xn), andn’s value is 3.Table1 shows the maximum errors may be reached by simulated experimental data obtained with Lagrange interpolation. In the table, Fz represents the vertical load of tire, and slip represents the slip rate.Table 1
The maximum errors of simulated experimental data by Lagrange interpolation method.
F
z (kN)
Slip (%)
[
0,2
,
6,8
]
[
8,12,15,20
]
[
20
,
50,80,100
]
2
127.5880 N
46.3727 N
204.7502 N
4
452.0734 N
94.3467 N
271.1843 N
6
1097.03 N
129.7200 N
264.4633 NFigure2 shows the results of the data of the second kind of tire with the vertical load of 2000 N in the same processing mode.Figure2 shows that, for the data of the second kind of tire, measured data tally with the experimental data obtained through simulation. The maximum absolute errors may be reached by the simulated data in the first, second, and third segments which are 0.1682 N, 0.0958 N, and 168.9035 N, respectively. The data’s maximum relative error with real experimental data is 1.08% and average relative error is 0.13% which is very small.


## 2.2. Simulating Experimental Data Based on Hermite Interpolation

In practical problems, the Hermite interpolation [16–18] can not only reduce the error between the simulated value and the real value but also reduce the first-order derivative error and thus makes the curve of interpolation have the first-order smoothness. The estimation formula of error Rx of cubic Hermite interpolation is as follows:(6)Rx=fx-H3x≤h4384M4,where f(x) is the Magic Formula, H3(x) is the cubic Hermite interpolation function, h=max0≤i≤n-1⁡xi+1-xi,M4=maxx0≤x≤xn⁡f4x, andn’s value is 3.We make the Hermite interpolation processing for the data of the first kind of tire and calculate the maximum error which may be reached by the experimental data obtained through simulation. Table2 shows the results.Table 2
The maximum errors of simulated experimental data by Hermite interpolation method.
F
z (kN)
Slip (%)
[
0,2
,
6,8
]
[
8,12,15,20
]
[
20
,
50,80,100
]
2
26.5117 N
1.9490 N
91.3059 N
4
106.4405 N
4.0898 N
117.1631 N
6
224.9895 N
5.8826 N
142.3747 NComparing with the real measured value of tire, the simulated experimental data of the first kind of tire shows the maximum relative error of 3.69% and the average relative error of 0.82% which is small.The second kind of tire’s maximum errors may be reached in the first, second, and third segments which are 0.0048 N, 0.0056 N, and 12.7693 N, respectively. Comparing with real experimental data, the data’s maximum relative error is 2.05% and average relative error is 0.18% which is small.


## 2.3. Simulating Experimental Data Based on Spline Interpolation

The Spline interpolation [16–18] can meet the higher requirement on smoothness relative to the Hermite interpolation. The method has first-order and second-order continuous differential quotients in the joining point. The estimation formula of error Rx of cubic Spline interpolation is as follows:(7)Rx=fx-S3x≤eM4h4,where fx is Magic Formula, S3x is cubic Spline interpolation function, e=5/384,h=max0≤i≤n-1⁡xi+1-xi,M4=maxx0≤x≤xn⁡f(4)(x), andn’s value is 3.This paper makes the Spline interpolation processing for the data of the first kind of tire and calculates the maximum errors which may be reached by experimental data obtained through simulation, and Table3 shows the results.Table 3
The maximum errors of simulated experimental data by Spline interpolation method.
F
z (kN)
Slip (%)
[
0,2
,
6,8
]
[
8,12,15,20
]
[
20
,
50,80,100
]
2
138.8472 N
10.0133 N
406.1219 N
4
591.0820 N
20.4681 N
528.9091 N
6
1343.9231 N
28.6259 N
644.3526 NComparing with the real measured value of tire, the simulated experimental data of the first kind of tire has the maximum relative error of 5.78% and the average relative error of 0.55% which is small.The second kind of tire’s maximum errors may be reached in the first, second, and third segments which are 0.0197 N, 0.0238 N, and 60.3164 N, respectively. Comparing with real experimental data, the data’s maximum relative error is 1.08% and average relative error is 0.13% which is small.


## 2.4. Parameter Identification of Magic Formula Tire Model Based on Genetic Algorithm

According to the content above, the data of tire force obtained based on Lagrange, Hermite, and Spline interpolation almost agree with real data of tire force. Then, we choose the tire force data obtained based on the Lagrange interpolation and use the genetic algorithm [19–21] and consider RSQ as the objective function of algorithm iteration to identify the parameters in Magic Formula tire model. The genetic algorithm’s initial population is generated randomly, and the number of individuals NIND=200 and the hereditary algebra GEN=30. Table 4 gives the results of primary parameter identification of the first kind of tire.Table 4
Primary parameter identification results of the first kind of tire.
Vertical loadFZ (kN)
B
C
D
E
R
S
Q
2
0.1549
1.64
2202.97
0.5643
0.9998
4
0.1836
1.64
4238.65
0.5997
0.9997
6
0.2080
1.64
6099.19
0.5708
0.9996Figure3 shows that the algorithm can get the optimal solution in the 15th generation or so, and the computer consumes the time of about 10 s. Because simulated experimental data are close to real experimental data, the error of any point in the curve can reach the minimum, thus eliminating the possibility of outlier locally. To better compare the results of parameter identification, we compare all measured data (101 data types under the same kind of vertical load) of the first kind of tire with the estimated values of Magic Formula tire mode which have been identified. As shown in Figure 4, the error is very small.Figure 3
Evolution graph of parameter identification of the first kind of tire with genetic algorithm.Figure 4
Comparison of experimental values and predicted values of the first kind of tire.The values of primary parameter identification of the second kind of tire areB=0.0161,C=1.3635,D=3085.20,E=1.4189,RSQ=0.999995.


## 3. Parameter Identification of Magic Formula Tire Model Based on the Nonlinear Research of the Tire Force

### 3.1. Magic Formula Nonlinear Research

Set(8)F=arctan⁡Bx-EBx-arctan⁡Bx.IfC·F=±π/2, then, in formula (1), yMax=±D.Find the derivative of (1), and attain (9)y˙x=Dcos⁡C·F·C·F˙.Sety˙(x)=0; then cos⁡C·F=0 or F˙=0. Hence, only when C·F=±π/2 or F˙=0, y˙(x) has the extreme point.Setx=0, y˙x=0=BCD. D is the peak factor; when the slip rate x>0, assuming D is positive, the slope at the origin is significantly positive, so B has the same positive and negative characteristics as C. (10)f=Bx-EBx-arctan⁡Bx=1-EBx+Earctan⁡Bx.Find the derivative of (10), and attain(11)f˙=1-EB+E·BB2x2+1.Find the derivative of (8), and attain (12)F˙=1f2+1·f˙.If and only ifx=+∞, F˙=0, will the tire longitudinal force generate the extreme point within the scope of the slip rate [0,100]; then F˙=0; hence, only when C·F=±π/2, the extreme point must be the maximum point.IfB>0, then C>0; according to (11), f is monotonically increasing; if B<0, then C<0, and f is monotonically decreasing. f is the odd function, so the composite function F is also the odd function; moreover, F and f have the same monotonicity.B
,
C must be the same symbol; (8) and (10) are the even functions with respect to B, so (1) is the even function with respect to B,C, namely, yB>0,C>0=yB<0,C<0. Hence, only B>0 and C>0 will be discussed below.


#### 3.1.1. IfE≤1

f
˙
>
0, F˙>0, F is monotonically increasing, and x→+∞, F→π/2. If C·F=π/2, x∈0,100, then C>1, yMax=D. Assuming the slip rate corresponding to the maximal longitudinal force as xm, it is easy to attain (13)E=Bxm-tan⁡π/2CBxm-arctan⁡Bxm,B=y˙x=0CD=y˙x=0C·yMax.ForE≤1, as long as C and y˙x=0 are solved, xm and y˙x=0 values are measured; then the parameters of Magic Formula can be solved rapidly. Assume k=arctan⁡y˙x=0=arctan⁡BCD, due to the character of arctan function; when the independent variable is greater than 1, the scope of the dependent variable change is greatly narrowed, so that the independent variable ranges from 40 to 5000; the dependent variable, namely, the angle, is only increased by about 1°. Compared with a large number of tire data, the slope range of tangent line for the longitudinal force Magic Formula at the origin is [40,5729.6]; the corresponding angle range is [88.57°,89.99°], namely, k∈1.5458,1.5706.Combined with the simulated experimental data based on interpolation method mentioned above, as long as the size of longitudinal force from the experimental measured data conforms to the “low-high-low” situation,xm and yx=xm values can be precisely calculated. In the process of identification, select all k values by enumeration method, so that only C is parameter to be calculated, and other parameter values are attained by relation. Theoretically, under each load, as long as three experimental data types are met, the supersimplified parameter identification can be done, because the three data can clearly reflect all information of the tire longitudinal force, such as the values of zero-order, first-order, and second-order derivatives; the experimental data from interpolation simulation can well reflect these rules.If the measured experimental data cannot meet the “low-high-low” requirement or the value span of slip rate is larger, allxm and k values are selected by enumeration method; the parameters to be calculated are B and C; other parameter values can be attained by relation.


#### 3.1.2. IfE>1

Whenx=0, then f˙Max=B, and subsequently f˙ is monotonically decreasing; when x→+∞, then f˙Max→1-EB<0; hence, f˙ is positive firstly and then negative.Hence,F is firstly increased to the maximum and decreased; three cases are discussed respectively; when x∈0,100, discuss the relationship between C·F and π/2. In Case I, if C·FMax=π/2, then the solution process of the parameters is similar to that of E≤1; in Case II, if C·FMax>π/2, then there will two maximums, and significantly they do not conform to the actual situation of the tire force, so this case is eliminated; in Case III, if C·FMax<π/2, then yMax<D. Assume f˙=0, xm=1/E-1·1/B. The calculated parameter values are as follows: (14)E=1+1Bxm2,D=y˙x=0BC.Therefore, only B and C are the parameters to be calculated; for the calculation of other parameter values refer to the analysis of E≤1.


### 3.2. Supersimplified Parameter Identification of Magic Formula

In combination with the study of Magic Formula, based onB,C,D,E nonlinear relationship with respect to the tire force, the supersimplified parameter identification is done. Based on the data for the first and second kind of tire, only three measured data types are selected under each vertical load; the tire force of the measured data needs to meet “low-high-low” requirement; in order to correct the model and select any data point with larger slip rate, a total of four data types are measured, and the parameter identification is based on the genetic algorithm.The curve close to longitudinal force peak within small scope is approximate to parabola; so Hermite interpolation method is used to predictxm and yx=xm.Table5 shows the comparison between the predicted values and measured values of xm and yx=xm.Table 5
Comparison between the predicted values and measured values ofxmand yx=xm.
F
z (kN)
Predictedxm (%)
Predictedyx=xm (N)
Measuredxm (%)
Measuredyx=xm (N)
x
m error (%)
y
x
=
x
m error (%)
2
12.02
2202.75
11.9
2202.8
1.01
0.002
4
10.02
4234.44
10.3
4235.2
2.72
0.018
6
9.98
6090.57
9.2
6097.2
8
0.11For the first kind of tire data, the calculation formulas ofE>1 and E≤1 are used, respectively, for parameter identification; when E>1, there is no solution, so E≤1. k∈1.5458,1.5706, select 0.0001 as the precision to select all k values within the area, corresponding to the vertical load Fz=2,4,6 (kN) of the first kind of tire; the optimal value of k is 1.5690, 1.5700, and 1.5703. C is the sole parameter to be calculated by genetic algorithm; as for all the measured data, a total of 101 measured data, all RSQ values are greater than 0.9999; moreover, the optimizing time is about 5 s; the parameter identification results are shown in Table 6.Table 6
Identification result of the primary parameters for the first kind of tire.
Vertical loadFz (kN)
B
C
D
E
R
S
Q
2
0.1539
1.64
2202.80
0.5394
0.9999988
4
0.1808
1.64
4235.20
0.5669
0.9999893
6
0.2015
1.64
6097.20
0.5509
0.9999446For the second kind of tire, the calculation formulas ofE>1 and E≤1 are used, respectively, for parameter identification; when E≤1, there is no solution. When E>1, the parameter identification result is k=1.5566,B=0.0173,C=1.2547,D=3197.05,E=1.3545,RSQ=0.9987. xm predicted value is 96.99, the measured value is 97.0, and the error is 0.01%. yx=xm predicted value is 2402.42, the measured value is 2402.42, and the error is 0%.The identification results of two tires are consistent with that mentioned above by interpolation method; furthermore, RSQ values of all the measured data are higher.


## 3.1. Magic Formula Nonlinear Research

Set(8)F=arctan⁡Bx-EBx-arctan⁡Bx.IfC·F=±π/2, then, in formula (1), yMax=±D.Find the derivative of (1), and attain (9)y˙x=Dcos⁡C·F·C·F˙.Sety˙(x)=0; then cos⁡C·F=0 or F˙=0. Hence, only when C·F=±π/2 or F˙=0, y˙(x) has the extreme point.Setx=0, y˙x=0=BCD. D is the peak factor; when the slip rate x>0, assuming D is positive, the slope at the origin is significantly positive, so B has the same positive and negative characteristics as C. (10)f=Bx-EBx-arctan⁡Bx=1-EBx+Earctan⁡Bx.Find the derivative of (10), and attain(11)f˙=1-EB+E·BB2x2+1.Find the derivative of (8), and attain (12)F˙=1f2+1·f˙.If and only ifx=+∞, F˙=0, will the tire longitudinal force generate the extreme point within the scope of the slip rate [0,100]; then F˙=0; hence, only when C·F=±π/2, the extreme point must be the maximum point.IfB>0, then C>0; according to (11), f is monotonically increasing; if B<0, then C<0, and f is monotonically decreasing. f is the odd function, so the composite function F is also the odd function; moreover, F and f have the same monotonicity.B
,
C must be the same symbol; (8) and (10) are the even functions with respect to B, so (1) is the even function with respect to B,C, namely, yB>0,C>0=yB<0,C<0. Hence, only B>0 and C>0 will be discussed below.


### 3.1.1. IfE≤1

f
˙
>
0, F˙>0, F is monotonically increasing, and x→+∞, F→π/2. If C·F=π/2, x∈0,100, then C>1, yMax=D. Assuming the slip rate corresponding to the maximal longitudinal force as xm, it is easy to attain (13)E=Bxm-tan⁡π/2CBxm-arctan⁡Bxm,B=y˙x=0CD=y˙x=0C·yMax.ForE≤1, as long as C and y˙x=0 are solved, xm and y˙x=0 values are measured; then the parameters of Magic Formula can be solved rapidly. Assume k=arctan⁡y˙x=0=arctan⁡BCD, due to the character of arctan function; when the independent variable is greater than 1, the scope of the dependent variable change is greatly narrowed, so that the independent variable ranges from 40 to 5000; the dependent variable, namely, the angle, is only increased by about 1°. Compared with a large number of tire data, the slope range of tangent line for the longitudinal force Magic Formula at the origin is [40,5729.6]; the corresponding angle range is [88.57°,89.99°], namely, k∈1.5458,1.5706.Combined with the simulated experimental data based on interpolation method mentioned above, as long as the size of longitudinal force from the experimental measured data conforms to the “low-high-low” situation,xm and yx=xm values can be precisely calculated. In the process of identification, select all k values by enumeration method, so that only C is parameter to be calculated, and other parameter values are attained by relation. Theoretically, under each load, as long as three experimental data types are met, the supersimplified parameter identification can be done, because the three data can clearly reflect all information of the tire longitudinal force, such as the values of zero-order, first-order, and second-order derivatives; the experimental data from interpolation simulation can well reflect these rules.If the measured experimental data cannot meet the “low-high-low” requirement or the value span of slip rate is larger, allxm and k values are selected by enumeration method; the parameters to be calculated are B and C; other parameter values can be attained by relation.


### 3.1.2. IfE>1

Whenx=0, then f˙Max=B, and subsequently f˙ is monotonically decreasing; when x→+∞, then f˙Max→1-EB<0; hence, f˙ is positive firstly and then negative.Hence,F is firstly increased to the maximum and decreased; three cases are discussed respectively; when x∈0,100, discuss the relationship between C·F and π/2. In Case I, if C·FMax=π/2, then the solution process of the parameters is similar to that of E≤1; in Case II, if C·FMax>π/2, then there will two maximums, and significantly they do not conform to the actual situation of the tire force, so this case is eliminated; in Case III, if C·FMax<π/2, then yMax<D. Assume f˙=0, xm=1/E-1·1/B. The calculated parameter values are as follows: (14)E=1+1Bxm2,D=y˙x=0BC.Therefore, only B and C are the parameters to be calculated; for the calculation of other parameter values refer to the analysis of E≤1.


## 3.1.1. IfE≤1

f
˙
>
0, F˙>0, F is monotonically increasing, and x→+∞, F→π/2. If C·F=π/2, x∈0,100, then C>1, yMax=D. Assuming the slip rate corresponding to the maximal longitudinal force as xm, it is easy to attain (13)E=Bxm-tan⁡π/2CBxm-arctan⁡Bxm,B=y˙x=0CD=y˙x=0C·yMax.ForE≤1, as long as C and y˙x=0 are solved, xm and y˙x=0 values are measured; then the parameters of Magic Formula can be solved rapidly. Assume k=arctan⁡y˙x=0=arctan⁡BCD, due to the character of arctan function; when the independent variable is greater than 1, the scope of the dependent variable change is greatly narrowed, so that the independent variable ranges from 40 to 5000; the dependent variable, namely, the angle, is only increased by about 1°. Compared with a large number of tire data, the slope range of tangent line for the longitudinal force Magic Formula at the origin is [40,5729.6]; the corresponding angle range is [88.57°,89.99°], namely, k∈1.5458,1.5706.Combined with the simulated experimental data based on interpolation method mentioned above, as long as the size of longitudinal force from the experimental measured data conforms to the “low-high-low” situation,xm and yx=xm values can be precisely calculated. In the process of identification, select all k values by enumeration method, so that only C is parameter to be calculated, and other parameter values are attained by relation. Theoretically, under each load, as long as three experimental data types are met, the supersimplified parameter identification can be done, because the three data can clearly reflect all information of the tire longitudinal force, such as the values of zero-order, first-order, and second-order derivatives; the experimental data from interpolation simulation can well reflect these rules.If the measured experimental data cannot meet the “low-high-low” requirement or the value span of slip rate is larger, allxm and k values are selected by enumeration method; the parameters to be calculated are B and C; other parameter values can be attained by relation.


## 3.1.2. IfE>1

Whenx=0, then f˙Max=B, and subsequently f˙ is monotonically decreasing; when x→+∞, then f˙Max→1-EB<0; hence, f˙ is positive firstly and then negative.Hence,F is firstly increased to the maximum and decreased; three cases are discussed respectively; when x∈0,100, discuss the relationship between C·F and π/2. In Case I, if C·FMax=π/2, then the solution process of the parameters is similar to that of E≤1; in Case II, if C·FMax>π/2, then there will two maximums, and significantly they do not conform to the actual situation of the tire force, so this case is eliminated; in Case III, if C·FMax<π/2, then yMax<D. Assume f˙=0, xm=1/E-1·1/B. The calculated parameter values are as follows: (14)E=1+1Bxm2,D=y˙x=0BC.Therefore, only B and C are the parameters to be calculated; for the calculation of other parameter values refer to the analysis of E≤1.


## 3.2. Supersimplified Parameter Identification of Magic Formula

In combination with the study of Magic Formula, based onB,C,D,E nonlinear relationship with respect to the tire force, the supersimplified parameter identification is done. Based on the data for the first and second kind of tire, only three measured data types are selected under each vertical load; the tire force of the measured data needs to meet “low-high-low” requirement; in order to correct the model and select any data point with larger slip rate, a total of four data types are measured, and the parameter identification is based on the genetic algorithm.The curve close to longitudinal force peak within small scope is approximate to parabola; so Hermite interpolation method is used to predictxm and yx=xm.Table5 shows the comparison between the predicted values and measured values of xm and yx=xm.Table 5
Comparison between the predicted values and measured values ofxmand yx=xm.
F
z (kN)
Predictedxm (%)
Predictedyx=xm (N)
Measuredxm (%)
Measuredyx=xm (N)
x
m error (%)
y
x
=
x
m error (%)
2
12.02
2202.75
11.9
2202.8
1.01
0.002
4
10.02
4234.44
10.3
4235.2
2.72
0.018
6
9.98
6090.57
9.2
6097.2
8
0.11For the first kind of tire data, the calculation formulas ofE>1 and E≤1 are used, respectively, for parameter identification; when E>1, there is no solution, so E≤1. k∈1.5458,1.5706, select 0.0001 as the precision to select all k values within the area, corresponding to the vertical load Fz=2,4,6 (kN) of the first kind of tire; the optimal value of k is 1.5690, 1.5700, and 1.5703. C is the sole parameter to be calculated by genetic algorithm; as for all the measured data, a total of 101 measured data, all RSQ values are greater than 0.9999; moreover, the optimizing time is about 5 s; the parameter identification results are shown in Table 6.Table 6
Identification result of the primary parameters for the first kind of tire.
Vertical loadFz (kN)
B
C
D
E
R
S
Q
2
0.1539
1.64
2202.80
0.5394
0.9999988
4
0.1808
1.64
4235.20
0.5669
0.9999893
6
0.2015
1.64
6097.20
0.5509
0.9999446For the second kind of tire, the calculation formulas ofE>1 and E≤1 are used, respectively, for parameter identification; when E≤1, there is no solution. When E>1, the parameter identification result is k=1.5566,B=0.0173,C=1.2547,D=3197.05,E=1.3545,RSQ=0.9987. xm predicted value is 96.99, the measured value is 97.0, and the error is 0.01%. yx=xm predicted value is 2402.42, the measured value is 2402.42, and the error is 0%.The identification results of two tires are consistent with that mentioned above by interpolation method; furthermore, RSQ values of all the measured data are higher.


## 4. The Maximal Slip Rate, Maximal Tire Longitudinal Force Prediction, and Parameter Identification

If three groups of measured data with smaller spacing cannot meet “low-high-low” requirement, it is feasible to predictk and xm value and then identify the parameters. k∈1.5458,1.5706, select 0.0001 as the precision to select all k values within the area; in each section of k value, make xm vary from 0 to 100, and substitute it into the derivation formula of the parameter identification, and still select RSQ with respect to all the measured data as the objective function.WhenFz=2000 N, xm measured value is 11.9. Compared with the predicted value in Table 7, it can be seen that the slip rate value is closer to xm measured value; the prediction effect will be better. In Table 7, all the xm measured values are almost less than 20; then compared with the measured value near to the slip rate 20, xm value can be obtained accurately. The selection method of k value is in analogy with xm acquisition method; it can be obtained that k value is 1.5690 and equal to the optimal value of k mentioned in Table 7.Table 7
Predicted results of the first kind of tire under the vertical loadFz=2000 N.
Slip rate value (%)
x
m (%)
k
10,12,14
12
1.5690
20,22,24
12
1.5690
30,32,34
12
1.5690
40,42,44
12
1.5690
50,52,54
16
1.5681
60,62,64
17
1.5678
70,72,74
26
1.5653
80,82,84
14
1.5685
90,92,94
22
1.5664Table8 shows the predicted results of the second kind of tire under the vertical load Fz=2000 N.Table 8
Predicted results of the second kind of tire under the vertical loadFz=2000 N.
Slip rate value (%)
x
m (%)
k
10,12,14
90
1.5561
20,22,24
90
1.5561
30,32,34
92
1.5561
40,42,44
95
1.5563
50,52,54
96
1.5563
60,62,64
97
1.5562
70,72,74
97
1.5562
80,82,84
97
1.5565
90,92,94
97
1.5565According to Table8, xm measured value is close to 90; find the measured value close to the slip rate 90; then xm predicted value is 97, k=1.5565, which is almost equal to the measured value 97 and k value 1.5566 mentioned in Table 8.Afterxm and k values are obtained, combine with the measured data for parameter identification, and the process is as mentioned in Section 3.If only the slip rates=[0,10,20,30,40,60,80,100] is selected for the first kind of tire data, for convenience of study, only calculate the tire longitudinal force Fx = [0, 2191.82, 2121.42, 1983.58, 1869.11, 1707.55, 1603.00, 1530.81] when Fz=2000 N. By calculation, xm=12, k=1.5690, which are consistent with the previously measured value. In order to compare with the identification result, use the improved simulated annealing algorithm [22–24] for parameter identification.Table9 shows the results of two parameter identification methods.Table 9
Results of two parameter identification methods.
Identification method
B
C
D
E
Number of iterations
R
S
Q (8 measured values)
R
S
Q
Simulated annealing algorithm
0.4466
1.7303
2267.18
0.9094
364
0.9993
0.6970
Method proposed in this paper
0.1517
1.6624
2207.50
0.5811
3
0.9982
0.9999Figure5 shows the comparison between the predicted values and measured values by new method and simulated annealing algorithm.Figure 5
Comparison between the predicted values and measured values by new method and simulated annealing algorithm.It can be seen from Table9 that the slip rate spacing of the measured points is larger and the quantity is less, if the intelligent algorithm, such as simulated annealing algorithm, is used to guarantee the optimal objective function as the premise, namely, only guarantee the error of the measured points to be minimal; however, for the position beyond the measured points, the error is larger, and RSQ value with respect to a total of 101 measured points is only 0.6970; but the slip rate is less than 10%; RSQ value is only 0.4565. Furthermore, if the identification programs are executed for many times, many groups of parameter solutions will be obtained, and almost these solutions only meet the requirement of 8 measured values; but, in terms of all other measured values, the errors are larger. However, the new method proposed in this paper can accurately calculate the maximal longitudinal force and the corresponding slip rate; moreover, the number of iterations is few; the predicted values by identification are almost equal to all the measured values.


## 5. Conclusion

This paper proposes an interpolation method and a method based on the nonlinear research of the tire force for the parameter identification of Magic Formula tire model. The research shows that the Lagrange, Hermite, and Spline interpolation methods used in this paper can all simulate the experimental data of longitudinal force of tire accurately. According to the changing laws for tire longitudinal force of tire, we divide experimental data into three segments so that we can improve the accuracy of tire force data obtained with interpolation effectively. According to analysis on errors of different interpolation methods, Lagrange, Hermite, and Spline interpolation methods have different ranges of application and are suitable for longitudinal force of tire with simulated slip rates of 20%~100%, 0~20%, and 8%~20%, respectively. Research results indicate that with only 9 experimental data types under the same kind of vertical load, we can ensure the high accuracy of longitudinal force of tire obtained with interpolation. Based on the tire longitudinal force data obtained with interpolation, we can make the parameter identification of Magic Formula tire model with the genetic algorithm fast and accurately. In the parameter identification process, the number of iterations is 15 and the time consumed by computer is about 10 s. With the method, we not only solve the problem of low accuracy of parameter identification caused by the limited number of experimental data but also reduce the workload in the experiment.This paper makes the nonlinear research on Magic Formula, deduces the interactive relations ofB,C,D, and E in Magic Formula, and proposes two kinds of parameter identification models with E>1 anE ≤ 1, respectively, and a method predicting the maximum longitudinal force. Research results show that with only 3 experimental data types under the same vertical load, we can make the parameter identification fast and accurately. The method reduces the number of iterations of parameter identification significantly to only 3~15 iterations which take 3~5 seconds to produce the result on a desktop with Quad CPU 3.4 GHz, and the maximum longitudinal force predicted coincides with real experimental data.



---

*Source: 6924506-2017-09-28.xml*